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Abstract 

Accurate  tropical  cyclone  (TC)  intensity  estimates  are  best  achieved  from  satellite 
observations.  The  Advanced  Microwave  Sounding  Unit  (AMSU)  has  operated  since  1998  on 
polar-orbiting  environmental  satellites  and  is  able  to  measure  the  warm  temperature  anomaly  in 
the  upper  troposphere  above  a  TC's  center.  Through  hydrostatic  equilibrium,  this  warm  anomaly 
is  roughly  proportional  to  the  TC's  sea-level  pressure  anomaly.  Based  on  this  principle,  the 
Cooperative  Institute  for  Meteorological  Satellite  Studies  (CEMSS)  provides  near  real-time 
AMSU-based  estimates  of  TC  minimum  sea-level  pressure  (MSLP)  to  forecast  centers 
worldwide.  These  estimates  are  as  accurate  as  the  benchmark  Dvorak  technique,  but  are  subject 
to  error  caused  by  precipitation  effects  (primarily  brightness  temperature  reduction  by  scattering) 
on  the  AMSU  55  GHz  channels  sensitive  to  upper-tropospheric  temperature. 

Simulated  AMSU  brightness  temperatures  (TB's)  are  produced  by  a  polarized 
reverse  Monte  Carlo  radiative  transfer  model  using  representative  TC  precipitation  profiles. 
Results  suggest  that  precipitation  depression  of  high-frequency  window  channel  TB's  is 
correlated  with  depression  of  sounding  channel  TB's  and  can  be  used  to  correct  for  scattering 
effects  on  the  AMSU  channels  used  in  TC  intensity  estimates.  Analysis  of  AMSU  data  over  the 
tropical  oceans  confirms  this,  and  forms  the  basis  for  an  empirical  scattering  correction  using 
AMSU  31  and  89  GHz  TB's.  This  scattering  correction  reduces  CIMSS  TC  MSLP  algorithm 
RMS  error  by  10%  in  a  7-year,  497  observation  sample. 
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1.  Introduction 

Tropical  cyclones  are  among  the  most  deadly  and  damaging  of  natural  phenomena.  In 
extreme  cases,  such  as  the  Bangladesh  Typhoon  of  1970,  they  have  killed  300,000  people 
(WMO  Tropical  Cyclone  Programme  2000).  In  the  Western  world,  timely  warnings,  mass 
communication,  and  efficient  transportation  make  loss  of  life  much  smaller.  However,  rapid 
development  along  vulnerable  coastlines  has  caused  a  different  problem:  exponentially 
increasing  damage  costs.  Hurricane  Andrew,  the  costliest  natural  disaster  to  strike  the  U.S., 
caused  $25  billion  of  damage  to  south  Florida  in  1992  (Rappaport  1993).  U.S.  hurricane  damage 
in  a  typical  year  approaches  $5  billion  (Pielke  and  Landsea  1998).  Evacuation  costs  alone  can  be 
staggering,  often  estimated  at  near  $1  million  per  mile  of  affected  coastline  (Whitehead  2000). 

With  a  large  number  of  bases  and  deployed  forces  in  the  tropical  Western  Pacific  and 
Indian  Ocean,  as  well  as  the  Southeast  United  States,  the  Department  of  Defense  (DoD)  is 
extremely  vulnerable  as  well  to  fatalities,  damage,  and  disrupted  operations  caused  by  TC's.  To 
minimize  the  risk  of  damage  to  critical  ships,  aircraft,  and  hardware,  commanders  are  forced  to 
evacuate  them  from  installations  in  the  path  of  a  TC,  costing  several  hundred  thousand  dollars— 
and  several  days  of  lost  training-per  occurrence  per  installation  (Dumas  and  Tibbetts  2004). 

Because  of  these  significant  impacts,  affected  nations  maintain  TC  forecast  and  warning 
centers.  The  U.S.  National  Hurricane  Center  (NHC)  and  Central  Pacific  Hurricane  Center  join 
counterparts  from  Japan,  Australia,  Fiji,  India,  and  France  to  form  the  World  Meteorological 
Organization's  network  of  Regional  Specialized  Meteorological  Centers  (RSMC's),  each  having 
official  TC  warning  responsibility  for  a  portion  of  the  tropical  oceans.  DoD  also  maintains  the 
Joint  Typhoon  Warning  Center  (JTWC)  at  Pearl  Harbor,  Hawaii,  with  responsibility  for  timely 
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Forecast  Intensity  Error  +  24hr  Forecast  Intensity  Error  +  24hr 


Figure  1.1:  Comparison  of  TC  track  (top  row)  and  intensity  (bottom  row)  forecast  accuracy 
trends.  Joint  Typhoon  Warning  Center  Western  North  Pacific  basin  performance  is  in  the  left 
column;  National  Hurricane  Center  Atlantic  basin  is  in  the  right  column;  24-,  48-,  and  72- 
hour  forecast  accuracies  are  depicted  (from  http://www.nhc.noaa.gov/verification/  and 
http://www.npmoc.navy.mil/jtwc/climostats/Statclimo.html). 


TC  forecast  and  warning  support  to  U.S.  Government  assets  throughout  the  Pacific  and  Indian 
Oceans  (USCINCPAC  Inst  3 140. IX). 

Accurately  predicting  TC  effects  requires  correctly  forecasting  storm  track,  maximum 
sustained  wind  speed  ( Vmax ),  and  spatial  wind  distribution.  Numerical  weather  prediction  (NWP) 
model  improvements  and  focused  research  focus  have  produced  impressive  recent  gains  in  track 
forecast  performance.  But  as  shown  in  Figure  1.1.,  intensity  (in  the  form  of  Vmax)  forecast 
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performance  has  improved  little  over  the  past  decade.  Accordingly,  research  emphasis  has 
shifted  somewhat  in  recent  years,  increasing  effort  devoted  to  accurate  TC  intensity  and  wind 
field  structure  forecasts. 

TC  warning  centers  utilize  a  suite  of  NWP  models,  statistical  aids,  and  empirical  rules 
when  formulating  their  official  TC  track  and  intensity  forecasts.  At  present,  the  best  intensity 
forecast  tools  are  statistical,  using  multiple  regression  to  predict  future  intensity  from  a  selection 
of  current  TC  parameters  (DeMaria  and  Kaplan  1999;  Knaff  et  al.  2003).  For  intensity  forecasts 
out  to  72  hours,  the  leading  predictors  used  by  the  statistical  models  are  Vmax  and  its  12-hour 
trend.  Accurate  intensity  forecasts,  therefore,  depend  critically  upon  accurate  current  intensity 
observations. 

Tropical  cyclone  winds  are  driven  by  the  strong  horizontal  pressure  gradient  between  the 
environment  and  the  low  pressure  at  the  storm's  center,  so  TC  intensity  can  be  described  by 
either  Vmax  or  the  minimum  sea-level  pressure  {MSLP)  at  storm  center.  Neither  measure  is 
entirely  satisfactory,  though.  A  single  maximum  wind  speed  value  is  not  representative  of  the 
highly  asymmetric  3-d  wind  distribution  caused  by  storm  motion,  shear,  eddy-scale  variability, 
and  variations  in  storm  size.  While  MSLP  gives  a  first-order  estimate  of  the  maximum  storm 
wind,  it  does  not  describe  the  fine-scale  variability  in  the  horizontal  pressure  gradient,  which  is 
responsible  for  the  maximum  sustained  wind.  Several  attempts  have  been  made  to  rectify  these 
definitions  of  intensity  through  development  of  empirical  wind-pressure  relationships  or  standard 
radial  distributions  of  wind  and  pressure  (Atkinson  and  Holliday  1977;  Holland  1980).  To  first- 
order,  these  relationships  provide  a  conversion  between  and  MSLP  as  measures  of  intensity. 
For  the  purposes  of  this  study,  MSLP  will  be  used  for  intensity,  since  it  is  more  readily  observed 
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by  aircraft  reconnaissance  than  the  highly-localized  Vma*. 

Forecasters  and  researchers  use  three  primary  means  to  measure  MSLP.  The  most 
desirable  is  in  situ  measurement  via  aircraft  reconnaissance.  Since  the  end  of  World  War  II,  the 
U.S.  Air  Force,  U.S.  Navy,  and  National  Oceanographic  and  Atmospheric  Administration 
(NOAA)  have  performed  this  mission  in  the  Western  North  Pacific  (WNP),  Central  Pacific  (CP), 
Eastern  Pacific  (EP),  and  Atlantic  basins  (Weatherford  and  Gray  1988;  Martin  and  Gray  1993). 
Aircraft,  such  as  the  WC-130H/J  and  WP-3D  currently  in  use,  have  the  endurance  and  speed  to 
reach  storms  up  to  1000  miles  out  to  sea.  Once  there,  these  aircraft  can  loiter  long  enough  to 
locate  the  storm's  center  and  penetrate  it  multiple  times,  measuring  wind  and  temperature  aloft 
and  deploying  dropsondes  to  collect  vertical  profiles  of  wind,  temperature,  pressure,  and 
humidity  (Hock  and  Franklin  1999).  Unfortunately,  reconnaissance  operations  are  currently 
limited  to  the  western  part  of  the  Atlantic  basin,  with  occasional  CP  and  EP  missions. 

In  the  absence  of  reconnaissance,  the  only  in  situ  wind  and  pressure  observations  come 
from  manned  or  automated  weather  stations,  ships,  and  buoys.  Given  the  vast  size  of  the  tropical 
ocean  basins,  the  relatively  small  size  of  a  TC,  and  instrument  uncertainty  (or  damage)  caused  by 
the  extremely  high  winds  and  rough  seas,  the  number  of  usable  in  situ  surface  observations  of 
TC  intensity  is  extremely  small. 

With  the  relative  sparcity  of  in  situ  surface  and  aircraft  reconnaissance  observations,  TC 
forecasters  have  come  to  rely  most  heavily  on  the  third  source  of  TC  observations:  weather 
satellites.  Shortly  after  the  introduction  of  operational  geostationary  weather  satellites,  Dvorak 
(1973;  1975)  developed  an  empirical  technique  for  estimating  intensity  using  TC  features  in 
visible  satellite  imagery.  This  technique,  expanded  to  incorporate  features  observed  in  enhanced 
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infrared  (IR)  imagery  (Dvorak  1984),  has  become  the  standard  for  operational  TC  intensity 
estimates.  It  employs  pattern  matching,  identification  of  maximum  and  minimum  cloud 
temperatures,  and  measurements  of  the  degree  of  spiral  band  curvature,  along  with  intensity 
change  rules  and  constraints  to  produce  a  current  intensity  on  a  numeric  scale  ranging  from  0.5 
to  8.0.  Analysts  at  the  RSMC's,  JTWC,  the  NOAA  National  Environmental  Satellite  Data  and 
Information  Service  Satellite  Analysis  Branch,  and  the  Air  Force  Weather  Agency  use  the 
Dvorak  technique  to  produce  operational  MSLP  estimates  with  root  mean  square  error  (RMSE) 
of  approximately  10  hPa  (Olander  et  al.  2004).  An  automated  version  of  the  Dvorak  technique 
developed  by  the  University  of  Wisconsin-Madison's  Cooperative  Institute  for  Meteorological 
Satellite  Studies  (CIMSS)  achieves  similar  results  (Olander  et  al.  2004).  The  excellent 
performance  of  both  the  manual  and  automated  Dvorak  techniques,  combined  with  the  excellent 
update  frequency  (hourly  or  half-hourly)  of  geostationary  imagery,  justifies  their  widespread 
operational  use.  However,  the  Dvorak  technique  has  shortcomings.  The  manual  technique  is 
labor-intensive,  requiring  a  highly  trained  and  experienced  satellite  analyst.  Its  rules,  combined 
with  the  ambiguities  of  imagery  interpretation,  can  introduce  considerable  subjectivity.  Dvorak's 
constraints  on  intensity  change  rates  cause  estimates  to  lag  rapid  TC  intensity  changes  (Velden  et 
al.  1998).  Finally,  empirical  techniques  such  as  this  are  conditioned  to  mean  events  and  often 
perform  poorly  in  extreme  cases.  A  more  physical  and  more  objective  remotely  sensed  TC 
intensity  estimation  technique  would  therefore  complement  the  Dvorak  techniques  in  operational 
use. 

Passive  and  active  microwave  instruments  on  polar-orbiting  weather  satellites  offer  some 
alternatives.  The  thick  cirrus  cloud  that  often  canopies  the  center  of  TC's  is  nearly  transparent  at 
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microwave  frequencies  (3  to  300  GHz).  Rain  emits  strongly  at  the  lower  end  of  this  frequency 
band,  while  frozen  precipitation  scatters  strongly  at  the  high  end.  This  makes  the  convective 
clouds  in  TC  spiral  bands  and  eyewalls  very  distinct  in  imagery  from  passive  microwave  imagers 
like  the  Special  Sensor  Microwave  Imager  (SSMI),  Tropical  Rainfall  Measuring  Mission 
(TRMM)  Microwave  Imager  (TMI),  and  Advanced  Microwave  Scanning  Radiometer  (AMSR). 
This  enables  highly  accurate  TC  position  fixes  in  cases  where  IR  imagery  would  have  difficulties 
(Alliss  et  al.  1993).  Passive  microwave  imagers  can  also  infer  sea  surface  wind  speed  from  the 
increasing  microwave  emissivity  of  the  wind-roughened  ocean  surface.  Active  microwave 
scatterometers  like  the  SeaWinds  instrument  are  able  to  do  the  same  via  ocean-surface  radar 
reflectivity,  as  well  as  infer  wind  direction  by  observing  surface  reflectivity  from  different  angles 
(Katsaros  et  al.  2001;  Yueh  et  al.  2003).  Both  types  of  instrument  can  give  a  large-scale  picture 
of  TC  wind  field  structure,  but  their  relatively  coarse  spatial  resolution  and  constraints  on  the 
maximum  wind  measurable  by  these  techniques  prevent  either  active  radar  or  passive  imagers 
from  accurately  estimating  Vmax,  and  neither  type  is  capable  of  estimating  MSLP. 

Polar-orbiting  passive  microwave  temperature  sounding  instruments,  however,  are  well- 
suited  to  estimating  MSLP.  Hydrostatic  equilibrium  relates  surface  pressure  to  the  temperature 
of  the  column  of  air  above.  From  the  temperature  difference  between  columns  in  the  ambient 
environment  and  at  the  TC  center,  we  can  infer  the  MSLP  difference  from  the  ambient  sea-level 
pressure.  Kidder  et  al.  (1978)  first  demonstrated  this  principle  using  an  experimental  microwave 
temperature  sounder.  At  CIMSS,  Velden  and  Smith  (1983),  Velden  (1989),  and  Velden  et  al. 
(1991)  expanded  on  the  idea  using  the  Microwave  Sounding  Unit  (MSU),  the  first-generation 
operational  instrument.  Today,  both  CIMSS  and  the  Cooperative  Institute  for  Research  in  the 
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Atmosphere  (CIRA)  produce  real-time  TC  intensity  estimates  using  the  Advanced  Microwave 
Sounding  Unit  (AMSU)  (Brueske  and  Velden  2003;  DeMuth  et  al.  2004).  This  type  of  intensity 
estimate  is  completely  automated  and  objective  and  can  achieve  accuracies  superior  to 
operational  Dvorak  intensity  estimates  (Herndon  et  al.  2004).  While  microwave  sounding 
instrument  TC  intensity  estimation  techniques  are  promising,  they  have  their  own  shortcomings: 
The  AMSU  instrument  is  carried  on  NOAA  Polar-orbiting  Operational  Environmental  Satellite 
(POES)  spacecraft.  There  are  currently  3  operational  POES  spacecraft,  which  limits  data  refresh 
to  4-6  hours,  compared  to  an  hour  or  less  for  IR  imagery  from  geostationary  satellites.  While 
geostationary  satellites  transmit  instrument  data  to  ground  processing  centers  nearly 
instantaneously,  polar-orbiting  spacecraft  store  data  for  playback  to  ground  stations  once  per 
orbit.  Frequently,  though,  polar-orbiters  may  go  2  or  3  orbits  between  playbacks,  making  data 
unavailable  for  up  to  five  hours  after  it  was  collected.  (The  National  Polar-orbiting  Operational 
Environmental  Satellite  System,  NPOESS,  plans  to  improve  this  situation  by  using  a  distributed 
worldwide  network  of  ground  stations  when  it  becomes  operational  after  2010)  The  coarse 
instrument  resolution  of  microwave  temperature  sounding  instruments  (AMSU-A  ranges  from  48 
to  150  km)  causes  it  to  sub-sample  the  relatively  small  TC  warm  core.  The  effect  is  currently 
addressed  at  CIMSS  using  ancillary  estimates  of  TC  size  (Herndon  et  al.  2004).  And  finally, 
large  liquid  and  frozen  hydrometeors  in  the  concentrated  cumulonimbus  convection  near  storm 
center  can  scatter  upwelling  microwave  radiation,  decreasing  the  accuracy  with  which 
temperature  sounding  instruments  can  measure  the  TC's  warm  core  magnitude  and  estimate  its 
MSLP.  The  purpose  of  this  study  is  to  foster  better  TC  intensity  forecasts  by  improving  the 
accuracy  of  CIMSS'  microwave  TC  intensity  estimates  through  quantifying  and  correcting  for 
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these  precipitation  effects. 

Chapter  2  of  this  dissertation  introduces  background  information  on  the  tropical  cyclone 
warm  core,  passive  microwave  temperature  sounding,  existing  microwave  TC  intensity 
estimation  techniques,  and  the  technical  details  of  the  AMSU-A  and  -B  instruments.  Chapter  3 
discusses  how  precipitation  affects  microwave  radiative  transfer  and  introduces  a  conceptual 
framework  for  estimating  precipitation  effects  on  AMSU  sounding  channels.  Chapter  4  presents 
results  from  numerical  simulations  of  precipitation  effects  on  AMSU  sounding  channels.  Details 
of  the  Reverse  Monte  Carlo  radiative  transfer  model  developed  for  this  purpose  are  contained  in 
the  Appendix.  Chapter  5  analyzes  AMSU  observations  of  tropical  convective  precipitation  and 
discusses  an  empirical  method  for  correcting  its  effect  on  temperature  sounding  channels. 
Chapter  6  applies  this  correction  to  the  CIMSS  AMSU  TC  intensity  estimation  technique  and 
validates  its  performance.  Chapter  7  presents  concluding  discussion  and  outlines  a  future 
application. 
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2.  Background 

Before  Chapter  3  discusses  precipitation  effects  on  microwave  radiative  transfer,  Chapter 
2  summarizes  four  key  areas  of  background.  Section  2.1  discusses  the  tropical  cyclone  (TC) 
warm  core  and  its  relation  to  the  TC's  minimum  sea-level  pressure  ( MSLP ).  Section  2.2  outlines 
microwave  radiative  transfer  and  temperature  sounding  principles.  Section  2.3  reviews  previous 
efforts  to  remotely  senseTC  MSLP  or  maximum  sustained  wind  (Vmax)  via  the  TC's  warm  core 
using  microwave  sounding  instruments.  Finally,  Section  2.4  presents  the  detailed  characteristics 
of  the  Advanced  Microwave  Sounding  Unit  (AMSU). 


2.1.  The  Tropical  Cyclone  Warm  Core 

TC's  are  long-lived  oceanic  warm-core  cyclonic  circulations  that  extend  vertically 
through  the  full  depth  of  the  troposphere  and  radially  for  hundreds  of  km.  Their  energy  source  is 
latent  and  sensible  heat  transferred  from  ocean  to  atmosphere  by  evaporation  and  conduction, 
then  transported  to  the  upper  troposphere  by  vertical  motion  in  organized  deep  convection 
(Emanuel  2003;  Simpson  et  al.  1997).  The  weakest  TC's,  with  sustained  winds  of  less  than  35  kt 
(or  MSLP  greater  than  about  1005  hPa)  are  called  tropical  depressions.  Tropical  storms  have 
sustained  winds  of  at  least  35  kt  but  less  than  65  kt  ( MSLP  between  1005  and  987  hPa).  Storms 


Table  2.1:  TC  Intensity  Categorization 


Sustained  Wind  (kt) 

MSLP  (hPa) 

Tropical  depression 

<35 

>1005 

Tropical  storm 

35-64 

987-1005 

Cat  1  Hurricane 

65-82 

975-987 

Cat  2  Hurricane 

83-95 

965-975 

Cat  3  Hurricane 

96-113 

950-965 

Cat  4  Hurricane 

114-135 

925-950 

Cat  5  Hurricane 

>135 

<925 

10 


with  sustained  winds  of  65  kt  or  greater  ( MSLP  <  987  hPa)  are  called— depending  on  ocean 
basin— hurricanes  (Atlantic,  Eastern  Pacific,  and  Central  Pacific  basins),  typhoons  (Western 
North  Pacific  basin),  or  tropical  cyclones  (Indian  Ocean  and  Southern  Hemisphere).  The  Saffir- 
Simpson  scale  subdivides  hurricanes  into  five  categories  according  to  increasing  sustained  wind 
(Table  2.1). 

TC  genesis  is  still  poorly  understood  because  of  the  scarcity  of  in  situ  observations  of 
nascent  storms  and  the  difficulty  in  parameterizing  or  explicitly  modeling  the  convective-scale 
processes  crucial  to  storm  formation  (Ooyama  1982;  Tripoli  1992).  TC's  occur  preferentially  in 
regions  with  clusters  of  tropical  mesoscale  convective  systems  (MCS's),  such  as  tropical  easterly 
waves  (the  most  common  Atlantic  basin  TC  progenitor)  or  the  monsoon  trough  (the  source  of 
most  WNP  TC's)  (Nakayama  2001;  Lander  1994). 

These  MCS's  contain  multiple  organized  convective  cells,  each  of  which  undergoes  a 
formative,  mature,  and  dissipative  stage  (Leary  and  Houze  1979;  Leary  and  Houze  1979b). 
Strong  upward  transport  of  latent  and  sensible  heat  characterizes  the  formative  and  mature 
stages.  Dissipating  convective  cells  leave  behind  anvil  cirrus  cloud  where  ice  crystals  settle, 
melt  upon  reaching  the  freezing  level,  and  fall  as  stratiform  rain.  Latent  heat  release  in  the 
convective  columns  warms  the  upper  troposphere,  while  subsidence  in  the  stratiform  rain 
beneath  the  anvil  brings  cool,  dry  air  from  the  mid-troposphere  down  to  the  surface  (Zipser 
1969).  The  warming  aloft  and  cooling  near  the  surface  stretch  the  air  column,  create  a  potential 
vorticity  anomaly,  and  help  to  generate  a  vortex.  The  heat  carried  aloft  in  the  convective  plumes 
creates  an  environment  neutrally  stable  with  respect  to  moist  convection.  Due  to  the  neutral 
stability,  the  Rossby  radius  of  deformation  is  small  in  this  environment,  so  gravity  waves  cannot 
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efficiently  carry  away  this  heat  (Simpson  et  al.  1997).  The  rotation  of  the  developing  vortex  also 
helps  by  confining  the  heat  released  to  within  a  small  radius  near  the  circulation  center.  Through 
hydrostatic  balance,  this  trapped  upper-tropospheric  warm  air  creates  a  radial  surface  pressure 
gradient  conducive  to  low-level  convergence. 

Further  development,  and  later,  steady-state  intensity,  is  maintained  by  cooperative 
interaction  between  the  primary  (azimuthal)  and  secondary  (radial  and  vertical)  circulations.  The 
primary  circulation  consists  of  cyclonic  low-level  inflow  and  anticyclonic  upper-level  outflow 
(except  near  the  core  of  the  circulation,  where  the  strong  convective  updrafts  carry  cyclonic 
rotation  all  the  way  to  top  of  troposphere).  The  radius  where  the  strongest  cyclonic  circulation 
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Figure  2.1:  Radial-vertical  cross-section  of 
tangential  wind  from  18,000  composited 
rawinsonde  observations  of  248  WNP  storms 
from  1961-1970.  The  cyclonic  circulation 
increases  toward  smaller  radii  and  is 
maximized  at  the  top  of  the  boundary  layer, 
while  the  anticyclonic  circulation  increases 
with  height  and  radius.  From  Frank  (1977). 


occurs  is  in  the  eyewall,  the  ring  of  convection 
that  forms  at  a  radius  of  tens  of  km  from  the 
circulation  center.  The  level  of  strongest 
cyclonic  circulation  is  near  the  top  of  the 
frictional  boundary  layer  at  around  850  hPa. 
The  azimuthal  circulation  weakens  with 
increasing  radius  due  to  the  weakening 
horizontal  pressure  gradient.  It  decreases 
downward  from  the  top  of  the  boundary  layer 
due  to  surface  friction.  It  also  decreases 
upward  from  the  boundary  layer-and 
eventually  reverses  to  become  anticyclonic  at 
large  radii  in  the  upper  troposphere— due  to  the 
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decreasing  radial  pressure  gradient  with 
increasing  height  (Frank  1977;  Emanuel 
2003).  Figure  2.1  is  a  radial-vertical  cross- 
section  of  tangential  wind  from  a  composite  of 
WNP  typhoons. 

The  secondary  circulation  flows 
radially  inward  in  the  boundary  layer,  upward 
in  the  eyewall  and  convective  spiral  bands,  and 
outward  in  the  upper  troposphere.  The  inward 
flow  within  the  boundary  layer  is  caused  by 
frictional  convergence  of  the  cyclonic  low- 
level  primary  circulation.  On  this  leg,  an  air 


Figure  2.2:  The  TC  secondary  circulation, 
consisting  of:  1)  ascent  in  the  eyewall  and 
spiral  band  cumulus  convection;  2)  radial 
outflow  where  longwave  IR  radiation  reduces 
6e\  3)  moist  adiabatic  descent  through  the 
moist  cloud  layer;  4)  radial  inflow  in  the 
boundary  layer  where  sensible  and  latent  heat 
flux  increase  9e;  and  5)  forced  dry  descent 
within  the  eye.  From  Frank  (1977). 


keeps  the  parcel's 


parcel's  angular  momentum  decreases.  Sensible  heat  flux  from  the  ocean 


temperature  nearly  equal  to  the  sea  surface  temperature,  while  latent  heat  flux,  in  the  form  of 
moisture  evaporated  from  the  ocean  surface,  increases  its  equivalent  potential  temperature,  0e. 
When  the  parcel  reaches  the  eyewall  (or  a  spiral  band),  convergence  forces  it  upward  (Emanuel 
1988;  2003).  Condensation  produces  prolific  (-100  mm  hr'1)  rainfall.  Upon  reaching  the  upper 


troposphere,  the  bulk  of  the  ascending  mass  turns  and  flows  radially  outward  in  the  anticyclonic 


outflow.  The  warm  outflowing  parcels  aloft  radiate  longwave  infrared  (IR)  radiation  to  space, 
reducing  their  6e  and  causing  them  to  descend.  When  the  subsiding  parcels  reach  the  mixed 


layer  beneath  the  trade  inversion,  they  descend  moist  adiabatic  ally  to  the  surface,  where  they 
start  the  cycle  over  again.  Figure  2.2  depicts  the  TC  secondary  circulation. 
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Figure  2.3:  Schematic  TC  eye  and  eyewall 
circulation.  Forced  dry  descent  creates  the 
TC's  warm  core.  From  Willoughby  1988. 


Figure  2.4:  Vertical  cross-section  of 
temperature  anomaly  for  Hurricane  Cleo,  18 
Aug  1958.  This  cross-section  was  measured 
by  simultaneous  penetrations  by  three 
research  aircraft  flying  at  the  levels  indicted 
by  dotted  lines.  Cleo's  MSLP  was 
approximately  972  hPa  at  the  time.  From  La 
Seur  and  Hawkins  1963. 


Some  of  the  rising  air  reaching  the  top 
of  the  upward  branch  of  the  secondary 
circulation  is  inertially  confined  within  the 
eyewall,  creating  a  convergent  region  in  the 
upper  troposphere  above  the  core  of  the  TC. 
This  upper-level  convergence  into  the  storm's 
core,  in  concert  with  entrainment  of  lower- 
level  air  out  of  the  eye  and  into  the  eyewall 
convection,  creates  weak  subsidence  in  the 
core.  This  clears  the  eye  of  clouds  in  upper 
and  mid  levels  and  adiabatically  warms  the 
subsiding  air  (Willoughby  1998).  Figure  2.3  is 
a  schematic  of  the  eyewall  secondary 
circulation  and  Figure  2.4  is  a  radial  cross- 
section  of  the  strong  warming  created  by  the 
adiabatic  descent  as  observed  by  research 
aircraft  penetrations  of  Category  2  Hurricane 
Cleo  (1958). 

Since  vertical  accelerations  are  small 
everywhere  except  in  the  convective  updrafts, 
hydrostatic  equilibrium  (where  a  parcel's 
buoyancy  balances  its  weight)  holds  and  the 
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If  we  choose  a  high  enough  zud,  the  height  of  the  corresponding  pressure  surface,  pud ,  will 
be  undisturbed  by  the  TC  below  and  we  can  assume  a  constant  pud  and  zud  for  both  the  TC  and  its 
environment  (Kidder  et  al.  1978).  Then  we  can  determine  the  TC's  MSLP  and  the  environment's 
surface  pressure  penv  using  the  vertically  averaged  temperatures  from  the  eye  and  the  surrounding 
environment.  Combining  these  expressions  yields  MSLP  as  a  function  of  penv  and  the  TC's 
temperature  anomaly: 
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Penv  is  approximately  1000  hPa,  zud  will  be  on  the  order  of  15  km,  and  both  the  TC  and 
environment  vertically-averaged  temperatures  will  be  on  the  order  of  250K  (Velden  and  Smith 


1983).  With  these  scaling  approximations,  and  assuming  a  maximum  TC  temperature  anomaly 
on  the  order  of  10K,  the  linear  term  of  the  Taylor  series  expansion  contributes  tens  of  hPa  to  the 
TC's  MSLP  deficit,  and  the  quadratic  term  contributes  on  the  order  of  1  hPa.  Higher  order  terms 
contribute  negligible  fractions  of  an  hPa.  So  we  obtain  an  approximate  series  solution  for  MSLP 
in  terms  of  the  TC  warm  core  temperature  anomaly: 


MSLP  ~penv- CAT +  C2AT2 


where 
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2.2.  Microwave  Radiative  Transfer  Overview 

Passive  microwave  radiometers,  such  as  AMSU-A  and  -B,  are  good  tools  for  measuring 
the  TC  warm  core  temperature  anomaly.  These  instruments  sense  the  radiation  emitted  by  both 
Earth  and  atmosphere.  Some  of  this  radiation  travels  directly  from  the  emitting  source  to  the 
instrument,  and  at  some  frequencies,  a  significant  fraction  is  reflected  by  the  ocean  surface,  or 
scattered  by  precipitation  in  the  atmosphere. 

The  fundamental  quantity  used  in  radiative  transfer  is  radiant  intensity,  or  radiance,  the 
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amount  of  power  incident  on  a  unit  surface  area  per  interval  of  wavelength,  traveling  in  a 
direction  normal  to  the  surface,  and  originating  from  a  unit  of  solid  angle.  The  units  of  radiance 
are  W  m2  /urn1  sr'1  and  the  radiance  emitted  at  wavelength  A  by  a  blackbody  at  temperature  T  is 
given  by  the  Planck  function: 


B(A,T)=- 


2  he2 


exp 


he 


\kBAT  j 


h  =  6.626x1  O'34  Js 

23  JR'1 
c  =  2.998x10®  ms-1 


- — ^  ,  kB  =1.381x10 
-1 


The  Planck  radiance  at  all  wavelengths  increases  with  temperature  and  the  peak  of  the 
Planck  curve  moves  to  shorter  wavelengths  with  increasing  temperature.  The  Planck  curve  for 
the  sun,  at  about  6000  K,  peaks  in  the  visible  portion  of  the  spectrum,  while  the  earth,  at 
approximately  250  K,  has  a  curve  that  peaks  in  the  far-IR.  So  for  terrestrial  remote  sensing 
problems  in  the  visible  and  IR,  the  relationship  between  an  object's  temperature  and  its  Planck 
radiance  is  highly  non-linear.  The  microwave  portion  of  the  spectrum,  occupying  wavelengths 
between  1  mm  and  10  cm  (frequencies  between  3  GHz  and  300  GHz),  is  well  into  the  long- 
wavelength  tail  of  the  Planck  curve  for  emitters  at  terrestrial  temperatures.  For  this  part  of  the 
curve,  the  Rayleigh-Jeans  approximation  simplifies  the  Planck  function  and  defines  a  nearly 
linear  relationship  between  radiance  and  temperature  at  a  given  wavelength: 


B(X,T)*^LT 


The  brightness  temperature,  TB,  of  a  given  emitter  is  the  temperature  of  a  blackbody  emitting  the 
same  radiance.  When  the  Rayleigh-Jeans  approximation  is  valid,  as  it  is  in  terrestrial  microwave 
remote  sensing  applications,  TB  can  then  be  used  interchangeably  with  radiance.  An  object's 
emissivity,  e,  is  the  ratio  of  its  actual  temperature  to  its  brightness  temperature,  so  that  TB  =  eT. 
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To  understand  how  microwave  radiation  emitted  by  the  earth  and  atmosphere  travels 
upward  through  the  atmosphere  to  a  sensor  in  orbit,  we  begin  with  Beer's  Law,  which  describes 
how  the  intensity  at  a  given  wavelength,  /;.,  decreases  as  radiation  travels  a  path  between  points 
sj  and  $2  in  a  medium  that  absorbs  but  does  not  scatter: 

/(s2)  =  /(^)exp^- 

^1.^)  =  [A(^s 

=  exp(-  t(j„  j2)) 

pa  is  the  absorption  coefficient  of  the  medium  and  has  units  of  inverse  length.  For  an  individual 
atmospheric  constituent,  pa  is  the  product  of  its  density  and  its  mass  absorption  coefficient,  ka 
(units  of  length2  per  mass),  which  is  frequency-dependent.  The  absorption  coefficient  is 
additive,  so  pa  for  a  mixture  of  gases  is  the  sum  of  all  individual  Pd s.  r(S],s2)  is  the  optical 
thickness  between  sj  and  s2,  and  t(sj,S2)  is  the  transmittance,  the  fraction  of  radiation  originating 
at  sj  that  reaches  s2 .  The  optical  depth  at  a  height  z  in  the  atmosphere  is  the  optical  thickness 
from  the  top  of  the  atmosphere  vertically  downward  to  z: 

t(z)  =  Pa{z')dz' 

If  we  define  zenith  angle,  0,  as  the  angle  between  vertical  and  the  propagation  direction  of  a 
beam,  and  further  define  p  =  cos(O),  the  transmittance  between  optical  depth  r(z)  and  the  top  of 
the  atmosphere  is: 

t(z)  =  expf-^^ 

l  A 

While  absorption  will  decrease  a  beam's  intensity  as  it  traverses  a  medium,  thermal 
emission  by  the  medium  itself  will  add  intensity  (by  Kirchhoff  s  law,  absorption  equals  emission 
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in  local  thermodynamic  equilibrium,  so  we  can  use  the  absorption  coefficient  in  conjunction  with 
the  Planck  function,  B(X,T),  to  represent  emission  by  the  medium): 

~r-  =  "A  {s)l  x  (0 + A  (s)b{^,T(s)) 
as 

If  we  measure  the  optical  thickness  of  the  medium  backward  along  the  path  of  propagation  from 
the  sensor  toward  the  point  of  emission,  so  that  dr  =  -fie  ds,  we  can  transform  the  radiative 
transfer  equation  from  units  of  geometric  length  to  units  of  optical  thickness: 

^-=Ix{t)-B(XT(t)). 

dT 

By  multiplying  each  term  by  an  integrating  factor  e'T(e.g.,  Petty  2004)  we  can  solve  for  intensity 
at  the  top  of  the  atmosphere  (TOA;  where  r  =  0)  in  terms  of  intensity  originating  at  the  surface 
(where  r-r*): 

dt 

~{lx{T)e-')  =  -B(X,T{T))e-' 
dr 

WY* )  =  -  f  B(Z,  Tit^'dr 

Ix{t  =  t - 1 x {t  =  0)  =  -  [  B{A, T(t'))e-rdT' 

Itoax  =  hFcx^  +  f  B{A,T{?))e-rdT' 

Substituting  transmittance,  t  =  eT  and  dt  =  -ex  dr,  for  optical  depth  as  the  vertical  coordinate,  and 
using  the  Rayleigh-Jeans  approximation  to  exchange  TB  for  h;. 

TBtoa  =TBsfct*+\T{r)df 

Here,  t*  is  the  transmittance  from  the  surface  to  the  top  of  the  atmosphere. 
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Air  density  decreases  exponentially  with  altitude,  with  a  scale  height  of  H : 

/ 

p{z)  =  p0  exp 

v 

So  the  transmittance  from  height  z  to  the  top  of  such  an  atmosphere  is  given  by: 

fia{z)  =  K 

^{z)  =  P  ok  a  exP^~ 

=  p0kaH 


f  _z_A 
\  lij 


dz' 


t(z)  =  exp 


PokaH  f  z ' 

— exp - 

M  l  H, 


Using  this  expression  for  transmittance  in  an  exponential  atmosphere,  we  can  transform 
to  geometric  height  coordinates: 


TBT0A=TBsfct*  +  [T{z)^dz 


t(z )  =  exp 
dt  _  p0k 


PokaH 


P 


exp|~¥ 


dz  p 


-exp 


PoKH 

P 


exp 


' 

z 

H 

=  W(z) 


TBT0A=TBsfct*+%T{zMz)dz 


Note  that  the  surface  brightness  temperature,  TBSfC,  consists  not  just  of  Planck  emission  from  the 
surface  itself,  but  also  reflected  downwelling  atmospheric  emission  and  cosmic  background 
radiation.  The  downwelling  atmospheric  TB  can  be  derived  in  a  similar  manner  to  the  upwelling 
emission. 

W(z)  is  the  emission  weighting  function,  which  describes  the  contribution  of  each  layer's 
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temperature  to  the  upwelling  radiance  at  the 

dt 

top  of  the  atmosphere.  Since  W  (z)  =  — , 

dz 

layers  that  contribute  most  to  the  upwelling 
radiance  are  those  where  the  transmittance  to 
the  top  of  the  atmosphere  changes  most 
rapidly  with  height.  In  an  exponential 
atmosphere,  the  upper  levels  have  such  low 
density  that  t(z)  is  near  1  and  dt/dz  =  W(z)  ~  0. 
Deeper  in  the  atmosphere,  the  density  may  be 
large  enough  that  the  atmosphere  is  opaque,  so 
t(z)  is  near  0  and  again  dt/dz  -  W(z)  ~  0.  The 
layers  with  large  W(z)  are  those  which  absorb 
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Figure  2.5:  Air  density  (upper  left),  optical 
depth  (upper  right),  transmittance  (lower 
left),  and  weighting  function  (lower-right)  for 
an  atmosphere  with  a  surface  density  of  1 .0 
kg  m'3,  an  average  temperature  of  250  K,  and 
a  mass  absorption  coefficient  of  .0006  m2  kg' 
\  Transmittance  and  weighting  function  are 
plotted  for  nadir  incidence  (solid  curve)  and 
48°  incidence  (dashed  curve). 


radiation  from  below,  but  whose  own  emitted  radiation  isn’t  absorbed  as  strongly  by  layers 


above.  Figure  2.5  depicts  weighting  functions  for  two  different  zenith  angles  at  a  fixed  mass 


absorption  coefficient. 

The  shape  of  W(z)  depends  primarily  on  the  mass  absorption  coefficient,  ka(X),  and  zenith 
angle,  in  the  form  of  ft.  ka(X)  has  very  fine-scale  variability  in  the  IR  portion  of  the  spectrum, 
with  many  intermixed  vibrational  transitions  of  O2,  water  vapor  (WV),  CO2,  O3,  and  other 
constituents.  Radiative  transfer  computations  in  these  bands  require  careful  handling  of  the 
multitude  of  overlapping  bands.  In  the  microwave,  though,  there  are  only  3  significant  absorbing 
processes.  The  first  is  rotational  absorption  by  molecular  oxygen.  This  occurs  for  a  broad  band 
from  51-67  GHz  and  a  narrow  band  at  1 18.75  GHz.  Rotational  absorption  by  water  vapor  also 
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contributes  two  absorption  bands,  a  weak  band  at  22.24  GHz,  and  a  strong  one  at  183.31  GHz 
(Rosenkranz  1993).  The  third  process  is  continuum  absorption  by  water  vapor.  The  exact 
mechanism  of  WV  continuum  absorption  is  not  well  known,  but  may  be  caused  either  by  the 
sum  of  overlapping  tails  of  the  multitude  of  pressure-  and  Doppler-broadened  WV  absorption 
bands  from  the  IR  spectrum,  or  may  be  caused  by  anomalous  chains  of  WV  molecules  adhering 
to  each  other  (Petty  2004).  Whatever  the  specific  mechanism,  WV  continuum  absorption 
increases  with  increasing  microwave  frequency. 

O2  is  a  fixed  atmospheric  constituent.  In  the  well-mixed  lower  and  middle  atmosphere,  it 
comprises  a  constant  21%  of  air  by  volume.  WV,  on  the  other  hand,  makes  up  0  to  2%  of  air  by 
volume,  but  its  concentration  varies  strongly  in  both  the  vertical  and  horizontal.  The  total 
absorption  by  both  constituents  is  additive,  so  the  net  microwave  absorption,  particularly  near  the 
two  WV  rotation  bands  and  at  higher  microwave  frequencies  where  continuum  absorption  is 
strong,  is  dependent  on  the  WV  concentration.  Figure  2.6  depicts  the  zenith  microwave 
transmittance  for  polar,  standard,  and  tropical  atmospheres,  illustrating  the  striking  difference 
varying  WV  content  can  have  on  microwave  transmittance. 

Microwave  sounding  instruments  make  use  of  these  absorption  features  in  order  to 
sample  the  vertical  temperature  and  moisture  profile.  On  the  fringes  of  the  O2  and  WV 
absorption  bands,  ka~and  therefore  t(z)-change  rapidly  with  frequency.  Thus  the  weighting 
functions,  W(z),  for  frequencies  on  the  edges  of  these  bands  will  peak  progressively  higher  in  the 
atmosphere  as  frequency  approaches  the  central  frequency  of  the  absorption  band.  Temperature 
sounding  instruments  typically  use  the  60  GHz  O2  band-with  several  channels  closely  spaced 
between  50-58  GHz-while  moisture  sounding  instruments  use  the  183  GHz  WV  band.  As 


shown  in  Figure  2.6,  the  width  and  strength  of  the  WV  bands  are  strongly  modulated  by 
atmospheric  WV  concentration,  so  the  shape  of  near- 183  GHz  weighting  functions  will  depend 
strongly  on  humidity.  More  moisture  decreases  transmittance  and  causes  W(z)  to  peak  higher. 

Transmittance  depends  on  path  length  through  the  absorbing  medium,  so  as  scan  angle 
increases  away  from  nadir  and  the  slant  path  through  the  atmosphere  increases  with  it, 
transmittance  will  decrease  and  W(z)  will  peak  higher.  This  effect,  limb-darkening,  can  cause 
brightness  temperatures  to  drop  by  as  much  as  10  K  from  nadir  to  a  scan  angle  of  45°.  For 


Figure  2.6:  Microwave  transmittance  spectrum  under  varying  water  vapor  concentrations. 
The  60  and  118  GHz  O2  rotation  bands  are  unchanged,  due  to  the  fixed  O2  concentration. 
The  weak  22  GHz  and  strong  183  GHz  WV  rotation  bands,  increase  in  intensity  as  WV 
increases.  Most  significant,  though,  is  the  marked  increase  in  continuum  absorption,  which 
makes  the  atmosphere  nearly  opaque  to  frequencies  above  170  GHz  in  the  tropical 
atmosphere.  From  Grody  (1976). 
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frequencies  whose  weighting  functions  contain  a  significant  contribution  from  the  stratosphere 
(where  temperature  increases  with  height)  the  limb  effect  increases  TB  instead. 

At  frequencies  where  the  atmosphere  has  high  transmittance,  such  as  the  "window" 
regions  below  22  GHz,  between  22  GHz  and  50  GHz,  between  70  GHz  and  1 10  GHz,  and 
between  125  GHz  and  170  GHz,  a  downward  looking  instrument  will  be  able  to  see  through  the 
(mostly)  transparent  atmosphere  and  will  sense  the  upwelling  brightness  temperature  from  the 
surface.  This  has  3  components:  thermal  emission  by  the  surface  itself,  reflected  downward 
emission  from  the  atmosphere,  and  reflected  cosmic  background  radiation.  The  cosmic 
background  emits  at  a  blackbody  temperature  of  2.7  K,  so  for  most  purposes,  it  can  be  ignored 
compared  with  terrestrial  emission  sources,  which  all  have  blackbody  temperatures  of  100-300 

K.  The  downward  atmospheric  emission  is 
computed  in  a  manner  identical  to  upward 
emission,  but  the  weighting  function  in  this 
case  is  computed  by  integrating  upward  from 
the  surface  to  z  instead  of  downward  from 
TOA  to  z- 

The  reflection  of  downward 
atmospheric  emission,  thermal  emission  by  the 
surface,  and  scattering  by  precipitation 
(discussed  in  Ch.  3)  bring  polarization  into 
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Figure  2.7:  Surface  reflectivity  (left  column) 

and  emissivity  (right  column)  for  3 1 .4  GHz  play  for  microwave  remote  sensing.  At  these 

and  89.0  GHz  window  frequencies.  Solid 

line  is  for  vertical  polarization;  dashed  line  is  frequencies,  ocean  surface  emissivity  is  greater 
for  horizontal. 
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for  vertically-polarized  radiation  than  for  horizontal,  so  for  window  channels  that  can  see  to  the 
surface,  the  observed  TB  will  be  polarized.  The  degree  of  polarization  of  sea-surface  emission 
depends  on  scan  angle,  since  the  vertically-  and  horizontally-polarized  emissivities  have  different 
scan-angle  dependences.  Vertically-polarized  emissivity  increases  from  its  nadir  value  to  a 
maximum  at  the  Brewster  angle,  then  decreases  to  0  at  90°  incidence,  while  horizontally- 
polarized  emissivity  decreases  smoothly  from  its  maximum  at  nadir  (where  Sh  -  cv)  to  0  at  90° 
incidence.  Both  vertically-  and  horizontally-polarized  emissivity  increase  with  increasing 
frequency.  Figure  2.7  depicts  ocean  surface  emissivity,  s,  and  reflectivity,  r.  By  Kirchhoff  s 
Law,  r  =  1  -  s.  The  surface  emissivity  for  window  frequencies  at  scan  angles  of  <  45°  ranges 
from  around  0.45  to  0.75,  so  for  clear  scenes,  the  ocean  surface  brightness  temperature  will  range 
from  about  150  K  to  225  K.  Microwave  land  surface  emissivity  varies  with  vegetation  type  and 
soil  moisture,  but  in  general  is  much  higher,  and  less  polarized,  than  ocean-surface  emissivity. 
Typical  land  TB' s  are  around  280  K. 

Retrieving  atmospheric  temperature  profiles  from  multi-channel  microwave  TB' s  can  be 
done  statistically  or  physically.  The  statistical  approach  uses  a  large  database  of  co-located 
rawinsonde  and  microwave  observations  segregated  by  latitude,  surface  type,  and  time  of  year. 
Regression  analysis  produces  matrices  containing  coefficients  relating  each  microwave  TB  to  the 
temperature  at  each  desired  atmospheric  level  (Grody  1993;  Goldberg  1999;  Reale  2001).  The 
physical  approach  is  a  data  assimilation  problem.  A  radiative  transfer  model  (the  forward 
model)  predicts  the  expected  TOA  TB's  resulting  from  a  first-guess  temperature  profile  obtained 
from  climatology  or  a  numerical  weather  prediction  model  forecast.  The  temperature  profile  is 
iteratively  adjusted  using  adjoint  techniques  to  minimize  a  cost  function,  which  quantifies  both 


25 


the  deviation  of  the  derived  profile  from  the  first-guess  and  the  discrepancy  between  predicted 
and  observed  TB's  (Chahine  1970;  Smith  1970;  Rodgers  1976). 

2.3.  Previous  Microwave  Sounding  Instruments  and  TC  Intensity  Techniques 

Polar-orbiting  passive  microwave  temperature  sounding  instruments  have  been  in  use  for 
over  three  decades.  The  earliest,  the  Nimbus  Experimental  Microwave  Spectrometer  (NEMS), 
flew  on  the  Nimbus-5  research  mission  in  1972.  It  sampled  a  single  200  km  nadir-looking  spot 
at  three  temperature  sounding  channels  (53.65,  54.90,  and  58.80  GHz).  The  Scanning 
Microwave  Spectrometer  (SCAMS)  flew  experimentally  on  Nimbus-6,  launched  in  1975. 
SCAMS  sampled  two  window  frequencies  (22.235  and  31.65  GHz)  and  3  temperature  sounding 
frequencies  (52.85,  53.85,  and  55.45  GHz)  in  13  fields  of  view  (FOV)  across  each  cross-track 
scan  line,  with  resolution  ranging  from  145  km  to  360  km.  The  first  operational  temperature 
sounder  was  the  Microwave  Sounding  Unit  (MSU),  nine  copies  of  which  flew  on  NOAA-6, 
launched  in  1979,  through  NOAA-14,  whose  instrument  is  still  operational  today.  MSU 
employed  4  temperature  sounding  channels  (50.30,  53.74,  54.96,  and  57.95  GHz)  with  cross¬ 
track  resolution  ranging  from  110  km  to  270  km.  The  Special  Sensor  Microwave  Temperature 
(SSMT)  instrument,  a  contemporary  of-and  very  similar  to-the  MSU,  also  flew  on  the  Defense 
Meteorological  Satellite  Program  (DMSP)  polar-orbiting  constellation  beginning  in  1979.  The 
current  generation  sounding  instrument  is  the  Advanced  Microwave  Sounding  Unit  (AMSU), 
which  first  flew  on  NOAA-15  in  1998.  A  detailed  discussion  of  the  AMSU-A  and  -B 
instruments  is  contained  in  Section  2.4.  Future  MW  temperature  sounding  instruments  include 
the  Special  Sensor  Microwave  Imager/Sounder  (SSMIS;  more  detailed  discussion  in  Chapter  7), 
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first  launched  on  DMSP  F16  in  late  2003,  but 
not  yet  operational,  and  the  Conically- 
Scanning  Microwave  Imager/Sounder  (CMIS), 
to  fly  on  the  National  Polar-orbiting 
Operational  Environmental  Satellite  System 
(NPOESS)  constellation  beginning  around 
2010. 

As  discussed  in  Section  2.1,  the  TC 
warm  anomaly  is  strongest  in  the  upper 
troposphere.  Based  on  the  microwave 
temperature  sounding  considerations  in 
Section  2.2,  we  expect  microwave  frequencies  near  55  GHz  to  have  weighting  functions  that 
peak  near  this  level,  as  shown  in  Figure  2.8.  Each  generation  of  sounding  instrument  since 
SCAMS  has  included  at  least  one  channel  near  this  frequency,  affording  maximum  sensitivity  to 
the  TC  warm  core,  and  making  them  useful  for  estimating  MSLP  via  the  TC  warm  anomaly. 
Section  2.3  briefly  discusses  these  previous  efforts. 

Kidder  et  al.  (1978)  were  the  first  to  estimate  TC  intensity  using  a  microwave  sounder. 
They  derived  a  linear  MSLP  vs.  TB  anomaly  relationship  for  an  upper-tropospheric  sounding 
channel  (55.45  GHz),  and  obtained  root  mean  square  error  of  15  hPa  (Figure  2.9)  for  a  set  of  36 
Western  North  Pacific  TC  observations.  In  a  later  paper  (1980),  the  same  authors  derived  similar 
relations  for  the  radii  of  30  kt  and  50  kt  sustained  winds  using  TB  gradients.  They  noted  that 
poor  spatial  resolution  and  large  instrument  noise  created  the  bulk  of  estimate  uncertainty,  but 


Figure  2.8:  Radial-vertical  cross-section 
(left)  of  composite  typhoon  from  Frank 
(1977).  AMSU-A  weighting  functions 
(right)  at  nadir  (solid)  and  edge  of  scan 
(dashed)  from  Goldberg  et  al.  (2001). 
AMSU-A  Channel  7  matches  most  closely 
the  level  (approximately  300  hPa)  where  the 
TC  warm  anomaly  is  largest. 
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Figure  2.9:  SCAMS  55.45  GHz  TB  anomaly  observed  by  Kidder  et  al.  (1978)  for  Typhoon 
June  at  1426  UTC  on  19  Nov  1975  (left)  and  MSLP  vs.  SCAMS  55.45  GHz  TB  anomaly 
(right). 

also  speculated  that  hydrometeor  scattering  might  be  a  source  of  uncertainty. 

Velden  (1982),  Velden  and  Smith  (1983),  and  Velden  et  al.  (1984)  developed  a  similar 
technique  using  the  Laplacian  of  the  250  hPa  retrieved  temperature  field  from  the  MSU.  They 
placed  more  emphasis  on  precipitation  scattering  as  a  source  of  uncertainty,  and  manually  edited 
fields-of-view  (FOV)  which  were  obviously  affected  (Figure  2.10).  In  subsequent  studies, 
Velden  (1989)  and  Velden  et  al.  (1991)  achieved  RMSE  of  8  and  14  hPa  in  the  Atlantic  and 
Northwest  Pacific  basins,  respectively. 


Figure  2.10:  250  hPa  retrieved  temperature  field  observed  by  Velden  and  Smith  (1983)  for 
Hurricane  Harvey  at  1800UTC  on  15  Sep  1981  (left).  Note  the  fields  of  view  near  storm 
center  manually  removed  due  to  precipitation  attenuation.  T250  Laplacian  vs.  MSLP  (right). 


A  significant  leap  in  instrument  capability  came  in  1998  with  the  launch  of  the  first 
AMSU-A  and  -B  suite  on  NOAA-15.  Brueske  (2001)  and  Brueske  and  Velden  (2003) 


developed  an  MSLP  estimation  technique  at  CIMSS  using  AMSU-A  channel  7  (54.94  GHz, 
sensitive  to  250  hPa  temperature)  TB  anomaly,  along  with  Goldberg's  (2001)  correction  for  limb 


effects  and  Merrill's  (1995;  Van  Burgel  1999) 
physically-based  warm  anomaly  retrieval 
scheme.  Brueske  and  Velden  achieved  5  hPa 
RMSE,  but  for  a  small  (n  ~  50)  sample  (Figure 
2.11),  and  established  CIMSS'  automated  real¬ 
time  intensity  bulletins.  Kabat  (2002)  added 
AMSU-A  TB8  (55.5  GHz,  150  hPa)  anomaly 
as  an  intensity  predictor,  partly  to  improve 
performance  in  cases  where  TB7  is  attenuated 


M3U»(hP») 


Figure  2.11:  AMSU-A  raw  and  retrieved 
TB7  anomaly  vs.  MSLP  for  the  Atlantic  and 
East  Pacific  basins  for  the  2000  TC  season 
(Brueske  and  Velden  2003). 
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by  scattering.  With  the  addition  of  corrections  for  sub-sampling  and  non-standard  environmental 
pressure,  current  performance  of  the  CIMSS  AMSU-based  method  ranges  from  6  hPa  RMSE  in 
the  Atlantic  to  10  hPa  RMSE  in  the  Indian  Ocean  (Herndon  et  al.  2004).  However,  the  CIMSS 
technique  does  not  yet  account  for  precipitation  attenuation  effects. 

Spencer  and  Braswell  (2001)  developed  an  AMSU-A  Vmax  estimation  technique  using  the 
radial  gradient  of  AMSU-A  TB8.  They  recognized  the  importance  of  hydrometeor  scattering  and 
accounted  for  it  by  including  spatial  gradients  of  both  TB4  and  an  AMSU-A  window  channel 
scattering  index  (discussed  in  Chapter  3)  among  their  predictors.  They  achieved  5  m  s'1 
(equivalent  to  about  10  hPa  for  strong  storms  or  5  hPa  for  weak  storms)  RMSE  (Figure  2.12). 
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Fig.  2.12:  (Left)  AMSU-A  radial  TB  gradient  by  channel  for  3  different  TC's.  Positive  TB 
gradient  for  TB6  through  TB9  indicates  the  TC  warm  core  while  TB1  through  TB5  are 
dominated  by  hydrometeor  scattering  effects.  (Right)  Estimated  Vmax  vs.  aircraft- 
reconnaissance  measured  Vmax.  From  Spencer  and  Braswell  (2001). 
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Zhu  et  al.  (2002;  2004)  and  Demuth  et 
al.  (2004)  use  3-D  fields  of  AMSU-A  retrieved 
temperatures  to  derive  the  corresponding  3-D 
pressure  and  wind  fields  using  hydrostatic  and 
gradient  wind  balance,  with  boundary 
conditions  obtained  from  NWP  analyses.  Zhu 
et  al.  employed  these  temperature  fields  to 
improve  NWP  forecast  initialization.  Demuth 
et  al.  employed  their  fields  to  derive  statistical 
TC  intensity  predictors,  from  which  they 
produce  real-time  TC  intensity  estimates, 
complimenting  the  estimates  from  CIMSS. 

Zhu  et  al.  address  precipitation  effects  with  a 
physical  retrieval  scheme  whose  forward 
model  accounts  for  emission  and  scattering. 
Demuth  et  al.  use  a  statistical  temperature 
retrieval,  with  a  climatological  first  guess. 
They  correct  the  TB's  used  in  the  retrieval  for 
precipitation  effects  using  a  two-step 


Figure  2.13:  Non-scattering-corrected  (left 
column)  vs.  scattering  corrected  (right 
column)  radial-vertical  cross-section  of 
radially-averaged  temperature  anomaly  (top 
row),  surface  pressure  (middle  row),  and 
radial-vertical  cross-section  of  radially- 
averaged  tangential  wind  (bottom  row). 
Anomalously  cold  lower-tropospheric 
sounding  channel  TB's  resulting  from 
precipitation  attenuation  cause  cold 
anomalies  in  the  low-level  retrieved 
temperature  field,  producing  a  spurious  high 
pressure  area  and  associated  wind  field 
perturbations.  From  DeMuth  et  al.  2004. 


procedure:  First,  they  use  the  linear  relationship  between  AMSU-A  cloud  liquid  water  estimates 


(Weng  et  al.  2003)  and  TB  depression  to  correct  the  instrument-resolution  TB  field.  Second,  they 


define  gridded  temperature  fields  more  that  0.5  K  colder  than  their  neighbors  as  affected  by  ice 
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scattering,  and  correct  those  fields  by  smoothing  (Figure  2.13).  DeMuth  et  al.  achieved  9  hPa 
RMSE. 

2.4.  The  Advanced  Microwave  Sounding  Unit 

AMSU-A  and  -B  are  cross-track  scanning  microwave  temperature  and  moisture  sounding 
instruments,  built  to  fly  on  the  NOAA  Advanced  TIROS-N  class  of  spacecraft  that  operate  as  the 
U.S.'s  Polar-orbiting  Operational  Environmental  Satellite  (POES)  constellation.  Together  with 
the  High-resolution  Infrared  Sounder,  HIRS/3,  they  comprise  the  Advanced  TIROS  Operational 
Vertical  Sounder  (ATOVS)  system.  The  HIRS  has  finer  vertical  resolution,  but  since  most 
clouds  are  opaque  in  the  IR  spectrum,  it  can  only  be  used  to  produce  temperature  soundings  in 
clear  scenes.  For  microwave  instruments  like  AMSU-A  and  -B,  cirrus  clouds  are  effectively 
transparent  and  other  clouds  are  only  weakly  absorbing  (discussed  in  more  detail  in  Chapter  3). 
This  enables  AMSU-A  to  produce  temperature  under  cloudy  conditions,  but  with  less  vertical 

Thermal  Control 


Figure  2.14:  NOAA-KLM  series  polar- 
orbiting  spacecraft.  AMSU-A  1,  -A2,  and  -B 
units  are  circled  on  the  Earth-facing  side  of 
the  spacecraft.  Image  from 
http://www.cira.colostate.edu/ 
ramm/hillger/NOAA-K_sketch.gif. 


resolution  than  HIRS  clear  soundings.  The 
AMSU-B  is  used  to  produce  humidity  profiles. 

AMSU-A  itself  consists  of  two 
separate  modules.  AMSU-A1  contains 
antennas  measuring  TB's  in  the  50-60  GHz  O2 
band  and  the  89.0  GHz  window  frequency, 
while  AMSU-A2  contains  antennas  for  the 
23.8  and  31.4  GHz  window  frequencies.  Both 
AMSU-A  and  AMSU-B  utilize  parabolic 
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reflectors,  which  are  rotated  360°  around  an  axis  to  reflect  radiation  into  a  detector.  When 
viewing  the  earth,  AMSU-A  steps  through  30  beam  positions  in  6  seconds,  then  sweeps  past  cold 
space  and  warm  blackbody  views  in  2  seconds  to  provide  cold  and  warm  calibration  limits. 
AMSU-B  completes  3  scans  for  each  AMSU-A  scan,  and  sweeps  continuously  across  the  scan 
line,  rather  than  stepping  and  staring.  AMSU-A's  half-power  beam-width  is  3  1/3°,  which  from  a 
nominal  spacecraft  altitude  of  833  km,  yields  a  nearly  circular  48  km  FOV  for  beam  positions  15 
and  16  (1  2/3°  each  side  of  nadir).  The  FOV  grows  to  a  150  km  x  80  km  ellipse  at  beam 
positions  1  and  30  (48  1/3°  each  side  of  nadir).  AMSU-B's  half-power  beam  width  is  1.1°,  giving 
it  3  times  finer  horizontal  resolution  than  AMSU-A.  The  total  swath  width  of  both  instruments  is 
approximately  2100  km. 

The  AMSU-A  1,  -A2,  and  -B  detectors  are  sensitive  to  a  single  polarization-horizontal  or 
vertical-but  due  to  the  rotation  of  their  reflectors,  the  polarization  plane  of  incoming  radiation  to 
which  each  is  sensitive  rotates  with  scan  angle.  The  vertical  polarization  plane  is  defined  by  the 
local  zenith  direction  and  the  instrument  line  of  sight.  The  horizontal  polarization  plane  is 
perpendicular  to  the  vertical  plane.  Channels  vertically-polarized  (polarization  angle  of  90°)  at 
nadir  (AMSU-A  channels  1-4,  7,  and  15;  AMSU-B  channels  16-19)  have  a  polarization  angle  of 
(90°  -  0)  at  scan  angle  0,  while  horizontally-polarized  channels  (AMSU-A  channels  5,  6,  and  8- 
14;  AMSU-B  channel  20)  have  a  polarization  angle  of  0. 

The  instrument  measures  antenna  temperature,  the  total  radiance  detected  by  the 
reflector's  main  lobe  and  side  lobes.  This  radiance  includes  components  emitted  by  the  earth  and 
the  cosmic  background,  and  emitted  or  reflected  by  the  spacecraft  itself.  Mo  (1999)  modeled 
these  contributions  and  found  they  cause  the  antenna  temperature  to  vary  smoothly  by  1-2  K 


Pressure  (mb) 
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Figure  2.15:  AMSU -A  temperature  sounding 
weighting  functions  at  near-nadir  (fields  of 
view  15  and  16;  solid  curves)  and  edge  of 
scan  (fields  of  view  1  and  30;  dashed  curves). 
From  Goldberg  et  al.  2001. 


across  the  scan  line.  Removing  the  error 
introduced  by  the  cold  space  and  spacecraft 
contributions  results  in  the  brightness 
temperature. 

The  radiometric  accuracy  specification 
(expressed  at  noise-equivalent  temperature 
uncertainty,  NEAT)  varies  with  channel,  but  in 
general  becomes  less  stringent  with  increasing 
frequency.  For  AMSU-A  tropospheric 
temperature  sounding  channels,  NEAT  is  0.25 
K. 

Table  2.2  summarizes  the  frequencies, 
radiometric  accuracy,  polarizations,  and 
primary  operational  application  for  each 


AMSU-A  and  -B  channel.  Figure  2.15  depicts  the  weighting  functions  for  each  AMSU-A 
temperature  sounding  channel. 
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Table  2.2:  AMSU-A  and  -B  instrument  characteristics 


Channel 

Frequency  (GHz) 

NEAT 

(K) 

Pol 

Application 

1 

23.8 

.30 

90° -0 

Cloud  liquid,  precip 

2 

31.4 

.30 

90° -0 

Cloud  liquid,  precip 

3 

50.3 

.40 

90° -0 

Cloud  liquid,  precip 

4 

52.8 

.25 

90° -0 

Tropospheric  (sfc)  temp 

5 

53.596 

.25 

0 

Tropospheric  (700  hPa)  temp 

.  6 

54.4 

.25 

0 

Tropospheric  (400  hPa)  temp 

5  ^ 

54.94 

.25 

90° -0 

Tropospheric  (250  hPa)  temp 

55.5 

.25 

0 

Tropospheric  (150  hPa)  temp 

1  9 

57.290344 

.25 

0 

10 

57.290344  ±  .217 

.40 

0 

11 

57.290344  +  .3222  +  .048 

.40 

0 

12 

57.290344  +  .3222  +  .022 

.60 

0 

Stratospheric  temp 

13 

57.290344  ±  .3222  ±  .010 

.80 

0 

14 

57.290344  ±  .3222  ±  .0045 

1.20 

0 

15 

89.0 

.50 

90° -0 

Cloud  liquid,  precip 

16 

89.0 

1.0 

90° -0 

Cloud  liquid,  precip 

S  17 

150.0 

1.0 

90° -0 

Ice  scattering 

w  18 

183.3  ±  1.0 

1.1 

90° -0 

Upper  tropospheric  humidity 

<  19 

183.3  +  3.0 

1.0 

90° -0 

Middle  tropospheric  humidity 

20 

183.3  +  7.0 

1.2 

0 

Lower  tropospheric  humidity 

Notes: 

1.  90°  polarization  is  vertical,  0°  is  horizontal 

2.  0  indicates  scan  angle  from  nadir  (ranges  from  0°  to  50°); 

3.  From  NOAA/NESDIS  NOAA-KLM  Instrument  Guide,  2000 
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3.  Precipitation  Effects  in  the  Microwave  Spectrum 

The  goal  of  this  study  is  to  correct  for  precipitation  effects  on  AMSU-A  sounding 
channel  brightness  temperatures.  This  chapter  applies  Mie  scattering  theory  to  hydrometeors  and 
the  microwave  spectrum,  examines  previous  studies  of  microwave  precipitation  effects,  and 
presents  a  1-dimensional  conceptual  model  demonstrating  the  feasibility  of  the  approach  to 
precipitation  correction  developed  in  later  chapters. 

3.1.  Mie  Theory 

In  order  to  understand  how  hydrometeors— liquid  or  frozen  water  suspended  in  the 
atmosphere-interact  with  microwave  radiation,  we  need  to  understand  the  theory  governing  the 
interaction  of  electromagnetic  waves  and  particles  they  encounter.  When  such  an  encounter 
occurs  the  incident  wave  can  create  dipole  moments  within  its  constituent  molecules,  which  will 
in  turn  radiate  their  own  electromagnetic  fields.  The  resulting  field  within  and  outside  the 
particle  is  the  sum  of  the  incident  field  and  the  re-radiated  fields  from  all  the  dipoles  comprising 
the  particle.  Some  of  the  molecules  within  the  particle  will  absorb  energy  without  re-radiating, 
and  interference  among  the  fields  from  the  dipoles  that  do  radiate  and  the  incident  field  will 
cause  energy  to  be  radiated  preferentially  in  certain  directions.  So  particles  absorb  some  of  the 
incident  energy,  and  scatter  the  rest  in  directions  different  from  the  incident  direction. 

Lorentz,  Mie,  and  Debye  each  arrived  separately  at  similar  analytical  solutions  for 
absorption  and  scattering  in  the  special  case  of  spherical  symmetry  (Bohren  and  Huffman  19S3). 
Today,  most  refer  to  this  solution  as  Mie  theory.  More  complex  particle  shapes  require  a  model 
of  the  individual  dipoles  in  the  particle  and  explicitly  account  for  the  interaction  of  each  re- 
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radiated  field.  The  discrete  dipole  approximation  (Draine  and  Flatau  1994)  and  finite  difference 
time  domain  technique  (Liou  2002)  are  examples  of  numerical  techniques  now  possible. 

Mie  theory  begins  with  the  vector  wave  equation,  which  satisfies  Maxwell's  equations, 
and  governs  the  propagation  of  electromagnetic  waves.  The  wave  equation  is  transformed  into 
spherical  coordinates,  and  the  appropriate  boundary  conditions  for  spherical  geometry  are 
applied.  The  solution  for  the  resulting  fields  is  obtained  by  separating  variables  so  the  solution 
has  radial  and  zenith-angle  components.  Solutions  are  linear  combinations  of  a  series  of 
spherical  Bessel  functions,  an  and  bn,  and  angle-dependent  functions,  nn  and  rn.  nn  and  t„  are 
functions  of  the  scattering  angle,  0,  between  the  incident  and  scattered  direction.  an  and  bn  are 
complex  functions  of  the  particle's  refractive  index  and  size  parameter,  x.  x  is  non-dimensional 
and  expresses  the  particle's  radius,  a,  relative  to  the  wavelength  of  the  incident  radiation: 

2m 

X~~x 

By  summing  the  energy  entering  and  leaving  a  volume  surrounding  the  particle, 
expressions  can  be  obtained  for  the  extinction,  Ce,  and  scattering,  C5,  cross  sections  (units  of 
area)  of  the  particle.  The  corresponding  scattering  and  extinction  efficiencies,  Qs,  and  Qe,  are  the 
ratio  of  these  cross-sections  to  the  particle's  geometric  cross-sectional  area,  and  are  defined  in 
terms  of  the  spherical  Bessel  functions: 

e,=4E(2"+1)Rek+0 

x  n=\ 

a=4i(2»+i)kr+W!) 

x  n= 1 

The  single-scatter  albedo,  GJ,  is  the  ratio  of  the  scattering  efficiency  to  the  extinction 
efficiency.  The  asymmetry  parameter,  g,  is  the  expectation  value  of  cos(6).  Particles  that  scatter 
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isotropically  will  have  g  =  0  (however,  g  can  also  be  zero  even  for  non-isotropic  scattering); 
those  that  scatter  preferentially  in  the  forward  direction  will  have  values  of  g  between  0  and  1, 
while  those  that  scatter  more  radiation  backward  will  have  g  between  0  and  -1.  g  is  also  a  series 
solution  of  the  spherical  Bessel  functions  (Bohren  and  Huffman  1983): 


Collectively,  Qe,  G7,  and  g,  are  referred  to  as  a  particle's  Mie  properties,  as  they  describe 


how  much  incident  radiation  a  particle  extinguishes,  what  fraction  of  the  extinguished  radiation 
is  scattered,  and  the  degree  to  which  it  is  forward-  or  back-scattered.  Figures  3.1  and  3.2  depict 
Mie  properties  at  the  lowest  and  highest  AMSU-A  and  -B  frequencies  for  rain,  snow,  and  two 
different  densities  of  graupel. 


The  differing  behavior  of  the  curves  for 
different  hydrometeors  results  from  the 
marked  difference  in  refractive  index  between 
ice  and  water  in  the  microwave  portion  of  the 
spectrum.  Ice  has  a  very  small  imaginary 
component,  which  causes  its  Qe  and  Q,  to  be 
very  nearly  equal.  Thus,  ice  mostly  scatters 
microwave  radiation,  while  absorbing  very 
little  Water,  on  the  other  hand,  has  an 

imaginary  component  of  the  same  order  of 
magnitude  as  its  real  component,  so  it  both 
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Figure  3.1 :  Mie  properties  for  liquid  water, 
ice,  graupel,  and  melting  graupel  at  23.8 
GHz  (AMSU-A  channel  11. 
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Graupel  (50%  Hq.  35%  air,  15%  Ice) 
Graupel  (70%  air,  30%  ice) 


Figure  3.2:  Same  as  Figure  3.1,  but  for 
183.31  GHz  (AMSU-B  channels  18-20). 


absorbs  and  scatters  (57  <  1 ).  Graupel  is  a 
type  of  hydrometeor  composed  of  rimed  snow 
particles,  which  may  also  contain  water  if  in 
the  process  of  melting.  Determining  the 
complex  refractive  index  for  this  type  of 
hydrometeor  requires  a  dielectric  mixing 
formula,  such  as  that  of  Maxwell  and  Garnett 
or  Bruggeman  (used  in  this  study)  and  the 
relative  mass  fraction  of  each  substance  in  the 
particle.  Figure  3.3  depicts  the  complex 
refractive  index  across  the  microwave 


spectrum  for  pure  water,  pure  ice,  and  two  densities  of  graupel. 


Frequency  (GHz) 

- Liquid  real 

"■  Liquid  imaginary 

-  Ice  real 

---  Ice  imaginary 


Frequency  (GHz) 

* — •  Graupel  (50%  liq,  35%  air,  15%  Ice)  real 

- Graupel  (50%  air,  35%  liq,  1 5%  ice)  imaginary 

- Graupel  (70%  air,  30%  ice)  real 

-  *  -  Graupel  (70%  air,  30%  ice)  imaginary 


Figure  3.3:  Complex  indices  of  refraction  across  the  microwave 
spectrum  for  pure  liquid  (black  curves,  left),  pure  ice  (gray  curves, 
left),  melting  graupel  (black  curves,  right),  and  low-density  graupel 
(gray  curves,  right). 
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The  variation  in  each  extinction  curve  with  hydrometeor  size  falls  into  three  broad 
regimes.  For  small  particles  whose  size  parameter  is  much  less  than  unity,  Qe  increases 
approximately  in  proportion  to  x,  while  G7  increases  with  x3.  This  is  the  Rayleigh  regime.  In  this 
regime,  absorption  depends  no  on  particle  size,  but  only  on  the  total  mass  of  absorbers  present, 
as  is  the  case  with  gaseous  absorption.  Cloud  droplets  have  small  enough  size  parameters  in  the 
microwave  spectrum  that  scattering  is  negligible  and  absorption  is  approximately  proportional  to 
cloud  liquid  water  density.  For  small  precipitating  drops  and  crystals,  scattering  starts  to  become 
significant,  but  in  the  Rayleigh  regime,  scattering  is  symmetric  forward  and  backward  and  g  ~0. 

As  the  size  parameter  approaches  unity  (in  the  microwave  spectrum,  this  includes 
hydrometeors  with  radii  from  0.25  to  2  mm),  scattering  and  extinction  efficiencies  peak.  This  is 
the  Mie  regime,  which  applies  to  large  hydrometeors.  The  extinction  efficiency  in  this  regime  is 
actually  greater  than  1  (it  can  in  fact  reach  ~  4)  due  to  diffraction  of  radiation  missing,  but 
passing  close  by,  the  particle.  In  the  Mie  regime,  the  asymmetry  parameter  approaches  one  as 
scattering  becomes  sharply  forward-peaked. 

Very  large  size  parameters  (greater  than  approximately  104)  comprise  the  geometric 
optics  regime,  where  a  particle  is  so  large  compared  to  the  incident  wavelength  that  it  no  longer 
needs  to  be  treated  as  a  collection  of  individual  dipoles.  An  incident  beam  can  be  traced  through 
it  as  a  ray,  obeying  Snell's  laws  for  each  reflection  and  refraction  at  the  boundaries.  In  the 
geometric  optics  regime,  scattering  and  extinction  efficiencies  approach  and  oscillate  (due  to 
interference)  about  limiting  values.  Extinction  efficiencies  exceeding  unity  and  large  forward 
scattering  peaks  in  the  Mie  and  geometric  optics  regimes  are  caused  by  diffraction  slightly 
distorting  the  paths  of  photons  that  pass  near,  but  do  not  contact,  the  particle. 
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The  path  traced  out  by  an  electromagnetic  wave’s  electric  field  vector  as  it  oscillates 
defines  the  wave's  polarization.  Maxwell's  equations  are  satisfied  by  waves  with  linear  or 
elliptical  polarization.  A  beam  of  radiation  may  in  turn  be  composed  of  waves  with  a  preferred 
polarization  or  a  collection  of  waves  with  random  polarizations,  in  which  case  the  beam  is 
unpolarized.  The  polarization  state  can  be  described  by  the  Stokes  vector,  /,  which  is  composed 
of  4  scalar  elements,  7,  Q,  U,  and  V.  I  describes  the  total  intensity  of  the  beam;  Q  the  intensity 
having  vertical  ( Q  >  0 )  or  horizontal  (Q  <  0)  polarization  relative  to  the  plane  of  scattering;  U 
the  intensity  having  +/-  45°  oblique  polarization;  and  V  the  intensity  with  either  right-hand  or 
left-hand  circular  polarization.  Each  interaction  the  wave  undergoes  can  be  described  by  a 
Mueller  matrix,  M,  which  relates  the  scattered  Stokes  vector,  to  the  incident  Stokes  vector,  /: 

'/■I  rs„  sa  s„  s14y/' 

Q  _  1  *^21  ‘-*22  *^23  ^24  Q 

U'  s3l  s32  s33  s3A  u 

“^42  *^43  ‘^44_V^y 

In  the  special  case  of  a  spherical  particle,  the  Mueller  matrix  becomes: 

'/■'  \su(e)  sl2{0)  o  o  Y/> 

Q  =J_  S12{0)  Sn{0)  0  0  Q 

U'  k2r 2  0  0  S33(0)  S34(0)  U 

y'J  [  0  0  -s34(0)  s33{0)\[v; 

where  Sn,  Sn,  S33,  and  S34  are  all  functions  of  the  scattering  zenith  angle,  6,  the  spherical  Bessel 
functions,  an  and  bn,  and  the  angle-dependent  functions,  nn  and  rn.  It  is  important  to  note  here 
that  scattering-even  by  spherical  particles-can  cause  polarization  (the  polarization  of  sky  light 
is  an  example). 
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The  scattered  intensity  at  a  given  scattering  angle,  as  determined  by  the  Mueller  matrix,  is 
related  to  the  incident  intensity  by  the  phase  function,  P(6).  The  phase  function  is  in  turn 
determined  by  elements  Sn  and  Sn,  as  well  as  the  intensity  and  degree  of  horizontal/vertical 
polarization  of  the  incident  beam.  Since  the  computation  of  these  matrix  elements  requires  a 
series  solution,  computing  them  for  all  scattering  angles  can  be  computationally  expensive. 
Instead,  the  asymmetry  parameter  can  be  used  with  one  of  a  variety  of  approximate  phase 
functions  to  predict  scattered  intensity  vs.  scattering  angle.  The  most  common  of  these  is  the 
Henyey-Greenstein  phase  function: 


1-g2 

(l  +  g2  -2gcos((9))2 

Chapter  2  dealt  with  one-dimensional  radiative  transfer  in  a  non-scattering  medium.  For 
radiation  traveling  in  a  specified  direction  Si  in  a  medium  that  does  scatter  (so  that  the  extinction 
coefficient,  f}e,  equals  the  sum  of  the  absorption,  pa,  and  scattering,  jis,  coefficients),  we  need  to 
use  the  general  form  of  the  radiative  transfer  equation: 


Phg  {&) 


dl{Sl) 

ds 


=  - j3el{Sl)+j3aB(T )  +  A.  Jp(Q,  Sl')l{£l')dil' 


Using  the  definitions  of  the  single-scatter  albedo  and  optical  path,  this  reduces  to: 

=  -/(Q)+ (l  -  BT)B(D  +  —  f  P(il,  Sl')l(Sl')dSl’ 
dr  4  7C  l 

So  to  compute  the  radiative  properties  of  an  atmospheric  layer,  we  need  to  know  three  quantities: 
the  extinction  coefficient,  /?e;  the  single-scatter  albedo,  (U,  and  the  scattering  phase  function, 
P(£2,£2’),  which  relates  the  intensity  scattered  in  direction  Q  that  originated  from  direction  Si' . 
We  can  approximate  the  phase  function  using  g.  Multiplying  the  extinction  efficiency,  cross- 
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sectional  area,  and  number  density  of  a  given-size  hydrometeor  results  in  the  extinction 
coefficient  for  that  size  particle.  Summing  the  extinction  coefficients  for  all  sizes  of  the 
hydrometeor  size  distribution  yields  the  total  extinction  coefficient.  The  same  procedure  yields 
the  total  scattering  coefficient,  and  the  ratio  of  these  quantities  produces  the  single-scatter 
albedo.  The  asymmetry  parameter  for  a  collection  of  particles  is  the  sum  of  each  size’s 
asymmetry  parameter,  weighted  by  number  density. 

Since  actual  precipitation  size  distributions  can  very  widely,  we  arrive  at  the  number 
density  for  each  infinitesimal  increment  of  liquid  drop  or  ice  crystal  diameter,  D,  using 
empirically-derived  drop  size  distribution  models.  The  well-known  Marshall-Palmer  distribution 
for  rain  and  graupel  relates  rain  rate,  R,  to  number  density: 


N(D)  =  N0e~A(R)D  , 


N0  -  0.08cm-4 


A  (R)  -  41  cm" 


R 

mm  hr"1 


where,  in  this  case,  R  is  the  rain-equivalent  precipitation  rate.  Petty  (2001a)  introduced  a 
modification  to  both  distributions  by  dividing  the  exponential  factor  by  25.  This  enables  the 
effective  particle  diameter  to  be  scaled  up  or  down  readily  in  order  to  produce  more  realistic 
extinction  and  scattering  coefficients.  In  this  study,  as  in  Petty  (2001a),  3=2  is  employed  for  the 
snow  size  distribution.  The  rain  rate  can  be  derived  from  a  given  mass  density  of  hydrometeors 
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Figure  3.4:  Extinction  coefficient  (km'1;  left),  single-scatter  albedo  (center),  and 
asymmetry  parameter  (right)  obtained  from  Mie  calculations  at  AMSU-A  and  -B 
frequencies  for  various  rain  densities.  Particle  size  distributions  used  in  the  calculations 
were  obtained  from  a  Marshall-Palmer  distribution  modified  using  the  method  of  Petty 
(2001)  to  account  for  fall-speed  variation  with  air  density. 

once  the  fall  speed,  u  (a  function  of  hydrometeor  size),  is  known.  Usually  a  power-law  relation, 

u(D)  =  aDy,  is  assumed. 

Figures  3.4  through  3.6  depict  the  extinction  coefficient,  single-scatter  albedo,  and 


Precipitation  density  (10M  kg  m^3)  Precipitation  density  (10M  kg  m*- 3)  Precipitation  density  (lOM  kg  m*- 3) 


Figure  3.5:  Same  as  Figure  3.4,  but  for  snow  (modeled  as  ice  spheres).  Particle  size 
distributions  derived  from  Sekhon-Srivastava  size  distribution,  with  reverse-exponential 
decay  parameter  scaling  and  fall-speed  modifications  from  Petty  (2001). 
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asymmetry  parameter  vs.  precipitation  density  for  rain,  snow,  and  graupel  at  AMSU-A  and  -B 
frequencies.  From  these  figures,  we  can  see  that  rain  extinction  efficiency  increases  with 
frequency.  Its  single-scatter  albedo  does  also,  making  it  a  moderate  scatterer  at  high  frequencies. 
At  those  frequencies,  it  scatters  strongly  in  the  forward  direction  ( g  ~0.5);  at  temperature 
sounding  frequencies  it  scatters  isotropically,  and  at  low  frequencies  it  actually  back-scatters 
slightly. 

Snow  extinction  is  much  smaller  than  rain  at  all  frequencies,  and  only  becomes 
significant  at  high  frequencies.  Single-scatter  albedo  is  larger  than  for  rain,  due  to  ice's  very 
small  imaginary  index  of  refraction.  Snow  asymmetry  parameters  are  between  0  and  0.5,  so 
when  snow  scatters,  it  does  so  with  a  moderately-strong  phase  function  forward  peak. 

Graupel's  low  density  makes  it  considerably  larger  (nearly  a  factor  of  2)  for  the  same 
precipitation  density.  The  larger  geometric  size  increases  its  extinction  coefficient  relative  to 
rain  and,  especially,  snow.  The  larger  particles  scatter  more  strongly,  and  with  a  sharper  forward 
peak  (g  reaching  0.8  for  high  frequencies). 


Graupel 

*  •  •  183.31  GHz 
XXX  89.0  GHz 
5-.  +++  54  94  GHz 
□□□  50.3  GHz 
O  31.4  GHz 
4"  ooo  23.8  GHz 


m  ffl]S  A  |fi  S  t  §  a | 

0  1  2  3  4  5 


Precipitation  density  (10M  kg  m'VS) 


Figure  3.6:  Same  as  Figures  3.4  and  3.5,  but  for  graupel  (70%  air,  30%  ice).  Fall-speed 
modified  Marshall-Palmer  size  distribution  used  for  Mie  calculations. 
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3.2.  Previous  Microwave  Precipitation  Studies 

To  date,  the  bulk  of  existing  microwave  precipitation  research  has  been  focused  on 
quantifying  surface  precipitation  rates  or  inferring  precipitation  profiles  using  microwave 
instruments. 

The  earliest  work  in  this  area  occurred  in  the  mid-1970's,  after  the  first  research 
microwave  radiometers  were  flown.  Grody  (1976)  noted  that  a  pair  of  microwave  frequencies 
located  different  distances  from  the  22  GHz  WV  absorption  band  could  be  used  to  determine  the 
optical  depth  of  water  vapor  and  cloud  liquid  water  present.  He  presented  regression-based 
algorithms  to  retrieve  total  precipitable  water  and  cloud  liquid  water  using  the  Nimbus-E 
Microwave  Spectrometer's  22.235  GHz  and  31.40  GHz  channels. 

Wilheit  et  al.  (1977)  employed  19.35  GHz  TB's  from  the  Electrically  Scanning 
Microwave  Radiometer  (EMSR)  to  estimate  surface  rain  rate.  He  used  Mie  theory  as  described 
above  to  model  TB19  of  a  raining  oceanic  cloud  of  varying  thicknesses  and  with  varying  freezing 
levels,  both  with  and  without  scattering.  Using  this  information,  he  developed  a  relation  between 
TB19  and  over-ocean  surface  rain  rates  (up  to  20  mm  hr'1)  that  was  accurate  to  within  a  factor  of 
two.  A  number  of  variations  of  this  emission-based  rain  rate  estimation  technique  have  been 
developed  for  succeeding  instruments. 

Wu  and  Weinman  in  1984  used  the  Nimbus-7  Scanning  Multi-channel  Microwave 
Radiometer  to  study  ice  scattering  effects  on  37  GHz  TB  and  concluded  that  at  this  frequency 
and  above,  scattering  effects  become  significant  and  could  be  used  to  infer  surface  rain  rates. 
Spencer  (1986)  expanded  on  this,  developing  a  scattering-based  rain  rate  estimate  which 
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compares  dual-polarized  TB37  to  expected  values  from  a  non-raining  cloudy  or  clear  oceanic 
scene. 

The  introduction  of  the  Special  Sensor  Microwave  Imager  (SSMI),  flown  continuously 
by  the  Defense  Meteorological  Satellite  Program  (DMSP)  constellation  since  1987,  enabled  a 
great  deal  of  groundbreaking  work  in  precipitation  measurement.  SSMI  is  a  conically-scanning 
microwave  radiometer  with  channels  at  19,  22,  37,  and  85  GHz.  All  channels  are  dual-polarized, 
except  22  GHz,  which  is  vertically-polarized  only.  Instrument  resolution  varies  from  70  km  at 
19  GHz  to  16  km  at  85  GHz.  SSMI  was  the  first  instrument  to  employ  both  low-  and  high- 
frequency  window  frequencies  simultaneously,  and  led  to  a  number  of  novel  algorithms  for 
retrieving  water  vapor,  cloud  liquid,  precipitation,  and  ocean  and  land  surface  parameters. 

Grody  (1991)  developed  an  empirical  scattering  index,  which  uses  TB19V  and  TB22V  to 
predict  what  TB85V  would  be  in  the  absence  of  scattering  hydrometeors.  He  used  the  difference 
between  predicted  and  actual  TB85V  to  classify  precipitating  scenes,  but  not  to  estimate  rain  rate. 
Liu  and  Curry  (1992)  followed  a  similar  approach,  but  used  both  low-frequency  (19  GHz)  and 
high-frequency  (85  GHz)  TB  differences  as  rain  rate  predictors.  In  his  SSMI  rain  rate  algorithm. 
Petty  (1994)  used  an  85  GHz  scattering  index,  S,  to  determine  a  first-guess  rain  rate,  then 
iteratively  adjusted  it  until  the  predicted  19  and  37  GHz  normalized  polarization  index,  P,  which 
measures  rain  cloud  opacity,  matches  the  observed  values. 

Bennartz  and  Petty  (2001)  sounded  a  cautionary  note  about  scattering-based  precipitation 
retrievals.  They  used  SSMI  TB' s,  coupled  with  varying  ice  particle  size  distributions,  to  predict 
radar  reflectivity,  which  they  compared  with  actual  radar  observations.  The  results  indicate  a 
strong  dependence  of  TB  on  the  assumed  particle  size  distribution  used  in  their  Mie  calculations. 
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The  Tropical  Rainfall  Measuring  Mission  (TRMM)  was  launched  in  1997,  and  continues 
to  operate  today.  TRMM  contains  a  microwave  imager  (TMI)  nearly  identical  to  SSMI,  as  well 
as  an  active  microwave  instrument,  the  precipitation  radar  (PR).  Bauer  (2001)  used  the  new 
instrument,  TMI,  as  well  as  a  new  approach  to  measuring  precipitation:  principal  component 
analysis  (PCA).  PCA  identifies  patterns  of  covariance  in  multivariate  datasets  by  computing  the 
eigenvalues  and  eigenvectors  of  the  data's  covariance  matrix.  The  eigenvectors  span  the 
observation  space,  and  so  each  multi-variable  observation  can  be  reproduced  by  a  linear 
combination  of  the  eigenvectors.  The  coefficient  corresponding  to  each  eigenvector  is  its 
principal  component  (PC).  The" eigenvalue  corresponding  to  each  eigenvector  represents  the 
fraction  of  the  data's  variance  explained  by  that  eigenvector.  Bauer  related  PC's,  rather  than  a 
single  channel  or  subset  of  channels,  to  rain  rate,  and— surprisingly— found  that  scattering  effects 
on  TB85  reduced  its  usefulness  in  predicting  rain  rate;  excluding  TB85  from  the  data  set  prior  to 
PCA  improved  the  correlation  of  the  PC's  with  rain  rate. 

Petty  (2001a)  took  a  similar  approach,  using  PCA  with  SSMI  TB's  to  analyze  tropical 
stratiform  precipitation.  He  found  that  the  first  EOF,  whose  elements  were  all  the  same  sign, 
corresponded  to  overall  ice  scattering  intensity.  The  second  EOF,  with  TB19  and  TB22 
deviations  of  opposite  sign  from  TB37  and  TB85,  represented  the  multi-channel  rain  signature 
(i.e.  low-frequency  rain  emission  and  high-frequency  ice  scattering).  The  third  EOF,  dominated 
by  TB85  deviation,  accounted  for  anomalously  strong  85  GHz  ice  scattering.  Interestingly,  the 
magnitudes  of  the  components  of  EOF  1  follow  a  power-law  relationship,  allowing  a  frequency- 
dependent  prediction  of  scattering  attenuation. 
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A  fourth  form  of  rain  retrieval  capitalizes  on  recent  strides  in  fine-scale  numerical 
weather  prediction  (NWP)  model  resolution  and  microphysics.  An  NWP  model  is  used  to 
produce  temperature,  humidity,  and  precipitation  fields  at  cloud  or  sub-cloud  scales.  The 
meteorological  profile  at  each  horizontal  grid  point  is  then  used  as  input  to  a  radiative  transfer 
model,  typically  a  1 -dimensional  plane-parallel  Eddington  two-stream  approximation,  to 
compute  upwelling  microwave  brightness  temperatures.  The  grid-scale  brightness  temperature 
field  is  then  convolved  to  the  resolution  of  the  instrument  being  simulated.  The  simulation 
results  are  stored,  creating  a  library  of  brightness  temperatures  matching  different  precipitation 
profiles.  These  libraries  form  the  basis  for  Bayesian  precipitation  profile  retrievals,  which  infer 
the  profiles  from  the  library  profiles  best  matching  a  set  of  observed  brightness  temperatures 
(Kummerow  1996).  An  extensive  series  of  papers  spanning  the  1990's  and  first  half  of  this 
decade  have  focused  on  optimal  techniques  for  creating  these  "cloud-radiation  databases,"  and 
improving  NWP  model  microphysical  assumptions  to  better  match  observed  brightness 
temperatures  (Mugnai  et  al.  1990;  Smith  et  al.  1991;  Mugnai  et  al.  1993;  Smith  et  al.  1994; 
Panegrossi  et  al.  1998). 

In  the  same  vein,  Panegrossi  (2004)  focused  attention  on  mismatches  between  observed 
Tropical  Rainfall  Measuring  Mission  (TRMM)  Microwave  Imager  (TMI)  TB's  from  Hurricane 
Bonnie  (1998)  and  simulated  TB's  from  a  University  of  Wisconsin  Nonhydrostatic  Modeling 
System  (UW-NMS)  simulation  of  the  storm.  She  demonstrated  that  frozen  hydrometeors  can 
have  strong  scattering  effects  on  37  GHz  brightness  temperatures,  and  reduced  simulated  vs. 
observed  discrepancies  at  this  frequency  by  changing  the  model's  hydrometeor  categorization  of 
low-density  graupel  (essentially  rimed  ice  crystals)  to  snow,  leaving  the  graupel  category  for 
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hail.  She  also  showed  that  anomalously  warm  85  GHz  brightness  temperatures  could  be  caused 
by  emission  from  supercooled  water  droplets  lofted  into  the  upper  troposphere  by  deep  cumulus 
convection  and  warm  10  GHz  brightness  temperatures  by  overly  strong  rain. 

Just  as  SSMI  and  the  TMI/PR  combination  each  triggered  a  wave  of  rain  rate  and 
precipitation  profile  retrievals,  the  1998  introduction  of  AMSU-A  and  -B,  with  a  suite  of  window 
channels  well-suited  for  cloud  water  and  precipitation  profiling,  spawned  a  large  body  of  work 
on  microwave  precipitation  effects.  Grody  (1999)  adapted  his  SSMI  scattering  index  for  AMSU- 
A  use,  using  regression  analysis  to  develop  an  expression  for  non-scattering  TB89  using  TB23 
and  TB31.  Grody  then  used  an  observed  TB89  more  than  9  K  colder  than  the  predicted  value  to 
distinguish  areas  of  precipitation  over  both  land  and  sea. 

Ferraro  et  al.  (2000)  describe  the  algorithms  developed  for  operational  use  at  NESDIS. 
They  quantify  scattering  by  simply  subtracting  TB89  from  TB23.  Scattering  indices  >  3  K  are 
associated  with  stratiform  rain;  >  40  K  with  cumulonimbus  convection.  For  each  precipitation 
type,  Ferraro  et  al.  derive  rain  rate  using  a  power  law  with  different  coefficients  for  each  type  of 
precipitation.  A  similar  algorithm,  with  the  benefit  of  AMSU-B's  finer  spatial  resolution,  uses 
TB89  -  TB150  as  the  scattering  index.  Bennartz  et  al.  (2002)  took  a  similar  approach,  using 
AMSU-A  TB89  or  AMSU-B  TB150  vs.  AMSU-A  TB23  to  derive  categorical  mid-latitude 
precipitation.  Weng  et  al.  (2003)  describe  NESDIS'  operational  total  precipitable  water  (TPW), 
and  cloud  liquid  water  algorithms,  which  capitalize  on  the  larger  continuum  water  vapor 
absorption  with  increasing  frequency.  Finally,  Bennartz  and  Bauer  (2003)  compared  the  ice 
scattering  signature  at  85,  150,  and  183  GHz  to  each  frequency's  sensitivity  to  surface  and  water 
vapor  emission.  They  found  150  GHz  to  be  the  optimal  frequency  for  characterizing  scattering. 
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It  is  less  sensitive  to  surface  emission  than  85  GHz  and,  while  183  GHz  channels  are  very 
sensitive  to  scattering  (Bums  et  al.  1997),  150  GHz  is  less  sensitive  than  183  GHz  to  emission  by 
variable  water  vapor  content. 

While  the  need  for  accurate  precipitation  rate  and  profile  retrievals  has  motivated  the 
large  body  of  work  studying  precipitation  effects  on  window  frequencies,  comparatively  little 
has  been  done  to  study  how  precipitation  affects  AMSU-A  temperature  sounding  channels  ( TB4 
through  TB14\  52.8  to  57.29  GHz).  Reale  (2001)  and  Greenwald  et  al.  (2004)  recognize  that 
scattering  can  affect  sounding  channel  TB's  and  point  out  that  current  global  NWP  model  data 
assimilation  schemes  reject  radiances  from  precipitating  AMSU-A  fields  of  view.  Such  fields  of 
view  are  identified  by  cloud  liquid  water  content  >  0.2  mm.  Li  and  Weng  (2002)  found  that 
scattering  TB  depression  in  the  lowest  tropospheric  sounding  channels  (AMSU-A  channels  4  and 
5,  with  weighting  functions  that  peak  at  the  surface  and  near  700  hPa  respectively)  is  well- 
correlated  with  ice  water  path  (IWP)  and  somewhat  correlated  with  liquid  water  path  (LWP). 

A  known  temperature  profile  is  a  necessary  precursor  to  retrieving  moisture  profiles 
using  the  183  GHz  AMSU-B  channels  ( TB18  through  TB20),  since  the  water  vapor  content 
dictates  what  level  each  channel  peaks  at,  but  the  temperature  near  that  level  dictates  what  the 
actual  183  GHz  TB  will  be.  As  a  first  step,  then,  in  their  humidity  profile  retrieval  algorithm, 
Chen  and  Staelin  (2003)  required  temperatures  from  AMSU-A  channels  4  through  8. 
Recognizing  that  channels  4  and  5  could  be  strongly  affected  (several  K)  by  precipitation 
attenuation,  they  devised  a  scheme  to  predict  TB4  and  TB5  precipitation-related  perturbations  at 
AMSU-B  resolution  using  the  strong  scattering  sensitivity  of  AMSU-B  TB19  and  TB20.  They 
then  corrected  TB4  and  TB5  using  a  regression  relation  with  the  TB19  and  TB20  perturbations. 
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Most  recently,  Bauer  and  Mugnai  (2003)  capitalized  on  precipitation  attenuation  of 
sounding  channels  to  develop  a  precipitation  profile  algorithm  employing  60  and  1 1 8  GHz  TB 
from  an  aircraft-mounted  microwave  radiometer.  Their  Bayesian  approach  uses  a  profile  library 
developed  from  a  numerical  simulation  of  Hurricane  Bonnie  (1998)  and  60  and  118  GHz  TB' s 
predicted  using  a  two-stream  Eddington  radiative  transfer  model.  The  retrieval  is  based  on  the 
differential  scattering  susceptibility  of  60  and  118  GHz  channels  that  peak  at  nearly  the  same 
level. 

3.3.  A  1-d  Conceptual  Model  for  Precipitation  Effects  on  AMSU  Brightness  Temperatures 

As  we  have  seen,  theory  and  previous  studies  demonstrate  the  strong  susceptibility  of 
high-frequency  window  channels  to  ice  scattering.  Chapters  4  and  5  will  use  radiative  transfer 
simulations  and  a  sample  of  observational  AMSU  data  to  test  the  hypothesis  that  scattering 
depression  of  high-frequency  window  channels  (TB15  through  TB20)  can  be  used  to  quantify  the 
scattering  depression  of  55  GHz  temperature  sounding  channels  (TB5  through  TB8).  Before 
proceeding  to  those  steps,  this  section  presents  a  simple  1-D  conceptual  model  to  check  the 
feasibility  of  such  an  approach. 

Figure  3.7  depicts  the  one-dimensional  plane-parallel  model  geometry.  The  ocean 
surface  is  specular,  with  reflectivity  rsea  and  temperature  Tsea.  The  troposphere  is  divided  into 
two  layers:  the  upper  troposphere,  which  is  the  area  of  .primary  sounding  channel  sensitivity,  and 
the  lower  troposphere,  which  lies  below  the  bulk  of  the  sounding  channel  weighting  function. 

The  lower  tropospheric  layer  has  temperature  7)a  and  transmittance  tia.  Superimposed  on  the 
lower  tropospheric  layer  is  a  rain  layer  with  temperature  Train,  transmittance  tmin,  and  single- 
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Figure  3.7:  Geometry  of  the  1-D  plane-parallel  conceptual  model  for  hydrometeor  effects 
on  AMSU  window  and  sounding  channel  TB's.  The  brightness  temperature  observed  by 
the  instrument  is  the  sum  of  12  components  (see  text).  Solid  lines  denote  direct  emission 
paths  while  dashed  lines  denote  reflected  emission  paths.  T  indicates  a  layer's  temperature, 
t  its  transmittance,  and  GJ  its  single-scattering  albedo. 

scatter  albedo  GJrain.  Capping  the  rain  layer  is  an  ice  layer  with  temperature  Tice,  transmittance 
tice,  and  single-scatter  albedo  GJice.  Above  the  cloud  is  an  upper-atmosphere  layer  with  Tua  and 
tua.  The  cosmic  background  temperature,  Tcb,  is  neglected.  Two  radiance  streams— one  upward 
and  one  downward— are  tracked.  All  scattering  directs  intensity  from  one  stream  into  the  other. 
Multiple  scattering  is  neglected. 

The  predicted  TB  measured  by  the  instrument  is  composed  of  twelve  components: 

1.  Direct  upward  emission  from  the  upper-tropospheric  layer: 


tb,  =r_(i-»J 


(Note  that  for  a  non-scattering  layer,  such  as  the  upper  troposphere,  the  sum  of  absorption  and 
scattering  is  unity.  Kirchhoff  s  Law  requires  emissivity  to  equal  absorptivity,  so  the  emissivity 
of  a  non-scattering  layer  becomes  (1-t).) 

2.  Downward  emission  from  the  upper-troposphere  reflected  by  the  ice  layer  and  transmitted 
back  upward  through  the  upper-troposphere: 

TB,=T  (l-f  W.  (l -f.  V 

2  ua\  uar^ice\  ice)  ua 

3.  Downward  upper-tropospheric  emission  reflected  by  the  rain  layer  and  transmitted  back 
upward  to  space: 

TB3  =  Tua  (l  -  ha  he® rain  C1  “  train  heha 

4.  Downward  upper-tropospheric  emission  reflected  by  the  ocean  surface  and  transmitted  back 
upward: 

BB4  Bua  (l  ha^icehau)laheahahainhceha 

5.  Upward  ice-layer  emission  transmitted  to  space: 

TB5=Tice{l-tice){l-CUiceha 

(For  the  ice  and  rain  layers,  where  scattering  occurs,  the  layer's  emissivity  is  (l-t)(l-G7).) 

6.  Downward  ice  layer  emission  reflected  off  the  rain  layer  and  transmitted  upward  to  space: 

TB6  ~  Tice  (l  ~  tice  Xl  ice  rain  (l  “ 1  rain  ^iceha 

1.  Downward  ice  layer  emission  reflected  off  the  ocean  surface  and  transmitted  back  upward: 

TB1  =  Tice  (l  -  tice  Xl  ~®ice  }  raJiaheatJraJ.Jua 

8.  Upward  emission  from  the  rain  layer  transmitted  through  the  ice  layer  and  upper-troposphere: 
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9.  Downward  emission  from  the  rain  layer  reflected  by  the  ocean  surface  and  transmitted  back 
upward: 


TB9  =  Train  0  “  * rainfo  ~  ® rain 

10.  Upward  emission  from  the  lower  troposphere  transmitted  through  the  rain,  ice,  and  upper- 
troposphere  layers: 

TB io  =  TIa  (l  —  tla  )trainticetua 

1 1 .  Downward  emission  from  the  lower  troposphere  reflected  by  the  ocean  surface,  then 
transmitted  upward  through  all  layers: 

TBl  1  =  ' Tla  (!  ■ -  ha  VseahaKamhJua 

12.  Upward  emission  from  the  ocean  surface,  transmitted  upward  to  space: 

TBn  =  TSea  (1  ■ -  hea  XaKaJiJua 

(For  the  ocean  surface,  reflectivity  and  absorptivity  sum  to  unity  so  its  emissivity  is  (7  -  r).) 
The  total  brightness  temperature  measured  from  space  is  then  the  sum  of  these  twelve 


components: 


TB  = 


Tjl-tJ 


+  GJice(l-tice)tua 
+  ®rainb-train)t?Jua 
+  r  t2tz  t2 1 

'  *  sea*  la1  rain1  ice* ua 


+  Tice  (1  -  hce  Xl  -  ® ice  )|  +  ®  rain  (*  "  Kain  keh 

h 

1 


rain  \  rain  r  ice  ua 

+  r  tft2- 1 .  t 

sea  la  rain  ice  ua 


+  Train  “  train  X1  “  ® rain  keKa 


+  r  t,t  ■ 

sea  la  rain 


+  T,  i\-t,  )t  t-  t 

la\  la)  rain  ice  i 

+  T  (l  —  r  )t,  t  .  t-  t 

sea  v  sea  )  la  rain* ice  i 


1 


+  hea1  la 


(3.1) 
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Limiting  cases  of  clear  conditions,  an  opaque  lower  troposphere,  an  opaque  rain  layer, 
and  an  opaque  ice  layer  highlight  the  effects  of  water  vapor  and  precipitation  at  different 
frequencies.  For  the  limiting  case  of  an  atmosphere  with  no  precipitation  (tra,„  =  tice  =  1 ): 

TBnon-precip  =  +  Tla  (*  ~  4,  ha  l1  +  T sea1  la  ) 

+  Tsea(l-rseahtua 

For  window  channels,  where  tm  ~tia~l,  TB  -4  Tsea  (l  -  rsea )  and  the  observed  TB  will  be 

approximately  equal  to  that  of  the  radiometrically  cold  ocean  surface.  For  temperature  and 
moisture  sounding  channels,  where  tia  ~ 0  and  tua  is  small,  TB  ~^Tm( 1  -  tm )  +  TIa tm  and  the 

observed  TB  will  be  nearly  equal  to  the  upper-tropospheric  temperature,  with  a  small 
contribution  by  emission  from  the  lower  troposphere.  As  lower-tropospheric  humidity  increases, 
and  tia  goes  to  zero  (faster  at  high-frequencies)  due  to  WV  absorption,  TB  — >  Tla  and  the 

observed  TB  approaches  the  opaque  lower  troposphere's  temperature. 

For  an  atmosphere  with  an  opaque  rain  layer  and  no  ice: 


TB 


opaque  _  warm  _  rain 


Tjl-tJ{l+GTraintua) 

+  Tram(l-G7rain)tua 


In  this  case,  observed  TB  depends  only  on  upper-troposphere  transmittance  and  the  rain  layer's 
single-scatter  albedo.  The  Mie  calculations  in  Sec  3.1  show  that  Wram  is  about  0.25  for  low- 
frequency  window  frequencies,  and  increases  to  ~  0.5  at  high-frequencies.  tua  will  be  close  to  1 
for  window  frequencies  and  TB  — >  Train  (l  -  G7rain ) .  So  the  opaque  rain  is  radiometrically  cold  for 


high  frequencies  due  to  scattering,  which  limits  its  emission,  while  it  appears  warmer  at  low 
frequencies  where  it  emits  more  effectively.  The  low  tua  limits  the  rain  layer's  effect  on  55  GHz 
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sounding  channel  TB.  In  effect,  only  the  bottom  tail  of  the  55  GHz  weighting  function  extends 
down  into  the  rain  layer. 

For  an  atmosphere  with  an  opaque  ice  layer: 

opaque _ice  =  Tua  0-  “  fua  X1  +  ^weKa  )  +  Tice  ® ice  Ka 

At  low  frequencies,  ice  extinction  is  very  weak  and  it  would  be  extremely  difficult  for  an  ice 
layer  to  be  opaque.  At  high  frequencies,  ice  can  be  significantly  opaque  and  its  higher  GJ  makes 
it  appear  even  colder  at  these  frequencies  than  the  opaque  rain  layer.  In  fact,  if  the  upper- 
troposphere  is  completely  transparent  and  the  ice  layer  is  completely  opaque,  TB  will  approach  0. 

The  1-D  plane-parallel  geometry  greatly  exaggerates  precipitation  effects  in  comparison 
with  real-world  observed  TB' s,  because  the  full  instrument  field  of  view  is  never  completely 
filled  with  uniform  deep  convection  producing  opaque  rain  and  ice  layers  like  the  ones  modeled 
here.  Instead,  localized  areas  affected  by  rain  and  ice  emission  and  scattering  are  averaged  with 
large  areas  of  ambient  TB,  muting  the  effect  of  the  precipitation.  Precipitating  clouds  should 
appear  warmer  when  modeled  in  3-D  for  two  reasons:  Petty  (1994)  pointed  out  that  when 
viewed  at  oblique  angles,  the  warm  rain  of  a  convective  cloud  can  be  viewed  through  the  side  of 
the  cloud  rather  than  through  the  ice  layer  capping  it,  reducing  scattering  attenuation. 
Additionally,  ocean  surface  reflection  of  the  warm  rain  emission  can  occur  for  some  geometries, 
warming  the  radiometrically  cold  ocean  surface’s  TB  and  counteracting  some  of  the  TB  reduction 
caused  by  ice  scattering  in  the  scene. 

Table  3.1  summarizes  modeled  low-frequency  window  ( TB31 ),  temperature  sounding 
( TB55 ),  and  high-frequency  window  ( TB89  and  TB150)  brightness  temperatures  using  equation 
(3.1)  for  ten  different  cases.  Rain  and  ice  layer  transmittances  and  single-scatter  albedos  are 
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TB31  OO  TB31  (K) 

♦  ♦  ♦  With  precip  '  *  ♦  ♦  With  precip 

— -  Without  precip  ?- —  Without  precip 


Figure  3.8:  Modeled  TB89  and  TB150  vs.  TB31.  Dashed  curve  depicts  predicted  TB  without 
precipitation  emission  or  scattering.  Solid  curve  depicts  modeled  TB  including  precipitation 
effects.  Table  3.1  summarizes  assumptions  and  results  for  each  of  the  10  cases  depicted  here. 

inferred  from  the  precipitation  Mie  properties  presented  in  Section  3.1.  The  rain  and  ice 

transmittances  for  higher  frequencies  are  defined  in  terms  of  the  31  GHz  transmittance,  using  the 

ratio  of  each  frequency's  extinction  coefficient  to  its  31  GHz  counterpart.  The  temperature  of 

each  layer  is  fixed.  For  window  frequencies  ( TB31 ,  TB89,  and  TB150 )  the  upper  troposphere 

layer  is  transparent,  while  for  the  temperature  sounding  channel  ( TB55 )  its  transmissivity  is  0.1. 

Sea  surface  reflectivity  is  set  to  0.55  for  TB31,  decreasing  to  0.35  for  TB150. 

Three  cases  contained  no  precipitation,  but  steadily  increasing  humidity,  represented  by 
decreasing  window  channel  tta.  This  increased  TB31  slightly  (135  to  187  K),  while  increasing 
TB89  (180  to  254  K)  and  TB150  (195  to  270  K)  more  significantly. 

Adding  a  warm  rain  layer  of  progressively  greater  opacity  continued  to  warm  TB31 
slightly  (195  to  209  K)  while  decreasing  both  the  high  frequency  window  channel  TB's  due  to 
scattering  ( TB89  drops  from  254  to  169  K  and  TB150  drops  from  270  to  160  K).  TB55  also 
begins  to  decrease  (by  as  much  as  1.5  K)  as  rain  intensity  increases,  due  to  the  rain  layer 
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replacing  warm  lower-tropospheric  emission  with  radiometrically  colder  rain.  Finally,  adding  an 
ice  layer  atop  the  warm  rain  layer,  and  increasing  its  opacity,  tremendously  reduced  the  high- 
frequency  window  TB's  (to  144  K  for  TB89  and  to  99  K  for  TB150 ),  and  continued  reducing 
TB55  (for  a  total  reduction  of  about  2.0  K). 

TB31  exhibits  very  little  effect  from  either  rain  or  ice,  so  TB89  and  TB150  scattering 
effects  can  be  visualized  by  plotting  them  against  TB31.  Figure  3.8  depicts  TB89  and  TB150  vs. 
TB31  for  non-precipitating  cases  (dashed  line)  and  precipitating  cases  (points).  Figure  3.9  in  turn 
depicts  the  strong  correlation  between  sounding  channel  brightness  temperature  reduction  due  to 
precipitation,  ATB55,  and  window-channel  TB  precipitation  reduction. 


Model  Sounding  vs:  Window  TB  Depression 
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Figure  3.9:  Modeled  TB55  depression  vs.  TB89  and  TB150  depression  for  the  10  cases 
summarized  in  Table  3.1. 


3.4.  Summary 

This  chapter  has  examined  precipitation  effects  throughout  the  microwave  spectrum,  with 
the  goal  of  determining  which  other  AMSU-A  or  -B  frequencies  may  be  useful  for  quantifying 
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precipitation  effects  on  the  AMSU-A  frequencies  near  55  GHz  used  to  measure  upper- 
tropospheric  temperature.  In  summary,  we  have  learned  that: 

•  Humidity  increases  window  channel  opacity,  especially  for  AMSU-A  TB1  and  TB15 
through  TB20,  where  humid  atmospheres  will  be  opaque 

•  Cloud  droplet  scattering  is  negligible  and  absorption  is  proportional  to  cloud  liquid  water. 
Cirrus  ice  crystals  are  small  enough  that  they  absorb  and  scatter  negligibly  in  the 
microwave  spectrum. 

•  Rain  emits  at  low  frequencies  (TB1  and  TB2)  and  scatters  at  high  frequencies  (TB  J5 
through  TB20) 

•  Snow  does  not  significantly  affect  TB1  or  TB2  unless  it  is  very  thick,  in  which  case  it  will 
scatter  weakly.  Snow  can  scatter  TB15  through  TB20  significantly. 

•  Graupel  will  scatter  strongly  at  all  frequencies,  but  especially  for  TB15  through  TB20. 

•  The  atmosphere  is,  in  general,  mostly  transparent  for  window  frequencies  ( TB1  and  TB2 
and  TB15  through  TB17 ),  so  they  are  sensitive  to  hydrometeors  throughout  the  vertical 
column.  Thus,  ATB  for  these  frequencies  is  closely  related  to  the  column-integrated 
liquid  and  ice  water  path.  For  temperature  sounding  ( TB5  through  TBS)  or  moisture 
sounding  ( TB18  through  TB20)  channels,  where  the  atmosphere  is  opaque  due  to  O2  and 
water  vapor  absorption,  respectively,  ATB  is  only  sensitive  to  hydrometeors  that  impinge 
significantly  on  the  weighting  function  for  that  channel. 

The  goal  for  Chapters  4  and  5  is  to  develop  an  empirical  scattering  correction  to  AMSU- 
A  sounding  channel  TB's  based  on  the  observed  depression  of  high-frequency  window  channel 
TB  s.  A  low-frequency  window  channel,  which  is  only  weakly  affected  by  scattering,  will  be 
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used  to  predict  the  non-scattering  high-frequency  window  channel  TB.  Principal  component 
analysis  will  also  be  explored  as  a  multi-channel  tool  for  predicting  sounding  channel  effects. 
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4.  Simulations  of  Precipitation  Effects  on  AMSU-A  and  -B  Channels 

Chapter  2  noted  that  the  current  CIMSS  TC  intensity  algorithm  infers  MSLP  from  the 
warm  anomaly  in  TB7  and  TB8.  The  goal  of  this  and  the  succeeding  chapter  is  to  develop  a 
correction  for  precipitation  effects  on  these  channels.  The  conceptual  model  in  Chapter  3 
showed  how  deep  rain  layers  impinging  on  these  channels'  weighting  functions  can  cool  TB  by 
replacing  radiometrically  warm  emission  from  the  humid  lower  tropospheric  with 
radiometrically  colder  (owing  to  its  ~  0.3  single-scatter  albedo)  rain  emission  and  reflected  cold 
upper-tropospheric  TB' s.  Frozen  precipitation  particles  in  cumulus  towers  will  affect  TB7  and 
TB8  more  strongly  than  rain,  by  virtue  of  their  greater  proclivity  to  scatter  (GJ~  0.9  for  graupel 
and  ~  0.5  for  snow). 

Previous  studies  of  precipitation  effects  on  microwave  radiation  have  shown  the  utility  of 
high-frequency  window  frequencies  (89  GHz  and  150  GHz)  for  quantifying  scattering  (Grody 
1991;  Petty  1994;  Grody  1999).  The  conceptual  model  from  Chapter  3  also  showed  that  the 
scattering-induced  TB  depression  at  these  frequencies  correlates  with  the  depression  of  sounding 
channel  TB's,  despite  the  fact  that  window  channels  respond  to  column  integrated  liquid  and  ice 
quantities,  while  sounding  channels  are  only  affected  by  liquid  and  ice  within  the  influence  of 
their  weighting  functions.  The  purpose  of  this  chapter  is  to  verify  that  this  relationship  holds 
under  more  realistic  atmospheric  profiles  of  humidity,  cloud,  rain,  snow,  and  graupel,  at  actual 
AMSU-A/B  frequencies  and  scan  angles,  and  with  a  more  rigorous  treatment  of  absorption  and 
scattering. 
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4.1.  Model 

Previous  studies  of  microwave  precipitation  effects  (Mugnai  et  al.  1990;  Smith  et  al. 

1991;  Mugnai  et  al.  1993;  Smith  et  al.  1994;  Kummerow  1996;  Panegrossi  et  al.  1998;  Bauer  and 
Mugnai  2003;  Panegrossi  2004)  have  coupled  fine-scale  numerical  weather  prediction  (NWP) 
model  output  with  a  radiative  transfer  model  (RTM)  using  1-D  plane  parallel  geometry,  and 
(typically)  the  Eddington  two-stream  approximation,  which  considers  a  single  upward  and 
downward  intensity  stream,  each  at  a  fixed  incidence  angle,  and  uses  a  scattering  phase  function 
simplified  to  for  the  limited  number  of  possible  scattering  angles. 

In  the  mid-1990's,  Petty  and  other  authors  introduced  Reverse  Monte  Carlo  (RMC) 
techniques  to  microwave  radiative  transfer  modeling.  The  Appendix  discusses  details  of 
previous  RMC  models,  as  well  as  the  one  developed  for  this  study.  Briefly,  RMC  models  trace 
photons  backward  from  a  sensor  through  a  medium,  following  them  as  they  are  scattered 
(perhaps  multiple  times),  and  ultimately  absorbed.  RMC  models  make  use  of  the  reciprocity 
theorem,  which  postulates  that  scattering  affects  photons  identically  regardless  of  which 
direction  they  are  traveling.  The  same  phase  function  describes  the  probability  of  a  forward¬ 
traveling  photon  scattering  into  a  given  direction  or  the  probability  that  a  scattered  backward¬ 
traveling  photon  came  from  a  given  direction.  Similarly,  Kirchhoff  s  Law  allows  us  to  consider 
the  point  at  which  a  backward-traveling  photon  is  absorbed  to  be  the  point  at  which  the  photon 
would  have  been  emitted  had  it  been  traveling  in  the  forward  direction. 

The  nature  of  RMC  models  enables  any  arbitrary  3-D  geometry  of  absorbers  and 
scatterers  to  be  constructed,  which  eliminates  the  need  for  1-D  plane-parallel  geometry  and 
limited  beam  angles.  The  flexibility  and  fidelity  of  RMC  techniques  to  the  physical  processes 
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being  simulated  comes  at  the  price  of  very  large  computational  expense.  Each  photon  must  be 
tracked  backward  through  possibly  many  scattering  events.  Each  of  those  events  requires 
generating  a  random  path  length  to  the  next  extinction  event,  randomly  determining  whether  the 
extinction  event  results  in  scattering  or  absorption,  randomly  determining  the  scattering 
direction,  then  performing  the  associated  propagation  direction  transforms,  position  updates,  and 
absorption  point  book-keeping.  A  large  number  of  photons  is  generally  required  for  the  RMC 
results  to  converge.  Even  then,  some  amount  of  random  noise  will  still  exist  in  the  results.  But 
as  computer  speed  and  memory  capacity  continues  to  rapidly  improve,  RMC  techniques  are 
becoming  more  practical.  The  reason  for  using  in  this  application  is  to  replicate  scattering 
processes  as  faithfully  as  possible  and  to  allow  generalization  to  3-D  geometry. 

As  discussed  in  Chapter  3,  the  Mie  extinction  and  scattering  properties  for  a  collection  of 
particles  at  a  point  in  space  can  be  derived  from  the  precipitation  densities  at  that  point.  We  need 
to  convert  these  densities  to  rain  rates  in  order  to  use  the  Marshall-Palmer  or  Sikhon-Srivastava 
empirical  size  distributions  to  determine  the  number  density  at  each  particle  size.  From  the 
particle  size  spectrum,  we  can  use  Mie  theory  to  compute  extinction  efficiency,  single-scatter 
albedo,  and  asymmetry  parameter  or  the  Mueller  matrix  elements.  The  Appendix  discusses  the 
five-step  process  of  translating  a  meteorological  profile  to  Mie  properties  and  then  conducting 
the  RMC  radiative  transfer. 

4.2.  Data  and  Methodology 

The  goal  of  this  simulation  study  is  to  simulate  the  AMSU-A  and  -B  TB's  produced  by  a 
number  of  prescribed  dry,  humid,  warm  rain,  tropical  stratiform,  and  tropical  cumulonimbus 
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precipitation  profiles.  Using  the  modeled  TB's,  we  will  compare  sounding  channel  precipitation 
attenuation  with  window  channel  precipitation  attenuation.  The  attenuation  is  quantified  by 
comparing  the  simulated  TB's  at  each  channel  with  TB's  modeled  with  rain,  snow,  and  graupel 
densities  set  to  0.  Petty  (2001a)  constructed  a  1-D  cloud  model  for  use  with  a  RMC  model  to 
study  precipitation  effects  at  Special  Sensor  Microwave  Imager  (SSMI)  frequencies  and  fixed 
view  geometry.  The  clear/dry,  humid,  warm  rain,  and  tropical  stratiform  profiles  used  in  this 
study  are  created  using  Petty's  parametric  rain  cloud  model,  which  uses  19  user-defined 
parameters,  along  with  detailed  aggregation  and  conversion  processes  to  generate  realistic  rain, 
snow,  and  graupel  profiles. 

Petty's  model  considers  only  droplets  falling  through  still  air,  and  is  not  intended  for 
convective  situations  with  strong  updrafts,  where,  for  instance,  supercooled  liquid  or  large  rimed 
graupel  particles  can  be  lofted  to  the  upper  troposphere.  To  obtain  realistic  tropical  convective 
profiles,  output  from  a  Mesoscale  Model  5  (MM5)  simulation  of  Hurricane  Bonnie  (1998)  was 
used.  This  particular  simulation  was  originally  run  at  the  NOAA  Environmental  Technology 
Laboratory  with  the  purpose  of  producing  simulated  radiances  for  an  experimental  geostationary 
infrared  sounding  instrument.  The  simulation  ran  for  96  hours  at  3-hour  time  steps  from  0000 
UTC  23  Aug  98  to  0000  UTC  27  Aug  98.  The  model  domain  is  a  297  x  249  horizontal  grid  with 
6.7  km  spacing  and  59  sigma  levels  between  reference  pressure  levels  of  1000  hPa  and  10  hPa. 
The  grid  covered  roughly  19°  to  37°  N  latitude  and  63°  to  82°  W  longitude  off  the  southeastern 
United  States.  Version  3.4  of  the  MM5  was  used,  with  initialization  data  from  the  National 
Center  for  Environmental  Prediction  /  National  Center  Atmospheric  Research  (NCEP/NCAR) 
Reanalysis  Project  (NNRP).  Physics  package  options  selected  for  the  run  were:  the  Kain-Fritsch 
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cumulus  parameterization,  the  MRF  (Hong-Pan)  boundary  layer  scheme,  the  Reisner  II 
microphysics  scheme  (which  features  supercooled  water,  slow-melting  snow,  and  graupel 
riming),  and  the  cloud-radiation  scheme.  Precipitation  output  at  each  grid  point  consists  of  bulk 
densities  for  rain,  snow,  and  graupel  categories,  as  well  as  cloud  liquid,  cloud  ice,  and  specific 
humidity. 

19  different  profiles  were  modeled.  For  each  profile,  the  RMC  model  simulated  AMSU- 
A  TB1  through  TBS  and  TB15  for  instrument  fields  of  view  (FOV's)  1,  5,  8,  11,  and  15  and 
AMSU-B  IB  16  through  TB20  for  FOV's  1, 1 1,  22,  33,  45.  The  channel  selection  covers  all 
relevant  temperature  sounding,  moisure  sounding,  and  window  channels;  AMSU-A  stratospheric 
temperature  sounding  channels  ( TB9  through  TB14)  were  omitted.  The  FOV  selection  is  spread 
evenly  across  scan  angles  from  edge  of  scan  (approximately  48  1/3°  scan  angle,  producing  an 


Table  4. 1 :  Simulation  Cases 
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Case  4:  60%  Uniform  Relative  Humidity;  No  Precipitation 
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Figure  4.1:  (a)  Profiles  of  temperature,  water  vapor,  cloud  water,  rain,  snow,  and  graupel  used 
for  simulation  case  4.  (b)  Resulting  window  channel  weighting  functions  produced  by  the 
reverse  Monte  Carlo  radiative  transfer  model  for  AMSU-A  channels  2  (31  GHz)  and  15  (89 
GHz),  FOV's  1,  8,  and  15  and  for  AMSU-B  channels  17  (150  GHz)  and  18-20  (183  GHz), 
FOV's  1,  23,  and  45.  The  top  and  bottom  level  weights  indicate  the  fractional  TB  contribution 
of  the  cosmic  background  and  the  ocean  surface,  respectively. 
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Figure  4.2:  Same  as  Figure  4.1,  but  for  case  8. 
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Case  9:  Tropical  Stratiform  Rain 
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Figure  4.3:  Same  as  Figure  4.1,  but  for  case  9. 
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Earth-incidence  angle  of  approximately  60°)  to  nadir.  5000  photons  were  used  for  each  channel 
and  FOV  combination.  The  RMC  code  was  run  twice  for  each  profile;  once  with  precipitation 
included  and  once  with  rain,  snow,  and  graupel  removed  (water  vapor  and  cloud  were  left  in 
place).  The  precipitation  attenuation,  ATB,  for  each  channel  is  computed  by  subtracting  the  non¬ 
precipitating  TB  from  the  precipitating  TB. 

Table  4.1  summarizes  the  meteorological  conditions  for  each  of  the  nineteen  profiles. 
Cases  1  through  6  were  non-cloudy,  non-precipitating  scenes  with  temperature  profiles  from 
Petty  (2001a)  and  uniform  relative  humidity  throughout  the  depth  of  the  atmosphere  ranging 
from  0%  (case  1)  to  100%  (case  6).  The  top  half  of  Figure  4.1  depicts  temperature  and  water 
vapor,  cloud  liquid,  rain,  snow,  and  graupel  density  profiles  for  Case  4  (60%  RH).  Profile  7  is 
the  oceanic  warm  rain  profile  from  Petty  (2001a),  which  consists  of  a  somewhat  dry  atmosphere 
with  cloud  and  rain  confined  below  the  freezing  level.  Case  8  (Figure  4.2)  is  modified  to 
increase  the  rain  rate.  Profile  9  is  Petty's  tropical  stratiform  rain  profile,  which  has  a  thick  snow 
and  graupel  layer  extending  from  the  freezing  level  to  10  km  (Figure  4.3). 

Figure  4.4  depicts  cloud  liquid  water,  column  integrated  rain  water,  and  column 
integrated  graupel  for  the  MM5  Hurricane  Bonnie  simulation  at  1200UTC  25  Aug  98  and  1200 
UTC  26  Aug  98  (60  and  84  hours  after  forecast  initialization).  Profiles  10-14  are  from  1200 
UTC  25  Aug.  At  this  time,  the  simulated  hurricane  had  an  intensity  of  958  hPa.  Five  profiles 
are  arrayed  in  an  "x"  around  the  storm  center:  Case  10  is  the  southwest  point,  located  in  a 
relatively  clear  area;  case  1 1  is  located  in  intense  spiral  band  precipitation  to  the  northwest;  case 
12  is  in  the  eye;  cases  13  and  14  are  in  spiral  band  precipitation  to  the  southeast  and  northeast 
respectively.  Figure  4.5  depicts  the  profile  from  the  northwest  spiral  band  and  contains  the 
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Figure  4.4:  Column  integrated  cloud  liquid  water  (top),  rain  (middle),  and  graupel  (bottom) 
for  MM5  simulation  of  Hurricane  Bonnie  (1998)  at  1200  UTC  25  Aug  1998  (right  column) 
and  1200  UTC  26  Aug  1998  (left  column).  Bold  +'s  indicate  grid  points  selected  for 
simulation  cases  10  through  14  (25/12  UTC)  and  15  through  19  (26/12  UTC). 
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Case  11:  Bonnie  1200UTC  25  Aug  98  Northwest  Spiral  Band 
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Figure  4.5:  Same  as  Figure  4.1,  but  for  case  1 1 
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Case  15:  Bonnie  1200UTC  26  Aug  98  West  Spiral  Band 
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Figure  4.6:  Same  as  Figure  4.1,  but  for  case  15. 
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heaviest  precipitation  simulated  in  any  of  the  cases.  Profiles  15-19  are  from  Bonnie  at  1200 
UTC  on  the  26th,  when  the  simulated  Bonnie  had  improved  its  organization,  intensified  to  948 
hPa,  and  was  nearing  landfall.  These  profiles  are  arranged  in  a  "+"  around  the  storm.  Figure  4.6 
depicts  the  profile  from  Bonnie's  western  spiral  band,  and  is  representative  of  spiral  band 
precipitation. 

4.3.  Results 

Accompanying  each  of  the  profiles  in  the  top  halves  of  Figures  4.1,  4.2,  4.3,  4.5,  and  4.6 
are  the  AMSU-A  and  -B  window  and  moisture  channel  weighting  functions  produced  by  the 
RMC  model  from  the  profile  in  the  top  half  of  the  figure.  The  weighting  functions  depict  the 
fraction  of  photons  in  the  RMC  model  that  originate  at  each  level.  The  weighting  functions  also 
include  the  fraction  of  photons  emitted  by  the  ocean  surface  or  the  cosmic  background,  and 
include  those  weights  at  the  bottom  and  top  levels  of  the  function  respectively.  For  window 
channels  in  less  than  opaque  conditions,  the  sea  surface  and  cold  space  are  prime  contributors  to 
the  observed  TB,  and  the  proportions  of  each  correspond  to  the  sea  surface  emissivity  and 
reflectivity,  respectively,  matching  the  frequency  and  scan  angle. 

Figure  4.7  depicts  TB15  and  TB17,  the  high-freq  window  channels,  vs.  TB1  and  TB2,  the 
low-frequency  windows,  for  each  case  and  scan  angle.  Figure  4.8  does  the  same  for  TB18 
through  TB20,  the  moisture  sounding  channels.  In  each  plot,  the  modeled  TB  with  precipitation 
turned  off  in  the  RMC  model  is  depicted  as  a  dot,  while  the  corresponding  TB  with  precip 
included  is  plotted  with  a  circle.  The  dashed  arrow  tracks  TB  evolution  as  RH  increases  from  0% 
to  100%  in  the  cloud-  and  precipitation-free  cases  (1  through  6).  The  primary  effect  in  these 
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Figure  4.7:  Simulated  high  frequency  window  channel  TB's  vs.  Low  Frequency  window  channel  TB's.  AMSU-A  TB15  (left 
column)  and  AMSU-B  TB17  (right  column)  vs.  AMSU-A  TB1  (bottom  row)  and  AMSU-A  TB2  (top  row).  Dots  indicate 
TB's  modeled  with  precipitation  excluded.  Circles  indicate  precipitating  TB's.  Dashed  arrow  follows  modeled  TB's  as 
relative  humidity  increases  from  0%  to  100%.  Dashed  circles  indicated  modeled  TB's  from  given  types  of  precipitation. 
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cases  is  on  TB18  through  TB20,  for  which  the  atmosphere  rapidly  becomes  opaque  (between  0% 
and  20%  RH,  as  indicated  by  the  large  TB  jump  between  cases  1  and  2)  and  their  weighting 
functions  stratify  with  18  highest  (~  9km),  19  in  the  middle  (~  7  km),  and  20  lowest  (~  4  km). 

As  RH  increases,  all  three  of  these  channels'  TB  decrease  due  to  their  weighting  functions 
ascending  in  the  increasingly  opaque  atmosphere. 

TB1 7  behaves  similarly  to  TB20,  since  it  is  strongly  affected  by  water  vapor  continuum 
absorption;  it  rapidly  increases  with  increasing  RH  to  a  maximum  at  around  40%  RH,  then 
slowly  decreases  as  its  weighting  function  also  ascends.  The  atmosphere  is  remains  slightly 
transparent  to  the  other  high-frequency  window  channel,  TB15,  which  is  still  increasing  as  RH 
reaches  100%. 

There  is  an  important  difference  in  the  low-frequency  window  TB  response  to  RH 
increase.  TB1  (23.8  GHz)  is  near  the  weak  22  GHz  water  vapor  line,  so  it  is  more  susceptible  to 
water  vapor  emission  and  it  increases  much  more  with  increasing  RH  than  TB2  (31.4  GHz).  The 
high-frequency  and  moisture  sounding  channels  reach  their  peak  values  over  a  smaller  range  of 
TB2  variation,  then  saturate  and  remaining  relatively  flat  as  TB2  increases.  This  should  make 
TB2  better  for  predicting  no-scatter  window  and  moisture  channel  TB's  because  uncertainty  in 
TB2  will  have  little/no  affect  on  relatively  constant  high-frequency  TB  once  the  atmosphere  has 
become  opaque.  By  contrast,  TB1  and  TB15  or  TB17  increase  together  nearly  linearly,  so 
uncertainty  in  TB1  anywhere  along  its  range  of  values  could  introduce  large  uncertainty  in 
predicted  no-scatter  TB15  or  TB17  using  TB1. 

The  warm  rain  cases  (7  and  8),  by  virtue  of  their  relatively  low  RH  (TPW  close  to  the 
uniform  40%  RH  case)  and  the  thinness  of  the  rain  layer,  act  very  similarly  to  the  cases  with 
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increasing  relative  humidity.  They  cause  no  appreciable  scattering  effect  on  any  channel. 

By  contrast,  the  Petty  2001  stratiform  rain  case  (9)  introduces  significant  scattering.  It 
consists  of  a  "snow-generating  layer"  extending  from  10  km  down  to  the  cloud  top  at  7  km. 
Accretion  and  conversion  processes  within  the  cloud  then  generate  a  mixed  graupel  and  snow 
layer  extending  down  to  the  freezing  level  at  4.4  km.  Snow  melts  instantaneously,  while  graupel 
takes  0.5  km  to  melt,  creating  a  bright-band  (large  graupel  particles,  with  steadily  increasing 
water  content,  which  giving  them  very  large  extinction  coefficient  and  moderate  single-scatter 
albedo)  near  the  freezing  level.  Rain  extends  to  the  surface.  The  combined  net  effect  of  these 
precipitation  contributors  is  an  approximately  40  K  warming  for  the  low-frequency  window  TB's 
due  to  rain  emission.  TB18  is  unaffected  since  the  bulk  of  its  weighting  function  is  above  most 
of  the  snow  and  graupel.  TB19  is  slightly  affected,  while  TB20  and  TB17  (with  their  similar 
weighting  function  peaks)  are  hit  hardest;  both  suffer  a  70-80  K  depression.  TB15  is  changed 
little  because  the  increased  scattering  TB  depression  caused  by  the  ice-phase  precipitation  layer 
is  offset  by  emission  from  the  thick  rain  layer,  which  replaces  the  cold  sea  surface.  This  is  the 
first  case  to  exhibit  sounding  channel  TB  depression.  As  discussed  in  Chapter  3,  this  is  due  to 
the  reflecting  ice  layer  impinging  on  the  low-level  tail  of  the  sounding  channel  weighting 
function,  replacing  warm  lower-troposphere  emission  with  reflected  cold  upper-troposphere 
emission. 

Each  of  the  two  Hurricane  Bonnie  output  scenes  introduced  a  variety  of  profile  types. 
Both  had  a  relatively  clear,  but  humid,  profile  within  the  eye  and  in  one  or  two  of  the 
surrounding  points.  Both  scenes  also  produced  profiles  with  weak  to  moderate  spiral  band 
precipitation.  This  type  of  profile  produces  a  10-20  K  increases  in  TB2  due  to  rain  emission,  a 
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10-30  K  decreases  in  77?  75  and  TB17  due  to  scattering,  little  effect  on  TB18  due  to  the  height  of 
its  weighting  function,  modest  effects  on  the  remaining  moisture  channels  (~10  K  depression  in 
TB19  and  ~20  K  depression  in  TB20),  and  weak  low-  and  mid-tropospheric  temperature 
sounding  channel  TB  depression,  but  little  or  no  TB8  depression. 

The  northwestern  profile  from  the  1200  UTC  25  Aug  98  scene  (case  11)  produced  a 
profile  with  extremely  heavy  convective  precipitation  containing  snow  as  high  as  18  km  and 
graupel  as  high  as  15  km.  Thick  cloud  extends  from  15  km  down  to  the  melting  level  at  around 
7  km,  and  heavy  rain  lies  below.  In  this  case,  TB2  warms  to  about  240  K,  the  warmest  value 


seen  in  any  of  the  simulated  profiles,  because  bright-band  emission  and  opaque  rain  completely 
obscure  the  sea  surface.  All  of  the  high-frequency  channels  have  sharp  peaks  in  their  weighting 
functions  at  around  12-13  km  in  the  ice  layer,  and  all  of  their  TB's  converge  to  about  230  K. 


That  represents  TB  depressions  ranging  from 
10  K  (TB18)  to  60  K  (77?/ 7).  The  ice  layer 
extends  through  all  tropospheric  temperature 
sounding  channels'  weighting  functions, 
including  channel  8's,  so  not  surprisingly,  this 
case  produced  the  most  significant  sounding 
channel  ATB  of  the  ensemble  (0.5  K  for  TB8,  3 
K  for  TB7,  8  K  for  TB6,  and  18  K  for  TB5). 

Figure  4.9  depicts  correlation 
coefficients  relating  sounding  channel  ATB  to 
high-frequency  window  and  moisture  channel 


Figure  4.9:  Correlations  between  high- 
frequency  window/moisture  channel  ATB 
and  temperature  sounding  channel  ATB. 
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Figure  4.10:  Simulated  sounding  channel  scattering 
depression  (dependent  variable;  ATB8  through  ATB5  in 
rows  descending  from  top)  vs.  ATB15  (left  column) 
and  ATB17  (right  column).  Correlation  coefficients  for 
each  combination  are  shown  in  Figure  4.9. 
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MC  Sounding  vs.  Moisture  Ch  Depression 
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Figure  4.1 1:  Same  as  Figure  4.10,  but  with  AMSU-B  moisture  sounding  channel  scattering 
depression  (ATB18  through  ATB20 )  as  the  independent  variable.  Correlations  with  sounding 
channel  ATB  are  in  Figure  4.9. 
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ATB.  Correlations  are  uniformly  weak  for  ATB8,  and  increase  with  lower-level  channels. 

ATB18  and  ATB  19  perform  best  across  the  board,  with  correlations  approaching  0.9  for  ATB5. 
Figures  4.10  and  4.1 1  show  scattering  depression  for  each  of  the  four  sounding  channels  vs.  both 
high-frequency  window  channels’  (4.10)  and  all  three  moisture  sounding  channels'  (4.1 1)  ATB. 
Least-squares  best-fit  lines  are  superimposed. 

Bauer  (2001)  and  Petty  (2001)  demonstrated  that  principal  component  analysis  was 
useful  for  isolating  multi-channel  precipitation  effects.  Figure  4.12  depicts  the  first  empirical 
orthogonal  function  ( EOF1 )  of  the  standardized  ATB  for  the  combined  set  of  low-frequency, 
high-frequency,  and  moisture  channels  ( ATB1 ,  ATB2,  and  ATB  15  through  ATB20).  This  EOF 
explains  72%  of  the  variation  in  the  data  set.  The  sign  of  an  EOF's  elements  is  arbitrary,  since 
the  sign  of  the  corresponding  principal  components  is  chosen  to  produce  the  correct  physical 
result.  What  is  important  is  the  relative  magnitude  of  each  element,  and  the  sign  of  that 
individual  element  with  respect  to  the  rest.  For  this  EOF,  ATB1  and  ATB2  are  the  opposite  sign 
of  ATB15  through  ATB20.  This  is  consistent  with  rain  emission  warming  the  low-frequency 
window  TB's  and  scattering  reducing  the  high-frequency  TB's.  Figure  4.13  depicts  ATB  for  each 
sounding  channel  vs.  the  first  principal  component  of  the  combined  window  channel  ATB  set 
with  accompanying  correlations.  Results  are  similar  to  using  an  individual  high-frequency 
channel  ATB.  The  correlation  with  PCI  decreases  from  about  0.9  for  ATB5  to  about  0.3  for 
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Figure  4.12:  Leading  EOF  (left)  of  simulated  window  and  moisture  channel 
ATB.  Opposing  sign  of  TB1  and  2  vs.  TB15-20  is  consistent  with  rain 
emission  warming  TB1  and  2  while  scattering  reduces  TB 15-20.  Variance 
explained  by  the  6  leading  EOF's  of  window  and  moisture  channel  ATB. 
EOF1  explains  72%  of  the  multi-channel  variance. 
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Figure  4.13:  Simulated  sounding  channel  ATB  vs.  the  leading  principal  component  of 
window  or  moisture  channel  ATB  (corresponding  to  EOF1  in  Figure  4.12).  Correlations 
with  each  sounding  channel  ATB  (at  right)  are  comparable  to  correlations  achieved  using 


ATB  18  or  ATB  19  alone. 
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4.4.  Discussion 

The  flatter  curves  for  high-frequency  and  moisture  channel  TB  vs.  TB2  indicate  that  TB2 
is  probably  a  better  choice  than  TB1  to  use  for  predicting  the  no-scatter  high-frequency  TB's. 

TB1  responds  strongly  to  both  humidity  variation  and  rain  emission,  while  TB2  is  much  less 
sensitive  to  humidity.  Since  we  are  interested  only  in  precipitation  effects,  excluding  TB1 
humidity  variation  as  a  source  of  uncertainty  in  predicted  no-scatter  high-frequency  TB's  should 
increase  their  accuracy. 

The  best  correlated  high-frequency  ATB  for  each  sounding  channel  ATB  is  either  TB18  or 
TB19.  However,  it  is  significant  to  point  out  that  the  no-scatter  TB  values  are  known  exactly  in 
these  simulation  cases.  TB18  depressions  are  very  small  (10  K  at  most),  so  it  is  doubtful-given 
the  uncertainty  in  no-scatter  high-frequency  TB's  predicted  using  7B2-that  in  practice  the  no¬ 
scatter  TB18  can  be  estimated  with  that  degree  of  accuracy. 

Principal  component  analysis  clearly  highlights  the  multi-channel  scattering  signature, 
but  because  of  the  redundancy  of  information  provided  by  the  scattering-sensitive  high- 
frequency  channels,  the  first  PC  doesn't  perform  significantly  better  than  a  single  high-frequency 
channel  for  predicting  sounding  channel  ATB. 

These  results  suggest  that  correcting  TB8  for  precipitation  scattering  effects  is  going  to  be 
problematic.  Simulated  scattering  ATB8's  were  on  the  same  order  as  the  Monte  Carlo  statistical 
noise,  which  suggests  that  scattering  effects  on  channel  8  may  not  be  discemable  from  real-world 
sources  of  noise.  The  TB8  weighting  function  peaks  near  15  km  at  nadir  (higher  at  edge  of  scan), 
so  even  in  vigorous  eyewall  convection,  it  is  difficult  for  a  significant  number  of  frozen 
hydrometeors  to  impact  this  channel  significantly  (just  as  it  is  difficult  to  affect  TB18,  which  also 
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peaks  very  high  in  the  troposphere  in  humid  conditions).  A  promising  approach  may  be  to  use 
moisture  channels  that  peak  near  the  same  level  to  correct  individual  sounding  channels,  i.e. 
TB18  to  correct  TB7  and  TB8,  TB19  to  correct  TB6,  and  TB20  to  correct  TBS. 

The  simulations  in  this  chapter  support  the  hypothesis  from  Chapter  3— namely,  that 
precipitation  scattering  depression  of  sounding  channel  TB's  can  be  quantified  by  using  the 
depression  in  high-frequency  window  or  moisture  sounding  channels.  It  appears  that  TB2,  which 
is  only  weakly  affected  by  scattering,  will  allow  prediction  of  high-frequency  no-scatter  TB's, 
from  which  scattering  depression  can  be  determined.  Then  sounding  channel  ATB  can  be 
inferred,  and  corrected  for,  by  using  the  empirical  relation  between  high-frequency  and  sounding 
channel  ATB.  The  next  chapter  presents  analysis  of  observed  AMSU-A  and  -B  data  and 
develops  such  a  correction  technique,  which  is  then  used  in  Chapter  6  to  improve  TC  intensity 


estimates. 
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5.  Developing  an  Empirical  Scattering  Correction 

The  conceptual  model  presented  in  Chapter  3  suggested  using  a  low-frequency  window 

channel  to  predict  high-frequency  window  or  moisture  channel  non-scattering  TB,  in  a  manner 

/ 

similar  to  the  Grody  (1991, 1999)  scattering  index.  The  simulation  study  in  Chapter  4 
demonstrated  that  scattering-induced  high-frequency  window  or  moisture  channel  TB  depression 
(ATB)  from  predicted  no-scatter  values  is  proportional  to  sounding  channel  TB  depression  due  to 
precipitation.  The  simulation  results  also  suggest  that  TB2  would  be  a  better  low-frequency 
predictor  of  high-frequency  no-scatter  TB  and  that  ATB  18  or  ATB  19  might  be  the  best  predictors 
for  sounding  ATB. 

This  chapter  will  apply  those  concepts  to  a  large  set  of  observed  AMSU  data  to  verify 
them  in  practice  and  develop  an  empirical  correction  for  sounding  channel  precipitation  effects. 
Chapter  6  will  apply  the  correction  technique  and  validate  the  performance  of  scattering- 
corrected  TB' s  in  the  CIMSS  AMSU-based  TC  intensity  estimation  technique. 

5.1.  Data  and  Methodology 

Since  a  TC's  upper-tropospheric  warm  anomaly  may  be  on  the  order  of  10  K,  and  the 
simulated  scattering  effect  on  AMSU-A  sounding  channels  was  typically  on  the  order  of  0.1  to 
10  K,  it  will  be  difficult  to  isolate  the  TB  reduction  caused  by  scattering  if  observed  only  near 
TC  centers.  So  this  chapter  begins  with  an  analysis  of  AMSU-A  and  -B  data  from  tropical 
scenes  without  a  TC  present. 

For  this  analysis,  four  days  of  archived  AMSU  data  are  used:  01  Jan,  01  Apr,  01  Jul,  and 
01  Oct  2003.  Three  AMSU  instruments  were  operational  at  the  time,  flying  on  the  NOAA-15,  - 
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shape  for  each  AMSU-B  FOV  in  the  3x3  array  with  the  AMSU-A  3  l/3°-width  Gaussian  beam 
centered  at  the  middle  of  the  array  produces  each  AMSU-B  FOV's  relative  contribution  to  a 
synthetic  AMSU-A  resolution  FOV.  The  comer  points  of  the  3  x  3  AMSU-B  array  contribute 
0.096  each  to  the  convolved  TB,  the  edge  points  contribute  0.1 18,  and  the  center  FOV 
contributes  0.144.  For  all  further  analysis,  AMSU-B  data  is  first  convolved  to  AMSU-A 
resolution. 

5.1.2.  Limb  Correction 

Limb  darkening,  the  decrease  in  TB  with  increasing  scan  angle  for  tropospheric  sounding 
channels,  is  caused  by  the  longer  slant  path  traveled  through  the  absorbing  atmosphere.  This 
effect  can  reduce  observed  TB  by  as  much  as  10  K  from  nadir  to  edge  of  scan  as  the  weighting 
function  ascends  to  approximately  the  level  of  the  next  higher  channel  at  nadir  (for  example,  the 
channel  7  weighting  function-which  peaks  at  approximately  250  hPa  at  nadir-peaks  at  about 
150  hPa  at  edge  of  scan,  the  same  level  at  channel  8  at  nadir;  see  Figure  2.15).  In  order  to 
maintain  uniform  sensitivity  to  a  fixed  atmospheric  level  across  the  full  AMSU-A  scan,  a  limb 
correction  technique  must  be  employed,  and  the  CIMSS  AMSU-based  TC  intensity  algorithm 
follows  the  one  developed  by  Goldberg  et  al.  (2001).  For  each  channel  and  FOV,  a  synthetic 
limb-corrected  TB  is  produced  using  a  weighted  sum  of  the  channel  and  its  neighbors  above  and 
below.  Goldberg  et  al.  (2001)  derived  the  coefficients  of  the  weighted  sum  using  a  large  sample 
of  AMSU  data  stratified  by  latitude,  ocean  or  land  surface,  and  time  of  year.  The  standard 
deviation  of  the  resulting  limb-corrected  TB' s  is  within  the  instrument  noise.  The  limb- 
correction  coefficients  can  be  thought  of  as  constructing  a  synthetic  weighting  function  that 
maintains  its  approximate  shape  and  peak  height  across  the  scan  line.  By  using  limb-corrected 
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sounding  channel  TB's,  the  CIMSS  TC  algorithm  ensures  that  a  channel  remains  sensitive  to 
temperatures  at  a  nearly  constant  level,  regardless  of  scan  angle.  For  this  analysis,  using  limb- 
corrected  sounding  TB's  ensures  that  their  sensitivity  to  hydrometeors  will  remain  relatively 
constant  across  the  scan  line  as  well. 

Goldberg  et  al.  (2001)  were  also  able  to  produce  limb-correction  coefficients  for  the 
window  channels.  This  is  less  straightforward  than  for  sounding  channels  because  the  window 
channels  usually  have  a  strong  surface  contribution  and  each  has  a  different  response  to 
humidity.  The  authors  chose  channels  1,  2,  and  4  as  contributors  to  the  limb-correction  for  TB1 
and  TB2  and  channels  2,  3,  and  15  as  contributors  to  limb-corrected  TB15.  This  study,  however, 
intends  to  use  each  channel's  unique  response  to  moisture,  cloud,  rain,  and  ice.  Blending 
multiple  window  channels  by  way  of  limb  correction  will  mix  these  effects  to  different  degrees 
at  different  scan  angles,  complicating  the  task  of  quantifying  scattering. 

A  final  note  about  limb  correction:  the  coefficients  derived  by  Goldberg  et  al.  (2001), 
and  used  in  operations  by  the  CIMSS  TC  intensity  technique,  produced  a  large  amount  of  noise 
(1  to  2  K)  from  FOV  to  FOV  across  a  scan  line.  The  underlying  raw  brightness  temperatures  all 
vary  smoothly  across  the  scan  line,  indicating  that  the  noisiness  is  produced  by  the  limb 
correction  coefficients.  The  noise  in  channel  8  is  on  the  same  order  as  the  expected  TB 
precipitation  perturbation  predicted  by  the  simulation  results  in  Chapter  4.  For  this  reason,  new 
limb-correction  coefficients  were  derived  for  the  sounding  channels  using  this  data  set.  The 
resulting  coefficients  for  channel  8  are  considerably  smoother  and  reduce  the  FOV  to  FOV  noise. 
Those  limb  coefficients  are  used  to  limb-correct  AMSU-A  sounding  channels  for  all  analyses  in 
this  and  the  subsequent  chapter. 
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16,  and  -17  spacecraft  (NOAA-17's  AMSU-A  has  since  failed).  Each  day's  data  includes 
AMSU-A  and  -B  scan  lines  which  are  completely  over  ocean  and  confined  between  15°  S  and 
15°  N  latitude.  The  over-ocean  requirement  eliminates  any  effects  of  land-surface  emissivity  on 
the  results.  The  latitude  bounds  were  chosen  to  minimize  the  effect  of  tropopause  height 
variation  on  the  results,  particularly  for  AMSU-A  channel  8,  which  is  sensitive  to  temperature  in 
the  upper  troposphere  near  150  hPa.  The  full  sample  contained  375,325  individual  fields  of  view 
(FOV).  After  the  initial  analysis  using  the  TC-free  tropical  data,  a  second  analysis  was 
conducted  for  comparison,  using  30  FOV  x  30  FOV  boxes  centered  on  TC's  from  the  2003 
Atlantic  hurricane  season.  This  additional  data  set  contained  1 12,500  FOV's  from  125  different 
TC  scenes. 

5.1.1.  AMSU-B  to  AMSU-A  Convolution 

Since  AMSU-B  has  3  times  finer  resolution  than  AMSU-A,  some  method  needs  to  be 
used  to  convolve  AMSU-B  brightness  temperatures  from  their  native  resolution  to  the  same 
resolution  as  the  AMSU-A  FOV's  with  which  they  will  be  used.  Bennartz  (2000)  demonstrated  a 
convolution  method  using  a  Backus-Gilbert  technique  on  observed  AMSU-B  data.  The  resulting 
convolution  coefficients  account  for  both  the  instrument  resolution  differences  as  well  as 
asymmetries  in  instrument  response  across  the  scan  line. 

The  convolution  approach  used  for  this  analysis  is  simpler  and  does  not  account  for 
variation  in  antenna  response  and  co-location  across  the  instruments'  scan  lines.  AMSU  antenna 
patterns  consist  of  a  strong  main  lobe  centered  on  the  instrument  line  of  sight,  with  half-power 
beam-widths  of  3  1/3°  for  AMSU-A  and  1.1°  for  AMSU-B.  Each  AMSU-A  FOV  corresponds  to 
a  3  x  3  array  of  AMSU-B  FOV’s.  Convolving  the  2-dimensional  l.l°-width  Gaussian  beam 
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5.2.  Results 

Figures  5.1  through  5.4  are  images  of  AMSU  TB's  from  one  satellite  (NOAA-16)  on  one 
of  the  four  days  (01  Jan  03)  for  the  western  hemisphere.  Five  ascending  orbits  and  seven 
descending  orbits  throughout  the  day  are  merged  to  form  a  single  image,  which  accounts  for  the 
discontinuities  in  the  images.  A  scan  line  is  only  used  if  it  is  completely  over  ocean  (a 
requirement  driven  by  the  limb-correction  development  process),  so  large  swaths  of  data  are 
missing  around  South  America  and  Africa. 

Figure  5.1  depicts  AMSU-A  window  channels:  1,  2,  3,  and  15.  A  TB  gradient  from  nadir 
to  edge  of  scan  is  evident  for  each  channel.  It  is  caused  by  limb  effect  in  the  humid  tropical 
atmosphere.  TB1  and  TB15  are  very  sensitive  to  humidity,  so  strong  meridional  gradients  are 
observed  in  these  channels  approaching  the  convection  in  the  intertropical  convergence  zone 
(ITCZ).  TB2  is  relatively  unaffected  by  humidity  so  it  exhibits  a  very  sharp  TB  gradient  at  the 
edge  of  the  ITCZ  precipitation.  TB3  is  more  difficult  to  interpret  because  it  is  sensitive  to  lower- 
tropospheric  temperature,  in  addition  to  sea  surface  emission  and  humidity.  It  is  not  particularly 
useful  for  our  purpose  of  predicting  high-frequency  TB.  The  strong  mesoscale  convective 
system  (MCS)  affects  several  FOV's  near  7°  N  120°  W  and  offers  a  good  example  of  the  multi¬ 
channel  scattering  signature.  It  exhibits  a  strong  warm  signal  in  TB2,  which  contrasts  sharply 
with  the  strong  depression  in  TB15.  Discussion  of  the  following  images  will  focus  on  this  MCS. 

Figure  5.2  shows  limb-corrected  AMSU-A  temperature  sounding  channels.  TB4  has 
strong  surface  sensitivity,  so  rain  emission  blocks  the  radiometrically  cold  sea  surface  (as  well  as 
possibly  increasing  sea  sfc  emissivity  as  falling  rain  roughens  it)  The  MCS  causes  about  5  K 
depression  in  TB4.  The  TB5  weighting  function  peaks  near  700  hPa,  and  it  contains  a  small,  but 
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non-negligible,  surface  component.  Some  rain  warming  is  evident  in  this  channel,  as  well  as 
about  2  K  of  scattering  TB  depression  by  the  MCS.  TB6  is  sensitive  to  400  hPa  and  has  no 
surface  contribution.  Zonal  TB  variation  becomes  the  largest  source  of  variability  for  this  and 
higher  sounding  channels.  The  MCS  causes  about  1.5  K  scattering  depression.  TB7  (250  hPa) 
suffers  a  1  K  reduction  by  the  MCS.  Whatever  scattering  effect  is  present  in  TB8  is  lost  among 
about  the  approximately  0.5  K  of  noise  in  the  scene,  consistent  with  Chapter  4's  conclusion  that 
scattering  effects  on  TB8  will  be  hard  to  quantify. 

Figure  5.3  depicts  AMSU-B  window  channels  ( TB16 ,  89  GHz,  and  TB17, 150  GHz) 
convolved  to  AMSU-A  resolution.  While  they  use  the  same  frequency,  the  convolved  AMSU-B 
TB16  is  about  10  K  colder  than  AMSU-A  TB15  in  scattering  areas  due  to  better  ability  of 
AMSU-B  to  resolve  finer  details  that  get  smoothed  by  AMSU-A  resolution.  Even  after 
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Figure  5.3:  AMSU-B  window  channels  ( TB16  and  TB17  top  to  bottom)  convolved  to  AMSU- 
A  resolution. 
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convolution,  the  coldest  AMSU-B  FOV's  contribute  strongly  to  cooling  the  full  FOV.  The  TB17 
weighting  function  is  above  surface  in  humid  tropical  conditions,  so  the  scene  is  a  relatively 
uniform  280  K  background  with  cooler  TB's  in  precipitation  and  strong  TB  reduction  by  the 
MCS,  which  actually  consists  of  several  diffuse  cold  spots  in  TB17,  instead  of  single  sharp 
maxima  observed  in  other  window  and  sounding  TB's. 

AMSU-B  moisture  channels  ( TB18  through  TB20 ,  183.31  +/-  1,  3,  and  7  GHz)  are  shown 
in  Figure  5.4.  The  progressively  lower  weighting  functions  of  TB18,  TB19,  and  TB20  are 
evident  in  the  increasing  background  TB.  As  channels  penetrate  deeper,  they  experience  greater 
scattering  TB  depression.  TB18  is  reduced  the  least  (but  still  a  considerable  amount)  by  the 
MCS,  while  TB20  appears  similar  to  77377,  a  result  of  their  very  similar  weighting  functions  in 
humid  conditions.  An  important  feature  to  note  is  the  large-scale  variation  of  20-40  K  TB  in 
each  channel,  caused  by  large-scale  variation  in  precipitable  water  across  the  ocean  basins. 

Figure  5.5  contains  scatter  plots  of  each  high-frequency  window  or  moisture  channel  vs. 
TB2.  The  TB15  and  TB17  curves  qualitatively  match  the  simulation-predicted  results  (Figure 
4.7).  Each  follows  a  well-defined  logarithmic  curve  from  ( TB2 ,  TB15)  ~  (150  K,  200  K)  to 
(7772,  TB15)  ~  (200  K,  270  K)  as  relative  humidity  (or  slant  path)  increases.  TB15  then 
decreases  with  increasing  TB2.  The  simulation  results  showed  that  this  is  scattering-induced. 
Similarly,  TB17  undergoes  sharp  initial  rise  with  TB2  followed  by  quick  drop  as  scattering 
rapidly  becomes  important.  TB18  through  TB20  are  also  qualitatively  similar  to  the  simulation 
results.  Each  channel  begins  with  a  high  TB  (low  weighting  function  peak)  at  low  TB2.  TB18 
through  TB20  then  decrease  with  increasing  TB2  due  to  the  weighting  function  rising.  The  more 
sparsely  populated  triangular  region  of  depressed  TB  for  channels  18-20  appears  to  be  a 
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Figure  5.5:  AMSU-A  TB15  and  AMSU-B  TB17  through  TB20  (convolved  to  AMSU-A 
resolution)  vs.  AMSU-A  TB2  for  all  FOV's  used  in  the  tropical  ocean  scenes  shown  in  Figures 
5. 1-5.4. 

scattering  effect.  Note  that,  despite  little  scattering  effect  in  the  simulations,  TB18  does  exhibit 
some  depression. 

Figure  5.6  depicts  the  same,  except  for  75,319  FOV's  within  1000  km  of  a  TC  center. 
Land  was  not  filtered  out  of  these  scenes.  For  land  FOV's,  the  high  surface  emissivity  makes  all 
window  TB's  large.  Each  plot  has  a  cluster  of  clear  land  scenes  in  its  upper  right;  this  is  most 
distinct  in  TB15  and  TB17.  A  dense  line  of  points  connects  this  land  cluster  to  the  opaque 
humid/cloudy  atmosphere  point  ( TB2 ,  TB15)  ~  (200  K,  270  K).  Similar  dense  lines  of  points 
(most  visible  in  TB15,  TB17,  and  TB20 )  connect  the  land  cluster  and  the  opaque  humid/cloudy 
curve  to  a  point  near  ( TB2 ,  TB15)  « (260  K,  200  K)  or  ( TB2 ,  TB17)  =  (260  K,  140  K). 

There  are  similar,  but  less  distinct,  land  clusters  in  TB18  through  TB20  (all  of  whose 
surface  contributions  are  usually  weak),  along  with  triangular  regions  fimneling  points  toward  a 
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Figure  5.6:  Same  as  Figure  5.5,  but  for  2003  Atlantic  hurricane  scenes  (with  land  FOV's 
included).  Dashed  circles  highlight  multi-channel  TB's  from  clear  ocean  scenes, 
humid/cloudy  scenes,  clear  land  scenes,  and  FOV's  filled  with  convective  precipitation. 

Dashed  line  indicates  approximate  no-scattering  TB's. 

limiting  point  near  TB2  =  260  K.  The  maximum  depression  occurring  for  each  channel  near  TB2 

=  260  K  is  likely  caused  by  FOV's  completely  filled  with  deep  cumulus  convection.  Points  along 

the  denser  lines  are  clear  scenes  over  ocean  or  land  with  varying  fractions  of  convection  within 

the  scene.  Points  within  the  triangle  bounding  the  three  limiting  points  each  contain  some  mix  of 

land,  ocean,  clouds,  and  convection. 

Since  the  observed  TB  is  a  weighted  average  of  the  different  constituent  scenes  within  the 
FOV,  the  amount  of  convection  filling  a  scene  can  be  estimated  by  how  far  its  point  in  TB2-TBx 
space  is  pulled  away  from  the  non-scattering  line  and  toward  the  limiting  convective  point.  This 
is  essentially  what  the  ATB  from  the  simulations  in  Chapter  4  was  measuring.  By  measuring 
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how  far  the  high-freq  window  or  moisture  channel  TB  is  depressed  from  the  no-scatter  line,  we 
obtain  a  proxy  for  the  fraction  of  the  beam  filled  with  scattering  precipitation. 

The  data  samples  depicted  in  Figures  5.5  and  5.6  form  the  basis  for  the  following 


empirical  relations  for  the  predicted  high-frequency  non-scattering  TB  as  a  function  of  TB2: 


TB\5pred 


-  529.9  +  7.8234  x  TB2  -  0.01909  x  TB22 

270 +  —  (7B2-200) 

85 


[TB2 


TB2 < 200 K 
200K  <  TB2  <  285  K  ■ 
TB2 > 285 K 


TBHpred 


- 1387.6  + 19.68  x  TB2  -  0.05794X  TB22 
«  285 
TB2 


TB2  <  110K 
110K  <  TB2  <  285/sT  • 
TB2 > 285 K 


TB18pred=240K 

TB19prei=260K 

TB20pred=280K 


Sounding  channel  TB  depression  is  caused  by  rain  and  ice  that  fills  the  beam  both 
horizontally  and  vertically,  since  precipitation  that  doesn't  reach  the  level  of  sensitivity  of  a 


channel's  weighting  function  won't  affect  it  significantly.  So  each  channel  will  have  a  different 
slope  relating  its  scattering  ATB  to  that  of  a  high-frequency  window  or  moisture  channel.  Figure 


5.7  shows  scatter  plots  of  limb-corrected  TB7  and  TB8  vs.  ATB  15  and  ATB19.  Dashed  lines 
overlaid  on  the  plots  are  least-squares  best-fit  lines  for  points  with  high-frequency  depressions 
exceeding  one  standard  deviation.  As  predicted  by  the  simulation  results,  TB8  has  the  weakest 
scattering  reduction  signature.  The  bimodal  TB8  distribution  is  caused  by  seasonal  variation  in 
tropopause  height;  the  01  Jan  03  and  01  Apr  03  average  TB8  is  215.5  K  while  die  01  Jun  03  and 
01  Oct  03  average  TB8  is  216.5  K.  Table  5.1  summarizes  sounding  channel  vs. 
window/moisture  channel  depression  correlations. 
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Figure  5.7:  Limb-corrected  TB8  (top  row)  and  TB7  (bottom  row)  vs. 
TB15  scattering  depression  (left  column)  and  TB19  scattering 
depression  (right  column).  Dashed  lines  are  least-squares  best-fit  to 
points  with  ATB15  or  ATB19  <  -la. 


Table  5.1:  AMSU-A  limb-corrected  temperature  sounding 
channel  TB  correlation  with  window/moisture  channel  ATB 
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Figure  5.8:  Sounding  channel  ATB  vs.  ATB15  (left)  and  ATB19  (right)  for  2003 
Atlantic  TC  scenes.  Dots  denote  points  within  100  km  of  TC  center;  x's,  points  at  500 
km  from  TC  center;  and  +'s  1000  km. 
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The  nearly  linear  relationships  evident  among  limb-corrected  sounding  TB  and 
window/moisture  TB  depressions  should  make  it  straightforward  to  correct  for  scattering  by 
simply  adding  the  product  of  the  slope  and  the  window/moisture  channel  depression.  There  is  a 
complication,  however.  Figure  5.8  depicts  limb-corrected  sounding  channel  depression  vs.  TB15 
and  TB19  depression  at  three  different  radii  from  storm  center  in  the  2003  Atlantic  TC  scene 
data.  Points  at  500  km  and  1000  km  from  storm  center  have  nearly  identical  slopes  as  obtained 
from  the  tropical  ocean  scenes.  But  points  within  100  km  of  storm  center  appear  to  have 
significantly  larger  slopes,  especially  for  TB7  and  TB8.  For  all  radii,  there  is  a  large  cluster  of 
points  with  small  window/moisture  depression.  But  at  small  radii,  the  TC  warm  anomaly  is 
reflected  in  the  much  greater  spread  in  sounding  channel  TB.  It  is  this  population  which  appears 
to  have  the  greater  slope. 

There  are  at  least  two  possible  explanations  for  this  discrepancy.  The  first  is  simply  that, 
while  all  tropical  convection  is  weak  by  mid-latitude  standards,  it  is  slightly  stronger  in  TC 
eyewalls  and  spiral  bands  than  it  is  in  the  tropical  oceans  in  general.  Cecil  (2002)  used  the 
Tropical  Rainfall  Measuring  Mission  Microwave  Imager  and  Precipitation  Radar  to  show  that 
TC  eyewall  convection  tends  to  have  higher  radar  reflectivity  at  higher  altitudes  than  general 
oceanic  convection,  as  well  as  colder  85  and  37  GHz  polarization-corrected  temperature  (PCT; 
i.e.  greater  ice  scattering).  Figure  5.9  from  Cecil  (2002)  compares  the  radar  reflectivity  profiles 
in  different  tropical  convective  settings.  Saturated  parcels  start  their  ascent  in  either  the  TC 
eyewall  or  a  tropical  cumulonimbus  cloud  with  approximately  the  same  moisture  content-that  of 
a  saturated  parcel  at  a  sea  surface  temperature  near  300  K.  High-frequency  window  channels  are 
primarily  sensitive  to  the  vertically  integrated  water  and  ice  content-which  should  be 
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Figure  5.9:  Tropical  Rainfall  Measuring  Mission  (TRMM)  Precipitation  Radar  (PR) 
cumulative  reflectivity  profiles  for  (left  to  right)  TC  eyewalls,  TC  inner  rainbands,  TC  outer 
rainbands,  tropical  oceanic  convection,  and  tropical  continental  convection.  Bold  contours 
are  50th  and  99th  percentile  of  the  cumulative  distribution.  Dashed  line  is  overlaid  at  15  km 
alritude  (approximately  the  level  of  AMSU-A  TB8  weighting  function  peak).  From  Cecil 
(2002). 

approximately  the  same  in  either  eyewall  or  general  tropical  conditions— while  a  sounding 
channel  is  sensitive  only  to  the  scatterers  that  reach  the  level  to  which  its  weighting  function  is 
sensitive.  Said  differently,  the  mechanically-driven  ascent  in  the  TC  eyewall  may  be  able  to  loft 
the  same  liquid  and  ice  content  to  a  higher  altitude  than  the  weakly  buoyant  moist  ascent  in  a 
convective  hot  tower  in  the  inter-tropical  convergence  zone  (ITCZ),  thus  affecting  upper-level 
sounding  channels  more  for  the  same  window  channel  depression. 

A  second,  related,  explanation  is  based  on  the  difference  in  the  ambient  environment 
surrounding  the  convective  updrafts.  By  the  time  a  TC  forms,  the  radially  confined  heat  released 
by  the  organizing  convection  has  created  a  nearly  moist  lapse  rate  through  the  troposphere,  so 
convective  updrafts  in  the  eyewall  and  inner  spiral  bands  are  nearly  isothermal  with  the 
surrounding  environment  aloft.  Away  from  a  TC,  convective  updrafts  occur  in  an  unsaturated 
environment  where  parcels  are  radiating  to  space  and  sinking  subsiding  everywhere  outside  the 
convective  updraft  (thus  the  requirement  for  hot  towers  as  articulated  by  Simpson  et  al.  1998). 

So  the  mid  and  upper  levels  in  the  tropics  can  be  considerably  cooler  and  drier  than  the  warm, 
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saturated  environment  in  the  immediate  vicinity  of  a  TC  center.  The  convective  cell  affecting  a 
sounding  channel  is  averaged  with  the  TB  contribution  from  the  rest  of  the  3-D  volume  of  the 
FOV.  In  the  isolated  tropical  convective  cell,  the  updraft  core  may  be  substantially  warmer  than 
the  surrounding  environment,  so  the  TB  depression  caused  by  scattering  competes  with  warming 
contributed  by  the  updraft  itself.  For  convection  in  and  near  the  TC  eyewall,  though,  the 
updraft's  temperature  is  nearly  equal  to  its  surroundings,  so  scattering  can  depress  the  sounding 
channel  TB  without  competition  from  the  sensible  heat  added  by  an  isolated  hot  tower's  core. 
Figure  5.10  illustrates  conceptually  the  disparate  effects  of  scattering  precipitation  on  AMSU 
high-frequency  window  and  temperature  sounding  channels. 

Whatever  the  physical  mechanism,  it  appears  necessary  to  apply  a  larger  sounding  vs. 
window/moisture  TB  depression  slope  to  correct  for  scattering  in  proximity  to  a  TC  center. 


AMSU  TH7  FOV 


General 

Tropical 

Environment 


Net  effect: 
small  -  ATB 


TC  Eyewall 
Environment 


Net  effect: 
larger  -  ATB 
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Net  effect: 
~same  -  ATB 


Figure  5.10:  Conceptual  diagram  of  ambient  environment  modulation  of  sounding  channel 
scattering  TB  depression.  AMSU-A  sounding  channel  (e.g.  TB7)  FOV  (left  column)  and  high- 
frequency  window  channel  (e.g.  TB15)  FOV  (right  column)  for  isolated  tropical  convective  cells 
in  the  ambient  tropical  environment  (top  row)  and  eyewall  convection  (bottom  row).  Window 
channel  ambient  environment  (the  opaque  humid  lower  troposphere)  and  scattering  by 
hydrometeors  are  the  same  in  both  cases.  Sounding  channel  ambient  environment  is  cooler  for 
isolated  tropical  convection  so  warm  convective  core  counters  some  of  the  hydrometeor 
scattering  effect.  Sounding  channel  ambient  environment  and  convective  core  temperatures  are 
approximately  equal,  so  scattering  acts  unopposed.  Net  effect  is  larger  sounding  channel  ATB 
for  the  same  window  channel  ATB  in  TC  eyewall  convection. 
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Table  5.2  summarizes  the  slopes  for  the  general  tropical  convection  case  (weaker  slopes)  and  the 


Table  5.2:  Sounding  ATB  vs.  window/moisture  maximum  ATB  slope 
_ for  near-TC  cases  and  general  tropical  cases _ 


ATB  15 

ATB  17 

ATB  18 

ATB  19 

ATB20 

max 

general 

max 

general 

max 

general 

max 

general 

max 

general 

ATB  8 

.028 

.0060 

.020 

.0015 

.093 

.0029 

.030 

.0030 

.032 

.0012 

ATB7 

.058 

.013 

.029 

.0045 

.11 

.019 

.055 

.012 

.031 

.0058 

ATB6 

.070 

.032 

.072 

.012 

.13 

.038 

.042 

.017 

ATB5 

.086 

.067 

.043 

.025 

.14 

.069 

.065 

.048 

.052 

.031 

maximum  slope  observed  within  100  km  of  the  2003  TC  centers. 


The  correlation  coefficients  in  Table  5.1  indicate  that  ATB  18,  ATB  19,  and  ATB  15  are  best 


correlated  with  sounding  channel  depression.  These  three  predictor  channels  are  each  tested  for 
their  ability  to  improve  the  correlation  between  TB7  and  TB8  warm  anomaly  and  TC  MSLP.  For 
each  of  the  three  candidate  predictor  channels,  three  correction  methods  are  tested.  The  first 


correction  method  used  only  the  sounding  ATB  vs.  high-frequency  ATB  slope  derived  from  the 


non-TC  tropical  data  set.  The  second  method  used  the  maximum  slopes  from  Table  5.2  for 
FOV's  within  a  3  x  3  box  centered  on  a  TC's  estimated  location  and  the  weak  non-TC  slopes 


elsewhere.  The  third  method  used  the  average  of  the  weak  and  strong  slopes  for  the  3x3  storm- 
centered  box  and  weak  slopes  elsewhere. 

Testing  data  for  this  chapter  and  Chapter  6  was  a  set  of  497  AMSU  observations  of 
Atlantic  basin  TC's  from  1998  to  2004.  Each  AMSU  overflight  was  coincident  (within  +/-  3 
hours)  with  an  aircraft  reconnaissance  MSLP  observation,  which  serves  as  ground  truth.  For 
each  AMSU  overflight,  a  30  FOV  x  30  FOV  box  is  centered  on  the  storm's  position,  which  was 
interpolated  from  post-analyzed  best-track  data.  AMSU-B  TB's  were  convolved,  as  described  in 
Section  5.1,  and  stored  along  with  the  AMSU-A  TB's  for  the  30  x  30  FOV  array. 
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The  desired  scattering  correction  is  applied  (Figure  5.1 1  is  an  example  of  the  TBU-based 
scattering  correction  applied  to  limb-corrected  TB7  and  TBS)  and  the  TB7  and  TBS  warm 
anomalies  are  computed  as  follows: 

1.  Limb-correct  the  raw  sounding  channel  IB's 

2.  Compute  the  predicted  no-scatter  TB15  (or  TB18  or  TB19)  using  TB2 

3.  Compute  ATB15  (or  ATB18  or  ATB19) 

4.  If  ATB15  (or  ATB18  or  ATB19 )  is  less  than  zero,  multiply  it  by  the  appropriate  sounding 
channel  slope  and  add  the  absolute  value  of  the  result  to  the  limb-corrected  sounding 
channel  TB 

5.  For  each  FOV  adjacent  to  the  interpolated  storm  position: 

a.  Compute  the  environmental  brightness  temperature,  TBenv,  using  up  to  4  FOV's 
spaced  +/-  8  FOV's  in  the  along-  and  cross-track  directions  from  the  FOV  being 
considered 

b.  ATB7  or  ATB8  is  the  difference  between  the  selected  FOV’s  TB  and  TBenv. 

6.  Pick  the  largest  ATB7  or  ATB8  from  the  array  of  FOV's  surrounding  the  center. 

The  full  dataset  is  divided  into  a  249-sample  test  set,  used  to  develop  regression 
coefficients,  and  a  248-sample  validation  set.  Table  5.3  summarizes  correlations  between  ATB7 
or  ATB8  and  TC  MSLP  for  non-corrected  TB' s  and  each  of  the  nine  scattering  correction  methods 
tested. 


TableSJj^Con^arisonofTBTmidTOSAnornal^ConjelationwithTCMS^j^n^Z^S 


Scattering 

TB15 

TB18 

TB19 

correction 

method 

None 

general 

max 

average 

general 

max 

average 

general 

max 

average 

TB7  warm 
anomaly 
correlation 
TB8  warm 

.803 

.856 

.886 

.899 

.830 

.647 

.797 

.753 

.340 

.528 

anomaly 

correlation 

.872 

.889 

.899 

.903 

.875 

.646 

.810 

.865 

.569 

.720 

Surprisingly,  TB18-  and  7B79-based  scattering  corrections  reduced  the  correlations, 
dramatically  in  some  cases.  Variation  in  humidity  causes  considerable  variability  in  the  height  at 


Raw  AMSU-0  TB  eh;=17  Non- scatter- corrected  AMSU-A  TB  (K)  U_th  -  ?  Ch15  Scatter-corrected  AMSU-A  TB  (K)  k_ch 
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TB7  warm  anomaly  has  increased  by  approximately  0.5  K.  Scattering  effects  on  TB8  are  far  less  noticeable.  Note  also  that 
the  cool  surrounding  TB' s  (especially  in  the  northeast  spiral  band)  have  increased,  yielding  a  smoother  environmental  TB 
field. 
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Uncorrected  Warm  TB  (K) 


Uncorrected  Warm  TB  (K) 


Ch  7  Delta-TB  Scattering  Correction 


Ch  8  Delta-TB  Scattering  Correction 


Figure  5.12:  Channel  7  (left  column)  and  8  (right  column)  scattering  correction  effect  on 
environmental  TB  (top  row),  warm  TB  (middle  row),  and  ATB  (bottom  row)  for  248  Atlantic 
TC  cases  from  1998-2004. 
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which  the  TB18  and  TB19  weighting  functions  peak,  causing  a  wide  range  of  TB  variability  (20  - 
30  K)  not  attributable  to  scattering.  In  proximity  to  a  TC  center,  where  there  may  be  a 
considerable  area  of  near  saturation,  these  channels'  weighting  functions  would  rise  and  reduce 
TB18  and  TB19,  independent  of  any  scattering  effects.  This  may  be  causing  them  to  over-correct 
the  sounding  channel  TB's. 

TB15- based  scattering  corrections,  however,  improved  correlations  for  both  ATB7  and 
ATB8.  The  best  performance  is  obtained  using  the  averaged  weak  and  strong  slope  near  storm 
center  and  the  weak  correction  at  large  radii.  Since  the  scattering  correction  is  applied  to  all 
FOV's  in  the  30  x  30  array,  it  has  potential  to  affect  both  TBenv  and  the  warm  pixel  TB.  Figure 
5.12  depicts  pre-  and  post-correction  TBenv,  warm  TB,  and  ATB  for  channels  7  and  8.  Figure  5.13 


Figure  5.13:  Non-corrected  (left  column)  vs.  scattering-corrected  (right  column)  ATB7  (top 
row)  and  ATB 8  (bottom  row)  correlation  with  TC  MSLP  for  248  Atlantic  TC  cases  from 
1998-2004. 
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includes  four  scatter  plots  of  uncorrected  and  corrected  ATB7  and  ATB8  vs.  MSLP.  The 
reduction  in  scatter  is  especially  dramatic  for  ATB7. 

5.3.  Discussion 

The  results  presented  in  this  chapter  explain  the  physical  mechanism  underlying 
variability  in  TB2  and  high-frequency  window/moisture  channel  space.  They  demonstrate  that 
TB2  can  be  used  effectively  to  estimate  TB15  or  TB17  in  the  absence  of  precipitation  effects  (due 
to  the  TB18  through  TB20  variability  caused  by  humidity,  TB1  is  less  effective  at  pinning  down 
what  those  TB' s  should  be  in  the  absence  of  scattering).  ATB15,  ATB18,  and  ATB19  have  the 
most  robust  relationship  with  limb-corrected  sounding  channel  ATB,  and  of  those  three  channels, 
TBI5  is  the  only  one  that  improves  the  correlation  between  upper  tropospheric  warm  TB 
anomaly  and  TC  MSLP.  In  Chapter  6,  the  TB15  scattering  correction  method  will  be  applied  to 
CIMSS'  AMSU-based  TC  intensity  estimation  algorithm. 
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6.  A  Precipitation  Correction  for  the  CIMSS  AMSU-based  TC  Intensity 
Estimation  Algorithm 

Chapter  5  built  on  the  theoretical  basis  from  Chapter  3  and  the  simulation  study  results 
from  Chapter  4  to  develop  an  empirical  sounding  channel  precipitation  correction.  In  Section 
5.3,  it  was  shown  that  the  correction  substantially  improved  ATB7  and  slightly  improved  ATB8 
correlations  with  TC  MSLP.  This  chapter  will  apply  the  same  precipitation  correction  method  to 
the  current  CIMSS  AMSU-based  TC  intensity  algorithm  and  study  its  impact. 

6.1.  The  Current  CIMSS  Operational  TC  Intensity  Algorithm 

The  CIMSS  algorithm  is  one  of  two  AMSU-based  techniques  currently  providing  near 
real-time  intensity  estimates  to  TC  warning  centers.  The  other  is  the  CIRA  technique  described 
in  Section  2.3,  which  uses  gridded  AMSU  retrieved  temperature  fields. 

The  current  algorithm  builds  on  the  previous  work  of  Kidder  et  al.  (1978)  and  Velden  and 
Smith  (1983),  and  has  undergone  several  evolutionary  changes  since  its  inception.  In  its  original 
form,  the  algorithm  corrected  for  sub-sampling  using  Merrill's  (1995)  warm  anomaly  structure 
retrieval  initialized  with  a  storm  size  estimate  derived  from  the  Laplacian  of  AMSU-B  89  GHz 
TB's  (Brueske  2001;  Brueske  and  Velden  2003).  Kabat  (2002)  introduced  a  ATB8  regression 
relation,  in  addition  to  the  existing  ATB7  relation,  along  with  logic  governing  when  to  choose  a 
TB8  MSLP  estimate  over  a  TB7  estimate.  Herndon  et  al.  (2004)  replaced  the  warm  anomaly 
retrieval  and  TB16  size  estimate  by  using  an  external  TC  radius  of  maximum  wind  (RMW) 
estimate  from  a  TC  warning  center  as  a  proxy  for  the  horizontal  extent  of  the  warm  anomaly. 
Using  the  RMW  compared  to  the  known  FOV  size,  Herndon  et  al.  developed  an  MSLP  size  bias 


Ill 


Figure  6.1:  CIMSS  AMSU-based  TC  intensity  estimation  algorithm  flowchart.  The  dashed 
box  at  upper  right  indicates  the  point  at  which  precipitation  correction  would  take  place. 
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correction,  and  noting  that  the  warm  anomaly  vs.  MSLP  relationship  was  slightly  non-linear, 
introduced  two  regression  slopes  for  each  channel,  to  be  used  according  to  the  strength  of  the 
warm  anomaly.  Figure  1  is  a  flowchart  outlining  the  current  algorithm. 

The  two  most  significant  recent  changes  with  respect  to  precipitation  effects  were  the 
inclusion  of  ATB8  by  Kabat  (2002)  and  the  growth  of  its  influence  on  the  final  MSLP  estimate 
with  changes  implemented  by  Herndon  et  al.  (2004).  Currently  a  channel  7-based  MSLP 
estimate  is  only  chosen  in  weak  cases  ( ATB8  <  1.2  K)  where  channel  7  yields  a  lower  MSLP 
estimate  than  TB8.  As  noted  in  Chapters  4  and  5,  channel  8's  weighting  function  peaks  high 
enough  in  the  troposphere  as  to  make  it  only  weakly  susceptible  to  precipitation  effects,  and  to 
make  correction  problematic. 

Sources  of  error  in  the  current  algorithm  include:  instrument  noise  (tenths  of  a  K), 
bracketing  of  the  warm  core  between  adjacent  fields  of  view  (could  introduce  errors  comparable 
to  the  warm  anomaly  magnitude  itself);  sub-sampling  of  the  warm  core  by  the  large  instrument 
FOV,  especially  at  large  scan  angles  (again,  could  introduce  error  comparable  to  warm  anomaly 
magnitude),  and  precipitation  effects,  which  for  TB7  and  TB8  are  on  the  order  of  2  K  and  tenths 
of  a  K,  respectively.  The  presence  of  the  other  sources  of  error  imposes  a  limit  on  how  well  the 
algorithm  can  perform,  even  with  a  perfect  scattering  correction.  CIRA's  algorithm,  the  Vmax 
technique  described  in  Spencer  and  Braswell  (2001),  and  the  CIMSS  algorithm  have  all 
converged  to  very  nearly  the  same  accuracy  figure,  suggesting  that  each  technique  may  be 
squeezing  as  much  accuracy  out  of  the  current  instrument  as  possible. 
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6.2.  Precipitation  Correction  Methodology 

The  details  of  the  precipitation  method  were  discussed  in  Chapter  5,  but  are  outlined 


again  here: 

1.  TB2  is  used  to  predict  TB15  in  the  absence  of  scattering 


TBl5pred  = 


-  529.9  +  7.8234  TB2  -  0.019093  TB22  if  (TB2  <  200;  humid  scenes) 

270  + — (TB2  —  200)  if  (200  <  TB2  <  285;  cloudy  land  scenes ) 

85 

TB2  if  (TB2  >  285;  clear  land  scenes ) 


2.  If  TB15  is  depressed  below  this  predicted  value,  the  difference  is  multiplied  by  the 
empirical  sounding  ATB  vs.  ATB15  slope  found  in  Chapter  5  to  obtain  the  sounding 
channel  TB  correction. 

3.  For  AMSU  TB' s  within  2  FOV's  of  the  TC  center,  a  stronger  TB  correction  slope  is  used, 
which  averages  the  slope  found  Chapter  5  for  the  general  tropical  convection  cases  with 
the  slope  observed  in  TC  cases.  Table  6. 1  summarizes  these  sounding  TB  correction 
slopes  for  use  with  ATB15. 


Table  6.1:  Sounding  channel  TB  vs.  ATB  15 
precipitation  correction  slopes 


Channel 

General  slope 

TC  core  slope 

8 

.006 

.017 

7 

.013 

.036 

6 

.032 

.051 

5 

.067 

.075 

6.3.  Validation  data  set 

The  data  set  used  for  this  validation  effort  consists  of  497  AMSU  observations  of  TC's. 
Historical  storm  positions,  intensities,  and  RMW  at  six-hour  intervals  were  obtained  from  the 
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National  Hurricane  Center's  TC  Best  Track  dataset.  Using  archived  orbital  elements  for  the 
NOAA  satellites  and  the  best-track  storm  positions,  each  AMSU  overflight  of  a  known  TC  since 
June  1998  (when  the  first  AMSU  instrument  flew  on  NOAA- 15)  was  determined.  From  that 
large  listing  of  overflights,  times  within  3  hours  of  a  WP-3D  or  WC-130  aircraft  reconnaissance 
observation  of  the  storm  were  selected.  AMSU  data  for  the  selected  overflights  were  obtained 
from  the  AMSU  archive  at  CIMSS,  supplemented  by  data  obtained  from  NOAA's 
Comprehensive  Large- Array  Stewardship  System. 

Aircraft  reconnaissance  MSLP  is  taken  to  be  ground  truth  for  this  study,  but  it  has  several 
sources  of  error.  First,  the  MSLP  is  determined  either  from  a  dropsonde,  which  has  accuracy  of 
1  hPa  (Hock  and  Franklin  1999),  or  by  extrapolating  MSLP  from  the  temperature  and  pressure 
aloft,  which  can  have  substantially  larger  error.  Second,  the  sonde  is  dropped  as  nearly  as 
possible  to  where  the  crew  believes  the  storm  center  to  be,  but  this  position  may  be  in  error, 
causing  the  true  minimum  in  sea-level  pressure  to  be  missed.  Third,  since  the  AMSU  overflight 
and  reconnaissance  observation  are  not  exactly  coincident,  storm  intensity  change  during  the 
interim  may  also  add  error  on  the  order  of  a  few  hPa. 

The  precipitation-corrected  TB7  and  TB8  anomalies  will  have  a  different  MSLP 
relationship  than  their  uncorrected  counterparts,  so  the  MSLP  regression  coefficients  and  size 
bias  correction  coefficients  are  re-derived  after  scattering  correction.  To  avoid  over-fitting,  the 
data  set  is  divided  in  half,  with  the  249  odd-numbered  observations  used  to  obtain  the  regression 
coefficients  and  the  248  even-numbered  observations  used  for  validation. 
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6.4.  Cleaner  TB7  or  TB8  vs.  MSLP  Relationship  from  Scattering  Correction 

As  noted  in  Chapter  5,  scattering  correction  improves  the  correlation  between  TB7  or  TB8 
anomaly  and  MSLP,  so  before  studying  performance  of  the  current  CIMSS  algorithm,  we  note 
that  a  simple  linear  regression  of  TB7  anomaly  achieves  12.4  hPa  RMSE  (a  3.6  hPa 
improvement)  and  TB8  anomaly  achieves  12.2  RMSE  (a  1 .2  hPa  improvement).  As  Herndon  et 
al.  (2004)  noted,  AMSU  will  sub-sample  the  warm  anomaly  in  cases  where  the  RMW  is  small 
and  the  instrument  FOV  is  large;  this  becomes  a  significant  error  source  for  estimates  at  large 


Recon  MSLP  (hPa) 
♦  v*  No  scattering  correction 
ooo  Scattering  corrected 
—  Truth 


Error  (Estimate  -  Recon  MSLP,  hPa) 
No  scattering  correction 
Scattering  corrected 


TB8 


Recon  MSLP  (hPa) 


-  *  *  No  scattering  correction 
ooo  scattering  corrected 

—  Truth 


No  scattering  correction 
Scattering  corrected 


Figure  6.2:  Improvement  in  single-channel  size-bias  corrected  MSLP  estimates  after 
precipitation  correction  is  applied  to  AMSU  TB7  and  TB8.  RMSE  for  channel  7  estimates 
after  correction  is  9.3  hPa;  for  channel  8  it  is  9.0  hPa. 
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scan  angles.  If  a  multiple  regression  bias  correction,  using  RMW  and  FOV  is  applied  to  the 
single-channel  estimates,  RMSE  improves  to  9.3  hPa  (for  channel  7)  and  9.0  hPa  (for  channel  8). 
These  RMSE  values  are  comparable  to  non-scattering  corrected  RMS  errors  for  the  CIMSS 
algorithm.  Figure  6.2  includes  scatter  plots  of  AMSU  channels  7  and  8  MSLP  estimates  vs. 
aircraft-measured  MSLP  with  and  without  scattering  correction. 

6.5.  CIMSS  Algorithm  Performance  Improvement 

Figure  6.3  is  the  same  as  Figure  6.2,  but  for  the  CIMSS  two-channel  algorithm.  RMSE 
and  average  mean  error  improvement  (summarized  in  Table  6.2)  are  more  modest.  Two  reasons 
explain  why:  First,  as  noted  above,  as  the  algorithm  has  evolved,  it  has  become  very  TB8 
centric.  TB7  estimates  only  come  into  play  for  weak  TC's,  and  even  then,  only  when  channel  7 
estimates  a  lower  MSLP  than  channel  8  (which  sometimes  occurs  when  the  nascent  warm 
anomaly  is  confined  to  around  250  hPa  and  hasn't  yet  grown  upward  to  150  hPa).  As  the  results 
in  Chapter  4  and  5  indicate,  AMSU  channel  8  is  largely  immune  to  scattering  by  tropical 
convection,  which  rarely  reaches  the  15  km  height  where  the  limb-corrected  channel  8  weighting 
function  peaks. 

The  second  reason  for  the  modest  CIMSS  algorithm  performance  improvement  is  that  the 
current  size  bias  correction  is  already  indirectly  correcting  for  some  scattering  effects.  Herndon 
et  al.  (2004)  found  a  weak  bias  for  small  storms  and  large  scan  angles  (where  the  AMSU  FOV  is 
largest)  caused  by  sub-sampling.  So  when  the  algorithm  determines  that  the  warm  anomaly  size 
(as  represented  by  RMW)  is  smaller  than  the  instrument  FOV  size,  it  decreases  the  MSLP 
estimate  to  correct  the  weak  bias  caused  by  the  sub-sampled  warm  anomaly. 


Current  CIMSS  Algorithm 


Recon  MSLP  (hPa) 
*  *  •  No  scattering  correction 
ooo  Scattering  corrected 
—  Truth 


Figure  6.3:  Comparison  of  CIMSS  TC  intensity  algorithm  performance 
before  and  after  precipitation  correction.  The  histogram  at  bottom  compares 
errors  without  (background)  and  with  (foreground)  scattering  correction. 


Table  6.2:  Comparison  of  RMSE,  Average  Mean  Error 
(AME),  and  bias  of  CIMSS  TC  intensity  estimates  with  and 


without  precipitation  correction 


Method 

RMSE  (hPa) 

AME  (hPa) 

bias  (hPa) 

Non-corrected 

8.8 

6.8 

-0.01 

Precipitation 

corrected 

8.0 

6.1 

1.1 

n  =  248 

Figure  6.4  is  a  scatter  plot  of  RMW  vs.  MSLP  for  the  TC's  in  the  n  =  248  validation 
dataset  tested  here,  and  it  shows  the  strong  relationship  between  increasing  TC  intensity  and 
decreasing  storm  radius.  But  as  storm  intensity  increases  (and  intense  convection  is  localized  in 


the  relatively  small  eyewall),  precipitation  scattering  is  more  likely  to  significantly  weaken 
sounding  channel  TB's  near  storm  center.  So  sub-sampling  and  precipitation  scattering 
complement  each  other  and  both  conspire  to  weaken  the  warm  anomaly  measured  by  the 
instrument,  weakening  the  AMSU  intensity  estimate.  Thus,  an  empirical  bias  correction  that 
corrects  for  weak  estimates  when  RMW  is  small  will  also  remove  some  of  the  error  caused  by 


RMW  vs.  MSLP 


Recon  MSLP  (hPa) 

Figure  6.4:  TC  radius  of  maximum  wind  vs. 
MSLP  for  n  =  248  validation  sample. 


precipitation  scattering.  As  evidence  that  this  is 
in  fact  occurring,  Table  6.3  compares  the 
RMSE  improvement  of  the  size  bias  correction 
for  non-precipitation  corrected  vs.  precipitation 
corrected  TB7  and  TBS  MSLP  estimates.  The 
size  bias  contribution  is  reduced  by  about  1/3 
(for  channel  8)  and  about  1/2  (for  channel  7). 
This  is  significant  because  the  RMW  value  is 


Table  6.3:  Comparison  of  size-bias  RMSE  improvement  on  non- 
corrected  and  precipitation-corrected  AMSU  TC  intensity  estimates 


Size  bias 

Pre  size  bias  Post  size  bias  RMSE 

correction  correction  improvement 

RMSE  (hPa)  RMSE  (hPa)  (hPa) 


9.9  8.8  1.1 

8.8  8.0  0.8 


Method _ 

Non-corrected 
Precipitation  corrected 
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obtained  from  the  TC  warning  centers,  where  it  is  estimated  by  satellite  analysts,  and  is  therefore 
subject  to  significant  inaccuracy,  and  often  is  missing  due  to  high  analyst  workload  or  neglect. 

Table  6.4  compares  RMSE,  AME,  and  bias  changes  vs.  TC  intensity  before  and  after 
precipitation  correction.  Performance  for  950-975  hPa  storms  actually  deteriorated  slightly  due 
to  creation  of  some  outliers.  The  regression  fit  for  weak  storms  and  very  strong  storms,  however 
was  improved. 


Table  6.4:  CIMSS  Algorithm  RMSE,  AME,  and  bias  by  TC  intensity 


RMSE 

AME 

bias 

No  precipitation 

>1000  hPa 

7.6 

5.9 

-1.3 

correction 

975-1000  hPa 

8.2 

6.3 

-2.5 

950-975  hPa 

9.5 

7.5 

0.1 

<950  hPa 

10.5 

8.3 

5.7 

Precipitation 

>1000  hPa 

5.7 

4.3 

-0.9 

corrected 

975-1000  hPa 

7.2 

5.9 

0.0 

950-975  hPa 

10.0 

7.5 

2.9 

<950  hPa 

9.9 

7.8 

4.2 

n=248 

6.5.  Another  Approach  Enabled  by  Scattering  Correction 

In  Section  2.1,  the  expression  relating  MSLP  to  temperature  anomaly  requires  the 
assumption  that  the  TB7  or  TB8  anomaly  represents  the  mean  temperature  anomaly  through  the 
full  depth  of  the  troposphere.  This  assumption  was  necessary  because  of  precipitation  effects  on 
the  low-  and  mid-tropospheric  sounding  channels  ( TB5  and  TB6).  Surface  emissivity  also  affects 
the  lowest  channel,  TB4.  This  precipitation  correction  does  not  correct  TB5  or  TB6  enough  to 
add  additional  information  beyond  what  the  precipitation  corrected  TB7  or  TB8  currently 
contains,  but  the  corrected  TB's  are  usable  for  producing  a  vertically-averaged  TB  anomaly, 
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which  has  comparable  RMSE  and  AME  (9.4 
hPa  and  7.1  hPa)  to  single-channel  estimates. 
An  area  for  future  work  is  to  refine  this 
vertically-integrated  technique  and  explore 
whether  it  has  skill  in  situations  where  single¬ 
channel  estimates  are  prone  to  error. 

6.6.  Discussion 

Correcting  for  precipitation 
substantially  improves  the  TB  anomaly  vs. 
MSLP  relationship  for  channels  5  through  7 
while  slightly  improving  channel  8.  The 
current  CIMSS  algorithm,  by  evolving  to 
become  more  reliant  on  the  TB8  anomaly  and 
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employing  a  size  bias  correction  that 
indirectly  corrects  for  scattering-induced  weak 
intensity  estimates,  had  already  reduced  the 


Figure  6.5:  Precipitation  correction 
improvement  in  a  TC  intensity  estimation 
technique  using  a  vertically-integrated  TB 
anomaly  employing  channels  5  through  8. 


impact  of  precipitation  on  intensity  estimates.  Essentially  all  of  the  gross  outliers  in  the  data  set 


tested  resulted  from  small  storm  size  and/or  scan  geometry  effects  (sub-sampling  and 
bracketing).  The  primary  benefit  of  the  scattering  correction  is  to  reduce  some  of  the  spread 
across  the  spectrum  of  intensity,  rather  than  correcting  individual  gross  outliers.  A  secondary 
benefit  is  the  reduction  of  dependence  on  the  external  RMW  estimates. 
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Most  significantly,  this  correction  method  is  the  most  physical  of  the  methods  employed 
to  date.  Velden  and  Smith  (1983),  Velden  (1989),  and  Velden  et  al.  (1991)  manually  removed 
FOV's  that  qualitatively  appeared  to  be  precipitation-affected.  Spencer  and  Braswell  (2001) 
included  precipitation-sensitive  channels  as  multiple-regression  predictors,  but  this  method  does 
not  enable  a  true  determination  of  the  unaffected  TB  gradient,  the  true  physical  link  to  Vmax. 
Demuth  et  al.  (2004)  employ  the  most  physical  method  to  date,  using  AMSU  cloud  liquid  water 
and  cloud  ice  water  correlation  with  reduced  tropospheric  sounding  channel  TB.  But  as 
discussed  in  Chapter  3,  Li  and  Weng  (2002)  point  out  that  this  type  of  correction  is  effective 
mostly  for  the  lower-tropospheric  sounding  channels  (TB4  and  7715).  DeMuth  et  al.  (2004)  are 
left  to  correct  the  upper-tropospheric  sounding  channels  by  smoothing  the  Laplacian  of  the  TB 
field  when  they  detect  an  anomalously  cold  TB.  The  problem  here  is  that  a  warm  anomaly  FOV, 
which  is  still  warmer  than  its  surroundings,  but  not  as  warm  as  it  should  be  in  the  absence  of 
precipitation  scattering,  will  not  be  corrected. 

The  precipitation  correction  method  here  uses  AMSU  channels  insensitive  to 
tropospheric  temperature  to  quantify  precipitation  effects  on  sounding  channel  TB’s,  then 
corrects  those  TB's  before  applying  them  to  a  TC  intensity  algorithm.  As  a  result,  this  technique 
can  be  applied  to  any  chosen  algorithm  relating  sounding  TB  to  TC  MSLP,  and  even  other 


similar  instruments. 
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7.  Conclusions  and  Future  Work 

7.1.  Summary 

Hurricanes  are  significant  natural  disasters,  so  it’s  important  to  forecast  their  impacts 
correctly  for  both  civil  society  and  military  operations.  Forecasting  TC  impacts  requires 
accuracy  in  forecasting  TC  track,  TC  maximum  sustained  winds,  and  the  TC  wind  distribution. 
Recently  TC  track  forecasting  has  improved  significantly,  while  intensity  forecasting  has 
stagnated  because  the  processes  governing  it  are  still  poorly  understood  or  not  resolvable  by 
existing  NWP  models.  The  best  forecast  tools  are  statistical,  relying  on  current  TC  intensity  and 
intensity  trend.  In  situ  TC  intensity  obs  are  very  rare,  so  the  bulk  of  intensity  measurements  are 
remotely-sensed  by  satellite.  The  mainstay  Dvorak  technique,  which  uses  visible  and  IR 
geostationary  imagery,  can  be  subjective,  labor  intensive,  and  constrained  by  rules  that  limit  its 
ability  to  respond  to  rapid  storm  changes.  Complimenting  this  remote  sensing  technique  with 
one  that  is  more  physical  and  objective  would  benefit  forecasters.  In  this  spirit,  microwave 
temperature  sounding  instruments  have  been  used  for  nearly  30  years  to  estimate  TC  MSLP  by 
measuring  warm  anomaly  in  upper-tropospheric  TB  above  storm  center.  Until  the  AMSU 
instrument,  this  wasn't  feasible  operationally  because  of  poor  instrument  resolution  and  data 
refresh  rates.  Today  CIMSS'  AMSU  technique  performs  on  a  par  with  the  Dvorak  technique,  but 
still  suffers  from  uncertainty,  part  of  which  is  due  to  scattering  by  hydrometeors  in  eyewall  or 
spiral  band  convection.  This  study  used  theory  and  observational  data  to  develop  an  empirical 
correction  to  reduction  in  tropospheric  sounding  channel  brightness  temperatures  (7B's)  caused 
by  scattering. 
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Chapter  3  showed  that  AMSU-A's  low  frequency  window  channels  were  sensitive  to 
relative  humidity  changes  and  rain  emission.  High  frequency  channels  were  subject  to  rain 
emission,  but  also  scattering  by  large  liquid  and  frozen  hydrometeors,  especially  graupel. 

Because  the  atmosphere  is  transparent  in  the  window  bands  between  O2  and  water  vapor 
absorption  lines,  these  channels  are  sensitive  to  emission/scattering  caused  by  the  total  optical 
depth  of  hydrometeors  throughout  the  depth  of  the  troposphere,  while  sounding  channels  are 
only  sensitive  to  hydrometeors  impinging  the  channel's  weighting  function.  Nevertheless,  a 
simple  1-d  plane-parallel  conceptual  model  showed  that  window  channel  TB  perturbations 
caused  by  a  rain  and  ice  layer  correlated  with  ATB  in  an  upper-tropospheric  sounding  channel. 

Chapter  4  employed  a  unique,  polarized  reverse  Monte  Carlo  radiative  transfer  model, 
coupled  with  prescribed  profiles  of  temperature,  pressure,  humidity,  cloud  liquid,  rain,  snow,  and 
graupel  (some  based  on  Petty  1-d  rain  model;  some  obtained  from  a  fine-scale  model  simulation 
of  1998's  Hurricane  Bonnie)  to  verify  under  more  realistic  conditions  the  utility  of  the  high- 
frequency  window  TB' s  for  correcting  sounding  channel  TB' s.  AMSU-A  channel  2  TB  was 
found  to  be  a  good  predictor  of  high-frequency  window  channels  ( TB15  and  TBIT)  in  the 
absence  of  scattering.  The  depression  of  high-frequency  TB' s  from  their  modeled  no-scattering 
values  correlated  well  with  depression  in  tropospheric  temperature  sounding  channels. 

Individual  high-frequency  channels,  particularly  AMSU-B  moisture  channels,  performed  as  well 
as  the  first  principal  component  of  the  channel  ensemble. 

Chapter  5  analyzed  four  days  of  N15,  N16,  and  N17  AMSU  data  over  oceans  between 
15°  S  and  15°  N  latitude.  Using  this  data,  and  a  compilation  of  125  TC  scenes  from  2003,  an 
empirical  prediction  for  no-scatter  TB15  through  TB20  TB's  using  TB2  was  developed.  The 
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depression  from  the  predicted  no-scatter  high-frequency  TB  values  was  correlated  with  sounding 
channel  TB  depressions.  The  slope  relating  window  and  sounding  channel  ATB,  however,  was 
found  to  be  stronger  within  100  km  of  a  TC  center.  Two  possible  explanations  are  offered,  both 
centered  on  difference  between  convective  clouds  in  the  inter- tropical  convergence  zone  (ITCZ), 
where  the  slopes  were  derived  versus  in  the  TC  eyewall  or  inner  rainbands.  The  first  is  based  on 
Cecil's  (2002)  finding  that  eyewall  and  rainband  convection  produce  greater  radar  reflectivity 
aloft  than  ITCZ  convection.  Since  both  types  of  convection  start  from  similar  near-saturated 
surface  parcels  at  approximately  the  same  temperature  as  the  sea  surface,  there  may  not  be  much 
vertically-integrated  ice  or  water  content  difference  between  the  two  types,  and  so  little  window- 
channel  TB  difference,  but  the  greater  reflectivity  aloft  indicates  hydrometeors  are  lifted  higher 
in  eyewall  convection,  causing  greater  effects  on  upper-tropospheric  temperature  sounding 
channels. 

The  second  explanation  is  based  on  the  different  ambient  thermodynamic  environment  of 
the  two  settings.  The  warming  in  an  ITCZ  hot  tower  interior  competes  against  scattering  to 
influence  microwave  TB,  reducing  the  scattering  effect.  Eyewall  convection  is  occurring  in  a 
warm,  saturated  environment  and  is  mostly  mechanically  driven.  So  scattering  is  free  to  reduce 
the  microwave  TB  without  interference  from  warming  within  cloud.  Both  processes  might  be 
both  at  work  simultaneously.  The  net  result  is  that  for  a  given  high-frequency  TB  depression, 
hydrometeors  in  a  TC  eyewall  or  inner  spiral  band  are  more  effective  at  reducing  sounding 
channel  TB.  Empirical  scattering  corrections,  using  different  slopes  near  TC  center  vs.  at  larger 
radii,  were  tested,  and  TB15  was  found  to  perform  best  in  reducing  the  scatter  between  channel  7 
and  8  warm  anomaly  magnitude  and  TC  MSLP. 
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In  chapter  6,  an  n  =  497  sample  of  AMSU  overflights  of  TC's  between  1998  and  2004 
was  used  to  validate  the  performance  of  the  scattering  correction  when  applied  to  the  current 
CIMSS  TC  intensity  technique.  Limb-corrected  sounding  channel  TB's  were  first  corrected  for 
scattering  using  the  empirical  method  developed  in  Chapter  5.  Applying  this  correction 
improved  CIMSS  algorithm  RMSE  and  average  mean  error  by  about  10%.  Part  of  the  reason  for 
the  modest  improvement  is  the  CIMSS  algorithm's  strong  dependence  on  the  TB8  warm 
anomaly,  which  (at  150  hPa)  is  largely  immune  from  scattering  effects.  Additionally,  the 
algorithm  is  already  indirectly  correcting  for  scattering  through  its  size  bias  correction.  At  small 
storm  sizes,  sub-sampling  reduces  the  AMSU  field  of  view's  sensitivity  to  the  warming  present. 
At  the  same  time  more  convection  is  occurring  (since  strong  storms  tend  to  be  smaller)  near  the 
center  and  weakening  observed  TB  anomaly.  Both  effects  reduce  the  observed  warm  anomaly, 
so  correcting  for  the  observed  weak  bias  for  small  storms  corrects  for  scattering  at  the  same  time 
that  it  corrects  for  sub-sampling. 

Correcting  all  tropospheric  sounding  channels  for  precipitation  effects  improves  the 
lower  channels'  accuracy  to  the  extent  that  a  vertically-averaged  TB  anomaly  now  performs  on  a 
par  with  single-channel  TB7  or  TB8  and  may  do  a  better  job  of  fully  sampling  the  deep-layer 
thermal  anomaly  in  future  intensity  techniques. 

7.2.  Future  Work:  Adaptation  to  SSMIS 

The  next  generation  of  operational  temperature  sounding  instrument,  the  Special  Sensor 
Microwave  Imager/Sounder  (SSMIS)  was  launched  in  Fall  2003  and  is  currently  completing  its 
calibration/validation  period.  This  instrument  combines  the  conical  scan  geometry  and  high 
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resolution  of  the  earlier  SSMI  instrument,  with  the  temperature  and  moisture  sounding  channels 
of  the  AMSU  instrument. 

SSMIS'  conical  scan  geometry  will  benefit  TC  intensity  estimation  because  limb- 
correction  and  variable  field-of-view  sizes  across  a  scan  line  will  no  longer  be  sources  of  error. 

Its  higher  resolution  (37  km  vs.  48  to  150  km  for  AMSU-A)  will  help  to  reduce  sub-sampling 
problems. 

But  SSMIS  has  potential  to  be  more  vulnerable  to  scattering.  Its  suite  of  temperature 
sounding  channels  includes  equivalents  to  AMSU-A's  channel  6  and  channel  8,  but  not  channel 
7.  Channel  8  is  the  ideal  MSLP  estimation  channel  for  AMSU,  since  it  is  relatively  immune 
from  precipitation  effects  while  retaining  sensitivity  to  the  TC  warm  anomaly.  The  key,  though, 
is  that  at  large  scan  angles,  where  temperature  sounding  channel  weighting  functions  move  to 
higher  levels  due  to  limb  darkening,  AMSU  has  channel  7  available  to  blend  with  channel  8  via 
limb  correction,  and  keep  the  limb-corrected  channel  8  sensitive  to  150  hPa  temperature  across 
the  full  scan.  SSMIS'  scan  angle  is  fixed  at  45°,  where  its  channel  8  weighting  function  will  peak 
significantly  higher  than  150  hPa.  Its  next-lower  channel  is  equivalent  to  AMSU-A  channel  6, 
whose  weighting  function  peaks  at  approximately  250  hPa  at  45°  scan  angle.  This  is 
approximately  the  level  to  which  AMSU's  limb-corrected  channel  7  is  sensitive  and  precipitation 
effects  are  much  more  pronounced  than  they  are  at  150  hPa. 

This  means  that  SSMIS  will  likely  rely  on  its  channel  6  equivalent  to  estimate  TC 
intensity,  and  will  need  scattering  correction  to  achieve  accuracy  equivalent  to  the  AMSU-based 
technique.  Fortunately,  SSMIS  also  has  low-  and  high-frequency  window  channels  similar  to 
AMSU.  The  precipitation  correction  technique  developed  in  this  study— temperature  sounding 
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channel  correction  via  depression  of  a  high-frequency  window  channel  from  its  low-frequency 
derived  no-scatter  predicted  value— will  be  directly  applicable,  and  potentially  more  vital  to 
producing  useful  TC  intensities  using  this  future  generation  of  instrument. 
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Appendix:  A  Polarized  Reverse  Monte  Carlo  Radiative  Transfer  Model 

A.l.  Motivation 

This  appendix  describes  in  detail  the  polarized  reverse  Monte  Carlo  (RMC)  radiative 
transfer  model  (RTM)  used  to  accomplish  the  simulation  study  in  Chapter  4.  The  purpose  of  that 
chapter  was  to  determine  whether  scattering-induced  depression  in  high-frequency  Advanced 
Microwave  Sounding  Unit  (AMSU)  brightness  temperatures  (TB)  could  be  used  to  quantify 
scattering  effects  on  AMSU  temperature  sounding  TB.  This  was  accomplished  using  an 
ensemble  of  19  one-dimensional  profiles  with  known  temperature,  moisture,  cloud  liquid,  rain, 
snow,  and  graupel.  Nine  of  those  profiles  were  prescribed  using  a  1-D  model  developed  by  Petty 
(2001a).  The  remaining  ten  were  obtained  from  a  Mesoscale  Model  5  (MM5)  simulation  of 
Hurricane  Bonnie  (1998).  By  comparing  modeled  TB's  at  AMSU  frequencies  and  scan 
geometries  with  precipitation  included  and  removed,  the  precipitation  effect  on  each  channel  was 
determined  for  use  in  developing  an  empirical  scattering  correction. 

The  requirement  for  the  study,  then,  was  an  RTM  that  was  flexible  enough  to  allow 
different  forms  of  meteorological  variable  inputs,  would  easily  allow  changes  or  removal  of 
individual  moisture  or  precipitation  components  to  isolate  their  effects,  would  run  quickly 
enough  to  allow  multiple  cases  to  be  completed  relatively  quickly,  would  allow  a  number  of 
different  types  of  output.  A  large  number  of  existing  radiative  transfer  codes  exist  in  widespread 
use,  but  each  would  require  considerable  time  to  learn  and  adapt  to  this  purpose.  Additionally, 
the  author  wanted  to  avoid  depending  results  from  a  "black  box,"  of  whose  inner  workings  he 
had  limited  understanding.  Finally,  future  work  with  this  model  will  include  adaptation  to  3-D 
inhomogeneous  scenes.  Other  radiative  transfer  methods— two-stream,  successive  orders  of 
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scattering,  discrete  ordinates-are  not  readily  adaptable  to  geometries  other  than  1-D  plane- 
parallel.  Only  Monte  Carlo  methods,  as  described  in  section  A.2,  are  well-suited  for  this 
geometry.  Other  Monte  Carlo  microwave  radiative  transfer  codes  have  been  developed,  but  all 
are  developed  with  the  Special  Sensor  Microwave  Imager  (SSMI)  in  mind  and  are  limited  to  the 
special  cases  of  vertical,  horizontal,  or  no  polarization,  and  fixed  scan  angles.  For  these  reasons, 
the  author  opted  to  develop  a  new  RMC  RTM  able  to  handle  the  large  number  of  AMSU-A  and  - 
B  scan  angles,  frequencies,  and  polarizations,  accept  different  types  of  meteorological  profiles, 
allow  easy  exclusion  of  some  or  all  precipitation  constituents,  and  readily  adaptable  to  3-D 
geometry. 

A.2.  Previous  Monte  Carlo  Models 

Petty  (1994,  1994a)  was  the  first  to  introduce  RMC  techniques  to  the  application  of 
microwave  radiative  transfer  in  the  atmosphere.  He  described  how  a  photon  (or  packet  of 
photons)  traveling  from  an  emitting  source,  through  a  scattering  medium,  to  a  sensor  could  be 
just  as  easily  traced  backward  along  the  same  path  from  sensor  to  emitter.  The  reciprocity 
theorem  allows  the  same  scattering  phase  function  that  predicts  the  probability  of  a  forward¬ 
traveling  photon  scattering  into  a  given  direction  to  be  used  for  determining  the  probability  that  a 
backward-traveling  photon  was  scattered  from  a  given  direction.  Random  numbers  are  used:  1) 
in  conjunction  with  the  local  extinction  coefficient  to  generate  a  photon's  path  length  between 
extinction  events;  and  2)  to  determine  the  direction  of  scattering  when  a  scattering  event  occurs. 
Petty  gave  each  photon  an  initial  weight  (visualized  as  a  probability  of  continuing  to  exist)  and 
removed  a  fraction  of  the  weight  equal  to  (1  -  Gf)  at  each  extinction  event.  This  continues  for  the 
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photon  through  numerous  scattering  events  until  its  weight  is  reduced  below  some  threshold  (at 
which  point  it  is  considered  totally  absorbed)  or  it  leaves  the  domain.  At  each  point  where  a 
photon  encounters  extinction,  the  probability  density  it  sheds  in  the  event  is  accumulated  at  that 
point.  After  all  the  photon  packets  have  been  exhausted,  the  relative  weight  at  each  point  in  the 
3-D  domain  determines  how  much  the  physical  temperature  at  that  point  contributes  to  the  TB 
observed  by  the  sensor.  Petty  handled  polarization  by  treating  each  photon  packet  as 
permanently  polarized,  either  vertically  or  horizontally,  and  using  the  appropriate  polarized 
reflectivity  upon  encountering  the  ocean  surface.  Necessary  inputs  to  the  model  are  3-d 
domains  of  temperature,  extinction  coefficient,  single-scatter  albedo,  and  asymmetry  parameter. 

Petty  used  his  results  to  demonstrate  the  feasibility  of  RMC  methods  for  predicting 
realistic  microwave  TB's  in  highly  inhomogeneous  3-D  domains.  He  also  noted  two  phenomena 
that  would  be  impossible  to  model  in  plane-parallel  geometry:  First,  a  significant  TB 
contribution  comes  from  warm  rain  emission  observed  from  the  sides  of  clouds,  reducing  the 
scattering  impact  of  the  ice  layer  in  the  top  of  the  cloud,  which  would  always  be  encountered  if 
plane-parallel  geometry  was  used.  Second,  ocean-surface  reflection  of  warm  rain  emission  can 
contribute  significantly— and  increase— the  observed  TB. 

Roberti  et  al.  (1994)  also  developed  a  RMC  technique  at  about  the  same  time,  and  present 
a  detailed  discussion  of  their  technique.  Liu  (1996)  developed  a  novel  polarization  treatment  in 
his  backward-forward  Monte  Carlo  (BFMC)  model.  He  traced  photons  backward  unpolarized  to 
their  emission  sources,  then  forward  through  the  same  scattering  points,  as  a  pair  of  photons,  one 
vertically  polarized  and  one  horizontal.  At  each  scattering  event,  he  used  a  rotated  phase  matrix, 
which  would  produce  different  intensities  for  the  vertically  and  horizontally  polarized  incident 
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photons,  to  account  for  non-isotropic  scattering.  The  end-state  polarization  for  the  photon  was 
determined  from  the  ratio  of  the  vertical/horizontal  pair’s  intensities.  Roberti  and  Kummerow 
(1999)  concluded,  instead,  that  rigorous  polarization  in  a  Monte  Carlo  model  is  only  possible  in 
the  forward  direction.  Kim  et  al.  (2004)  compared  various  radiative  transfer  methods,  using  an 
RMC  technique  as  the  accuracy  benchmark. 


A.3.  Key  Assumptions  and  Approximations 

Several  assumptions  influence  the  results  from  this  RMC  implementation.  First,  the 
assumed  form  of  the  liquid  and  frozen  precipitation  size  distribution  determines  how  a  given 
precipitation  density  affects  the  Mie  properties  of  the  medium.  This  model  uses  the  Marshall- 
Palmer  (MP)  distribution  for  rain  and  graupel,  modified  per  Petty  (2001a)  to  account  for  fall- 
speed  variation  with  density,  and  its  effect  on  relating  precipitation  density  to  rain  rate.  The 
Sekhon-Srivastava  (SS)  distribution  is  used  for  snow,  with  the  exponential  decay  parameter, 
A(R)  scaled  by  2'2,  to  produce  more  realistic  snow  size  distributions,  again  per  Petty  (2001a). 

The  composition  of  graupel  will  affect  its  Mie  parameters  strongly,  by  changing  its 
complex  refractive  index,  as  well  as  its  geometric  size.  For  the  study  in  Chapter  4  were  obtained 
following  Petty's  methodology:  above  the  freezing  level,  graupel  is  70%  air  and  30%  ice.  When 
it  falls  below  the  freezing  level,  it  melts  smoothly  over  0.5  km  of  descent,  with  its  water  fraction 
increasing  linearly  from  0%  to  100%  with  the  remainder  composed  of  the  70%  air-30%  ice 
mixture.  The  complex  refractive  index  for  all  cases  is  computed  using  the  Bruggeman  formula: 
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where  sj  and  £2  are  the  dielectric  constants  of  two  mixed  substances,/;  is  the  fraction  of  the  first 
substance,  and  sav  is  the  dielectric  constant  of  the  mixture  (Bohren  and  Huffman  1983;  Petty 
2004).  For  melting  graupel,  sav  for  the  air-ice  mixture  is  again  averaged,  this  time  with  the 
appropriate  fraction  of  water. 

Mie  theory,  of  course,  introduces  a  very  important  assumption-spherical  particles. 
Aerodynamic  forces  make  falling  raindrops  oblate  and  give  them  a  preferred  orientation  and 
riming  and  accretion  of  frozen  particles  may  produce  pronounced  asymmetries.  Treatment  of 
such  non-spherical  particles  and  their  preferred  orientations  is  problematic,  and  well  beyond  the 
scope  of  this  study,  but  it  needs  to  be  noted  as  a  source  of  uncertainty. 

Finally,  polarization  is  handled  in  a  limited  fashion  in  this  implementation.  Photons  are 
initialized  with  a  linear  polarization  matching  the  polarization  angle  of  the  sensor,  which  rotates 
with  scan  angle.  A  rigorous  treatment  of  polarization  would  require  using  the  Mie  properties  at 
each  grid  point  to  compute  the  4x4  Mueller  matrix  elements  (Section  3.1),  then  using  that  matrix 
to  determine  the  probability  of  an  incident  polarized  photon  scattering  in  any  direction,  and  the 
polarization  it  would  assume  at  each  direction.  The  scattered  direction  would  then  be  chosen 
randomly  according  to  the  directional  probability  distribution.  This  is  computationally 
prohibitive,  particularly  in  3  dimensions. 

Instead,  photons  are  limited  to  linear  polarization.  At  each  level  in  the  1 -dimensional 
domain,  Mie  theory  is  used  to  compute  the  upper-left  diagonal  elements  of  the  Mueller  matrix. 
The  scattered  intensity  is  then  computed  for  each  scattering  angle,  6.  This  is  repeated  after 
rotating  the  incident  photon  through  0°  to  90°  about  the  incident  direction.  The  result  is  a  two- 
dimensional  array  of  scattered  intensities  for  all  combinations  of  scattering  angle  and  rotation 
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angle.  This  array  is  stored  for  each  frequency,  channel,  and  model  grid  point  before  beginning 
the  RMC  radiative  transfer. 

When  a  photon  is  scattered  by  the  RMC  model,  the  scattering  direction  is  randomly 
chosen  using  the  stored  probability  density  function.  The  incident  and  scattering  directions 
together  form  the  plane  of  scattering,  and  the  incident  photon's  polarization  state  is  computed 
relative  to  this  plane,  and  preserved  after  it  is  scattered.  A  rigorous  treatment  of  scattering,  by 
use  of  the  full  Mueller  matrix  would  allow  obliquely  polarized  incident  photons  to  produce 
elliptically  polarized  scattered  photons  and  vice  versa  (via  the  lower-right  comer  Mueller  matrix 
elements).  This  treatment  of  scattering  ignores  that  aspect.  However,  the  scattered  polarization 
state  is  of  third  order  importance;  of  first  order  is  whether  scattering  occurs  at  all,  determined 
together  by  the  extinction  coefficient  and  single-scatter  albedo,  while  the  asymmetry  parameter 
(operating  through  a  phase  function  approximation  or  embodied  in  the  truncated  Mueller  matrix 
used  here)  assumes  second-order  importance  by  determining  the  angular  distribution  of 
scattering  energy. 

A.4.  Methodology 

The  model  implemented  for  the  work  accomplished  in  Chapter  4  contains  five  separate 
modules: 

A.4.1.  Meteorological  Profile  Generation 

The  first  module  generates  the  temperature,  pressure,  and  air,  water  vapor,  cloud  liquid, 
rain,  snow,  and  graupel  density  for  each  1-D  level.  Drop  size  spectra  are  computed  using  the 
modified  MP  and  SS  formulas  for  user-defined  precipitation  size  (liquid  mass-equivalent  sizes 
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for  snow  and  graupel)  bins.  Two  versions  of  this  module  were  written.  The  first  replicates 
Petty's  (2001a)  1-d  parametric  cloud  model  for  a  regularly-spaced  vertical  grid  of  user-defined 
spacing.  This  version  generates  precipitation  rates  for  each  precipitation  species  at  each  level; 
the  MP  and  SS  distributions  are  used  directly  to  compute  size  spectra.  The  second  version  reads 
a  binary  MM5  model  output  file  and  produces  output  at  the  irregularly-spaced  heights  of  the 
sigma  levels.  The  MM5  uses  a  bulk  microphysics  scheme,  so  its  precipitation  output  is  a  mass 
density  for  rain,  snow,  and  graupel  at  each  grid  point.  These  are  converted  to  rain  rate,  using  fall 
speed  to  relate  the  two  quantities,  enabling  particle  size  spectra  computation  via  the  MP  and  SS 
distributions.  In  both  cases,  output  is  stored  in  text  format  for  use  by  succeeding  modules.  The 
profiles  in  Figures  4.1a,  .2a,  .3a,  .5a,  and  .6a  were  generated  by  this  module. 

A.4.2.  Mie  Property  Generation 

The  user  controls  the  instrument  characteristics-frequencies,  sidebands,  polarizations, 
and  scan  angles— used  by  this  module.  It  uses  these  to  compute  complex  refractive  indices  for 
rain,  snow  (modeled  as  a  solid  ice 
sphere),  and  each  unique 
composition  of  graupel  present 
(using  the  Bruggeman  mixing 
formula)  for  each  instrument 
frequency.  Next,  using  the  Mie 
solution  as  outlined  in  Bohren  and 
Huffman  (1983)  Chapter  4,  the 
module  computes  extinction  and 


Extinction  coefficient 


Single  Scatter  Albedo 


cloud  liquid  water 
hydrometeors 
■  total 


Figure  A.l:  Example  extinction  coefficient  (left)  and 
single-scatter  albedo  (right)  profiles  for  AMSU-A 
channel  8  for  the  profile  from  Chapter  4  case  9  (Figure 
4.3a).  Each  constituent's  contribution  to  total 
extinction  (bold  line)  is  included. 


135 


scattering  cross-sections  and  efficiencies,  and  Mueller  matrix  precursor  coefficients  used  in 
module  4  for  phase  function  generation  for  each  instrument  frequency  and  model  level.  Output 
arrays  are  stored  in  text  format  for  use  by  the  next  two  modules.  Figures  3.1  through  3.3  were 
created  using  output  from  this  module. 

A.4.3.  Extinction  Profile  Generation 


This  module  uses  the  extinction  and  scattering  coefficients  to  compute  the  extinction 
coefficient  and  single  scatter  albedo  for  each  instrument  frequency  and  model  level.  The 
Rosenkranz  (1998)  model  is  used  for  dry  air  and  water  vapor  absorption.  Cloud  water  absorption 

is  computed  using  the  Liebe  et  al.  (1991)  model.  Cloud 


ice  absorption  and  scattering  in  the  microwave  spectrum 
is  negligible  and  is  ignored.  At  each  vertical  level,  the 
total  extinction  and  scattering  coefficients  are  computed 
by  summing  each  size  bin's  extinction  or  scattering 


coefficient,  weighted  by  its  number  density.  After 
computing  a  layer's  total  absorption  (dry  air  +  water 
vapor  +  cloud  liquid  +  rain  +  snow  +  graupel)  and  total 
scattering  (rain  +  snow  +  graupel),  the  bulk  single-scatter 


Figure  A.2:  Example  phase 
function  for  AMSU-A  channel  2 
(top)  and  AMSU-B  channel  18 
(bottom)  for  a  model  level  with 
large  graupel  content  (and  hence  a 
large  asymmetry  parameter). 

Dots  are  for  a  photon  polarized 
nearly  vertically  (5°)  with  respect 
to  the  scattering  plane;  crosses  are 
for  a  photon  polarized  nearly 
horizontally  (85°). 


albedo  follows. 

AAA.  Phase  Function  Generation 

This  module  uses  stored  SI  and  S2  coefficients 
for  each  instrument  frequency  and  model  level  (see 
Bohren  and  Huffman  1983  Section  4.4.4)  from  the  Mie 
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Figure  A.3:  Flowchart  for  reverse  Monte  Carlo  radiative  transfer  model. 
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property  module.  In  a  manner  similar  to  that  used  to  compute  bulk  extinction  and  scattering 
efficiencies,  it  computes  effective  SI  and  S2  from  each  precipitation  size  bin's  number-density- 
weighted  value.  From  the  effective  S 1  and  S2,  the  upper-left  comer  Mueller  matrix  elements 
follow.  These  are  used,  as  described  in  Section  A.3,  with  an  incident  unit-intensity  photon 
rotated  through  all  rotation  angles  to  build  the  two-dimensional  scattering  direction  probability 
distribution.  For  each  scattering  angle  and  rotation  angle  bin,  the  scattered  intensity  (weighted 
by  the  solid  angle  subtended  by  the  bin's  angular  limits)  is  computed  and  normalized.  Each  solid 
angle  bin,  then,  contains  a  value  relating  the  scattered  intensity  in  that  bin  to  the  incident 
intensity,  or  equivalently,  the  probability  of  the  incident  photon  scattering  into  that  bin.  This 
scattering  direction  distribution  array  is  stored  as  a  text  file  for  each  instrument  channel  and 
model  level.  Figure  A.2  shows  an  example  scattering  probability  as  a  function  of  scattering 
angle,  6,  for  rotation  angles  of  5°  and  85°. 

A.4.5.  Reverse  Monte  Carlo  Radiative  Transfer 

The  final  module  reads  the  temperature,  extinction  efficiency,  single-scatter  albedo,  and 
scattering  direction  probability  density,  then  performs  the  RMC  radiative  transfer  for  a  user- 
specified  subset  of  instrument  channels  and  scan  angles,  and  number  of  photons.  Results  are 
depicted  graphically  and  output  in  text  format.  Figure  A.3  is  a  flowchart  outlining  the  RMC 
code. 

A.5.  Future  Work 

As  discussed  in  Section  A.3,  several  assumptions  are  inherent  in  the  design  of  the  RMC 
model.  Several  of  these— the  assumed  particle  size  distributions,  the  spherical  geometry  inherent 
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to  Mie  theory,  the  graupel  composition,  and  the  dielectric  mixing  formula— are  well  beyond  the 
scope  of  this  work  and  are  limitations  affecting  nearly  all  radiative  transfer  models. 

One  assumption,  though,  can  be  improved  by  further  work:  the  treatment  of  polarization. 
While  a  full  implementation  of  the  Mueller  matrix  would  be  computationally  prohibitive  for  a 
large  3-dimensional  domain,  a  parameterization  of  individual  matrix  elements  in  terms  of 
precipitation  density  may  be  possible.  This  would  enable  construction  of  the  matrix  elements  via 
a  look-up  table  instead  of  Mie  calculations  at  each  grid  point.  The  photon's  full  polarization 
state,  expressed  via  the  Stokes  vector,  would  then  be  tracked  through  the  photon's  path  to  its 
emission  point.  If  the  photon  encounters  the  surface,  the  Fresnell  scattering  matrix  would  be 
employed.  The  underlying  assumption  of  reciprocity  remains  in  place,  allowing  the  Mueller 
matrix  to  be  used  for  scattering  interchangeably  in  the  forward  and  backward  directions. 

The  results  from  Chapter  4  demonstrate  the  validity  and  utility  of  this  version  of  RMC 
radiative  transfer  model,  and  the  author  looks  forward  to  further  developing  it  and  employing  it 
for  additional  applications  in  the  future. 
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